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PREFACE 


The Digital Computer Laboratory of the University of 
Illinois has a three-fold program which is concerned with re¬ 
search in the field of digital computers, with maintenance for 
University use of a high-speed digital computer, and with teaching 
of design and use of digital computers. The preparation of pro¬ 
grams for the Illiac is the responsibility of the person who wishes 
his problem solved. But since each prospective user cannot always 
attend the courses on programming offered by the University, it 
has seemed desirable to help make it possible for him to learn 
the elements of Illiac use by himself. 

The result has been this book, first issued in March 195^-, 
which makes use of two years of experience in operating the Illiac 
for University research and several years of experience in the 
teaching of programming. Many of the chapters were used in pre¬ 
liminary form as supplementary material in Mathematics 3^5/ the 
basic course on programming, and the book is now used as standard 
course material. 

Like most of the work of the Digital Computer Laboratory, 
this is the result of a group effort by many different people. The 
chapters were written by S. Gill, R. E. Meagher, D. E. Muller, 

J. P. Nash, J. E. Robertson, T. Shapin and D.J. Wheeler. Other 
members of the Laboratory staff read the material and offered 
valuable criticism. The typewritten copy was prepared by Helen B. 
Clark and Muriel P. Fetzner; the drawings were made by George 
Ehrlich and Charles Breed. 
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CHAPTER 1 


INTRODUCTION 


In 1948 the University of Illinois began to look into the 
possibility of acquiring an automatic digital computer. When attempts 
to buy one or to have one built failed, the University decided to 
build one of its own, and the Digital Computer Laboratory was organized 
in February 194 9 . Shortly thereafter an agreement was made to build 
a second computer for the Army to be used at the Aberdeen Proving 
Ground. 

It was decided to build the two machines following the de¬ 
sign of a computer then being designed and built at the Institute for 
Advanced Study. While many modifications were later made, both 
machines built by the University owe much to the early designs ob¬ 
tained from Princeton. 

The two machines were built more or less together so that 
advantage could be taken of the savings obtainable by paralleling 
the design and construction work, but work on the Army machine, called 
the ORDVAC, was emphasized, and it was finished first. It passed its 
acceptance tests in February 1952 at Aberdeen and it has been in use 
there since that time. 

Work on the University computer (later named the Illiac) was 
completed in September 1952 and the computer was first made available 
for University use when classes began on Monday, September 22, 1952. 
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The Illiac is an automatic electronic digital computer. 

It is digital because it handles numbers as sets of digits 
which have discrete values, rather than as scale readings or measure¬ 
ments, which are continuously variable. Apparatus for handling digits 
is more complicated than that for handling continuous quantities, but 
it is capable of giving unlimited accuracy by using suitable numbers 
of digits. 

The Illiac is electronic. In the last ten years electronic 
circuits for storing, transmitting, adding, subtracting, multiplying 
and dividing numbers in digital form at extremely high speeds have been 
devised. The actual addition of two numbers in the Illiac t ak es only 
about 75 microseconds. 

Such speed is useless unless the machine can be made to go 
ahead on its own with may thousands of operations, without human inter¬ 
vention. The Illiac is therefore automatic, in the sense that it can 
be given orders telling it how to proceed, and will then act on these 
orders automatically. 

In common with many other computers of a similar type, the 
Illiac contains the following five essential features: 

(1) An arithmetic unit. 

( 2 ) A memory or store. 

( 3 ) Devices for the input and output of 
information (e.g. numbers) to and 
from the machine. 

(4) Means for the transfer of information 
between the various parts of the machine. 
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( 5 ) Means for the automatic control of 
the whole machine. 

The arithmetic unit carries out the individual arithmetical operations 
that mak e up every computation; it can be thought of as the electronic 
equivalent of a desk calculating machine. It is described in detail 
in Chapter 2. 

The memory is needed because, in any lengthy calculation, 
numbers produced at early stages of the calculation are frequently 
required to be used at later stages; they must therefore be recorded 
or '’remembered". The memory is capable of recording 102k numbers.'*' 

These can be recorded (i.e. transferred to the memory from the arith¬ 
metic unit) individually, as directed by the computer's control de¬ 
vice, and recalled again individually in a similar way. The memory 
may be thought of as 102k little boxes or locations, each accomodating 
one number, and labelled with the numbers 0 through 1023. The label 
of a location is called its address . A number in the memory is 
identified by the address of the location containing it. 

Information enters and leaves the Illiac coded in the form 
of a pattern of holes in punched paper tape; there is a tape reader 
for input, and an automatic punch for output. There is also a tele¬ 
typewriter which can be used to provide output from the machine 
directly in printed form. Several machines are available for pre¬ 
paring punched tape, copying it, comparing it, and producing printed 
versions of the information on it. 

1 

This is the high-speed memory. The magnetic drum furnishes additional 
storage. 
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The problem of controlling the whole computer has been 
solved by stipulating that every individual operation that occurs within 


the machine must be one of a certain set of specified permissible 
operations, and that no two such operations can occur simultaneously. 
Thus the design problem was reduced to that of engineering the various 
permissible operations and arranging for them to be executed in any 
desired sequence. It is up to the user of the Illiac to specify the 
sequence of operations or program , which the Illiac must execute to 
carry out his calculation. 

Each permissible operation can be specified in a concise 
coded form called an order . The correspondence between the set of 
permissible operations and the set of orders which specify them is 
called the order code of the Illiac. It is given in detail in Chapter 
3» A coded problem is called a program or routine . 

The machine’s control unit has the task of accepting orders 
one by one, and of causing the machine to carry out the operations 
specified according to the order code. If each order were taken by 
the control unit directly from a punched tape, then to make full use 
of the speed of the rest of the machine the tape would have to pass 
through the tape reader at about 200 miles per hour. Instead, the 
orders are recorded in the memory along with the numbers, so that the 
control unit merely has to take Its orders from the memory, which it 
can do electronically at high speed. This is made possible by coding 
each order to look like a number. To be more precise, orders are 


1-4 



stored in pairs, one pair to a memory location. The information 
contained in one memory location is often called a word, meaning 
either a number or an order pair. Of course, the more orders there 
are in the memory, the less room there is for numbers. Both orders 
and numbers are fed into the machine initially on punched tape. 

Normally orders are obeyed by the control unit in the 
sequence in which they are stored in the memory, e.g.: 

Left-hand order in location 6, 

Bight-hand order in location 6, 

Left-hand order in location J, 

Right-hand order in location 7 , 

Left-hand order in location 8, etc. 

Sometimes, however, this sequence is broken and the control unit 
starts over at some new position in the memory; this is called a 
transfer of control . There are special orders which cause this. 

There is also special provision for making a transfer of control 
depend on the value of some number obtained by the machine during 
the calculation. Thus the machine can be made to "choose" one of 
two or more alternative courses of action according to the way 
things happen to work out. 

If control is transferred a few locations back in the memory, 
the machine will repeat the operations specified by the intervening 
orders. It is possible to cause this repetition to occur any number 
of times, leading to a cyclic behavior of the machine. Practically 
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every calculation which, the machine performs contains several such 
cycles, often one inside another. In this way it often happens that 
the same order gets carried out many thousands of times, so that a 
few orders suffice to keep the machine busy for several minutes.. If 
each order in the memory were to be carried out once only, the Illiac 
would get through them all in a quarter of a second (even if the memory 
contained nothing but orders). In practice, calculations vary in 
duration from a minute to a few hours. 

The occurrence of cycles is one of the things that compli¬ 
cates the programming of a calculation. Another is the fact that, 
since orders are stored in the memory in the same form as numbers, 
they can be operated on and altered during the course of a calculation 
(at the behest of other orders) just as if they were numbers. All this 
makes possible some most interesting calculations; it can also 
programming difficult. 

Fortunately a coder can often, as described in Chapter 4, 
make use of bits of programming done by other people. Thus a typical 
program consists of a number of groups of orders, some written by the 
coder, others already available. The latter will be available in 
punched tape form, and can be copied mechanically onto the program 
tape along with the new orders. Tape preparation is described in 
Chapter 9° 

When the whole tape for a particular program has been prepared 
it can be placed in the tape reader of the Illiac. The Illiac reads 
the tape, forms the orders and numbers punched on it and stores them 
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In the memory. When the program is in the memory, the machine begins 
to execute the orders, continuing until it comes to some particular 
order which causes it to stop. If the programming is correct, this 
is the end of the calculation. If there is a mistake in the programming 
various things may happen; remedies are discussed in Chapter 8. 

Somewhere in the program will be some orders which cause 
the machine to punch some output tape. This carries the results of 
the calculation. The program may also contain orders causing the 
machine to read more input tape, carrying data for the calculation. 

The reading of most of the program tape is accomplished by 
the Illiac executing a particular set of orders called the Decimal 
Order Input (See Chapter 5) which is always punched at the beginning 
of every program tape and hence read into the machine before the rest 
of the tape. The Decimal Order Input not only assembles the program 
Inside the machine; it also makes certain modifications and conver¬ 
sions, so that the way in which orders are represented when punched 
is somewhat different from their final form in the memory. The 
object is to make programming easier. It is important to remember 
that the written form of an order and the form which it assumes in 
the memory are not the same thing. The relationship between the two 
is determined by the Decimal Order Input. 

Remaining chapters of this manual are devoted to: The 
arrangement of calculations so that all the numbers encountered 
are the right size (Chapter 6); ways of programming certain types of 
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simple tasks (Chapter 7 ); how to estimate the duration of a calculation 
(Chapter 10 ); an example of a typical, program (Chapter ll); and how the 
program library is organized (Chapter 13), Finally, concise descriptions 
of the principal contents of this collection are given,, 
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CHAPTER 2 


THE ARITHMETIC OF THE ILLIAC 

The construction of a computer involves a compromise 
between engineering economies on the one hand and ease and flexi¬ 
bility of use on the other. As a result, the details of the 
operations of arithmetic are often dictated by engineering design 
considerations. The following paragraphs describe the peculiarities 
of the representation of numbers and of the operations of arithmetic 
in the Illiac. 

2.1 REPRESENTATION OF NUMBERS. The simplest device for 
storage of numerical information is an electronic element having two 
stable states in which one state is called "zero" and the other state 
"one". One such element is capable of holding one binary digit. 
Registers composed of forty such elements are provided in the Illiac 
for representation of forty binary digit numbers. The Illiac is 
described as a binary parallel digital computer having a precision 
o ,f forty binary digits — roughly equivalent to twelve decima l 
digits. 

The Illiac is a fixed-point computers the location of 

y _ . _——— -— —-- “ " 

the fixed binary point between the first and second digits is in- 

U. -—-----—— 

dicated in Figure 2.1. 
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^/location of binary point 

3£ 


i 


v. 


sign digit 


-v- 

39 non-sign digits 


Figure 2.1 

The Forty Digits of a Word 


.. The l eftmost dig it of a number is the .sign digit. If the sign 
^ digit is zero, the number is pos itive or zero; if the sign digit is 
t r- one > the number is negative . For example, + 7/8 is represented in 
the Illiac as 


0.11100 ... 00 = + 1/2 + 1/4 + 1/8 = + 7/8 

Negative Numbers« Negative numbers are represented in 
Illiac as complements with respect to 2• The process of comple¬ 
mentation is carried out by forming the digitwise complement (re¬ 
placing ones by zeros and zeros by ones) and then adding a unit in 
the least significant (thirty-ninth) non-sign digit. As an e xam ple, 
- 7/8 is formed from + 7/8 by complementation as follows: 


0.11100 ... 00 
1.00011 ... 11 

+ 0.00000 ... 01 

1.00100 ... 00 


+ 7/8 

ones replaced by zeros and zeros by ones 
addition of 2 - ^ 

machine representation of - 7/8 


It is sometimes essential to distinguish between the 


2-2 


machine representation of a number and its arithmetic value. Suppose 
we have a number whose arithmetic value is x and whose machine re¬ 
presentation has the sign digit x^ and non-sign digits x i (i = 1,2 ,...,39). 
If x is positive} i.e., if x Q = 0, then 

39 -i 

x = Z 2 x . (2.1) 

i=l 1 


If x is negative (x Q = l), then the relationship is 


39 


,-i 


LiU-O 


CA-' 




x = -1 + Z 2 * x. .1 . | 

i-1 1 \ t 11 


( 2 . 2 ) 


Equations 2*1 and 2.2 may "be combined in the one equation 


X = -x 0 + 


39 -i 
Z 2 x. . 
i=l 1 


(2.3) 


If equations 2.1 through 2.3 are to hold, x must be restricted to 
the range -1 < x < 1. The programmer must keep in mind that the 
Illiac arithmetic unit is designed for numbers in the range -1 < x < 1. 
When an operation is performed which yields a number outside this 
range, an unwanted number within the range results. This effect is 
called overflow and, although certain operations are unaffected by 
overflow, computational errors usually result. 

2 .2 FUNDAMENTAL OPERATIONS OF THE ILLIAC ARITHMETIC UNIT. 
The following five operations are fundamental in the arithmetic unit 
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of the Illiac: 


c onrp 1 ementa t ion, 
addition, 

left shift, or multiplication by two. 


1 . 

2 . 

3 . 

^• right shift, or division by two, 

5 * clearing to zero. 

Complementation is executed in the Illiac in the manner 
described in the previous section by a circuit called a complement 
gate. The complement gate is actuated by a signal from the control 
and is capable of supplying either the number unchanged, or its 
digitwise complement. Addition is carried out by the adder, which 
is capable of forming the sum of two 4 o binary digit addends. 

A-J- e ft shjd^t is a displacement of the binary digits one 
digital position to the left and corresponds to a multiplication by 
two. ,There will be an overflow if the nu mber to be shifted is 
+ l/2 or if it has an absolute value greater than one-half. Th* 3 

----*-— 4 *- 

right shift is a displacement of the binary digits one digital 

position to the right and corresponds to a division by two. For 
example, 

0.0111 shifted left is 0.1110; 7/l6 x 2 = 7/8 

0.1110 shifted right is 0.0111; 7/8 -f 2 = 7/16 

1.0010 shifted right is 1.1001; -7/8 x 2 = -7/16 

1.0100 shifted left is 0.1000; -3/4 x 2 = l/2 because 

of overflow. 
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It should be noted that the sign digit is propagated when the right 
shift is executed. The Illiac arithmetic unit has two shifting 
registers, each capable of executing both left and right shifts. 

Clearing to zero involves setting all digits of a number 
to zero; the corresponding arithmetic value is zero. 

2„3 THE ILLIAC ARITHMETIC UNIT . The structure of the 
Illiac arithmetic unit is shown in Figure 2 . 2 . The arithmetic unit 
is composed of two shifting registers, the accumulator A and the 
quotient register Q, and one non-shifting register, the number 
register R^. Also required are the complement gate and the adder. 

The A and Q registers are the only registers to which the programmer 

O 

has direct access; the number register R” 1 is used to hold temporarily 
the numbers brought from the memory for arithmetic operations. It 
is essential that the programmer be familiar with the roles played 
by the A and Q registers in the operations of arithmetic; many pro¬ 
gramming errors arise from placing operands in, or removing results 
from, the wrong register. The functions of the registers during the 
operations of arithmetic are described in the sections which follow. 

2.4 ADDITIONS (order Ty pe L ). Before an addition in- 
structi on begins, the augen^ . lies in the accumulator register A. 

During execution of the addition instruction, the addend is trans- 

3 

ferred from a specified memory location to the number register R • 

The digits of R^ are then sensed through the complement gate unchanged, 
so that the addend forms one of the inputs to the adder. The augend 
in A is the second adder input. The adder forms the sum which is 
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transferred to the accumulator A, replacing the augend. The quotient 
...register Q is undisturbed by the addition instruction. 



Figure 2.2 

The Illiac Arithmetic Unit 

Two variants of addition are "hold add” and "clear add." 
?he "hold add" instruction leav es the augend in A undisturbed 
until the sum is f ormed by the adder. The "clear add" instruction 
clears A to zero initially, thus setting the augend to zero. The 

' Iclear add" instruction i s thus a transfer order which moves a 

^ ” “ _ _ __ _______________ 

number from 'a specified memory location to the accumulator A. 

2.5 SUBTRACTIONS (Order Type L). Subtraction in 
the Illiac arithmetic unit is performed by adding the complement 
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of the s ubtrahend to the minuend . Before the subtract instruction 
begins, the min uend lies in the accumulator A , and is used as one 
of the two adder inputs. The subtrahend is brought from a specified 
memory location to R 3 ; its complement is formed by the complement 
gate and is used as the second adder input. The adder thus forms 
the difference by forming the sum of the minuend and the complement 
of the subtrahend. The adder output is then transferred to the 
accumulator A, replacing the minuend. 

Either the "hold subtract" instruction or the "clear 
subtract" instruction can be used by the programmer. For the former ' 
instruction, the result of a previous operation is left in the accumu¬ 
lator A as minuend; for the "clear subtract," the minuend in A is 
set to zero, so that A contains th e negative of the number in a 
specified memory loc ation when the operation has ended . 

2.6 ABSOLUTE VALUE ADDITION AMD SUBTRACTION ( Order Type L) . 
It is possible to form the absolute value of the addend (or sub¬ 
trahend) after it has been transferred from a specified memory loca¬ 
tion to R 3 and before it is added to the augend (or subtracted from 
the minuend) in the accumulator. Ordinarily the setting of th e com¬ 
plement gate depends upon whether the instruction is an addition o r 

3 

a subtraction; for forming the absolute value of the numbe r in R , 
it is necessary to sense the sign digit of R 3 as wel l. For example, 
if the addition of the absolute value of a negative addend is requir ed, 
the sign digit of the addend in R 3 is sensed indicating that a 
complementation is necessary. The add instruction o rdinarily doe s 
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HQ^j^g^jjre ^mpleinentation; the net effect of the sensing of the 
add instruction and the R 3 sign digit is that the complement gate is 
/''get to form the complement of th e addend. 

2.7 INCREMENT ADD ORDERS (Order Type p). It was noted 
in section 2.1 that the operation of complementation is performed 
by first forming the digitwise complement of the number held in R 3 
and by then adding a unit in the least significant digital position 
of the adder. ^For the ordinary additio n and subtraction orders 
.-^° r ^ er t yP e k), the least significant digit insertion occurs only 
^when the complement gate is set to form the digitwise complement. 

In the Illiac,^ special increment add orders (type F) are provided. 
F or these orders, the __rel ationshi p between the setting of the com- 
pl ament gate and the in sertion of the least significant digit is 

reve rsed ._ Thus,, the "clear increment .add" instruction s ets the 

augend initially in A to zero and adds the addend from R 3 to 2~ 39 
^and places the result in A. Conversely, the "clear increment sub¬ 
tract" instruction sets the minuend in A to zero, replaces ones 
of the s ubtrahend by zeros and replaces ze ros of the subtrahend 
by ones, and tran sfers the digitwise complement of the subtrahend 
thus formed from the adder to A . Detailed descriptions of further 
orders of this type are given on page J-PR . 

2.8 ADD FROM Q (Order Type S) AND INCREMENT ADD FROM Q 
(Order Type K) INSTRUCTIONS For the add from Q and increment add 
from Q instructions, the addend or subtrahend is tranferred to the 
number register R 3 from the Q register, rather than from a 
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location. Otherwise, the K and S order types 


specified memory 
are the same as the corresponding F and L orde r types describe d 

ip sections 2.4 through 2.7. 

2.9- THE SHIFT ORDERS (Order Types 0 . l). Since the 

right and left shift operations are fundamental in the Illiac 

arithmetic unit, specific shift instructions are provided for 

the programmer. All digits in "both A and Q, except the sign digit 

a of Q# are shifted. A left shift (order type 0 ) of one digital 
H o 

position replaces 

V V V* • * ' a 37 ' a 38 ' a 39 in A 

and 

V V q 2 * ’ ’ * q 37' q 38 ' q 39 in Q 
by 

a l' & 2 . J a 3 > * * * ' a 38 ’ a 39 ' q l in A 

and 

V V q 3' • • • ' q 38 r q 39’ 0 in 

The right shift (order type l) of one digital position replaces 

V V a 2' * • • ' a 37’ a 38 " a 39 in A 

and 
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Iq* • • • t I37J ^38* q 39 ^ 

by 

a 0' a 0’ V ' * ‘ ' a 36' a 37' a 38 in A 

and 

q 0 > a 39/ q l' ’ * ' * q 36' q 37 # q 38 in 

From the behavior of the left and right shift we see 
that we may consider the shifting to take place in a single register 
AQ of 79 digits consisting of A followed by Q with q Q deleted. 

Thus the left shift replaces 

a 0> a p> a 2 > • • • > a 37 ? a 38^ a 39 ; * ■ * > q 37^ ^38* q 39 

by 

a l' a 2' a 3' • ' • ' a 38 # a 39' q l ; V q 3' * * * ' q 38' q 39' ° 

While the right shift replaces 

V V V * • * ' a 37' a 38 ' a 39 ; V V * * * ' q 37' q 3S' q 39 

by 

a 0' a 0 ’ a i* • • * * a ^ 6 > a 37 # a 38 J a 39' q i JI * ’ * * ^37* ^38' 

The number 1 n of shifts can be specified in the address 
digits of the shift order by the programmer. A left shift of n 


2-10 



digital positions replaces AQ by 2 n AQ. Similarly, a right shift of 
n digital positions replaces AQ by 2 n AQ$ n lies in the range 
1 < n < 63 . 

2.10 MULTIPLICATION (Order Type 7). Initially, the 
multipl ier lias-in the Quotient register Q. At the beginning of the 
multiplication instruction, the multiplicand is transferred from a 
location specified by the address of the multiply order into the 
number register R , where it remains throughout the multiplicati on. 

The multiplication then consists of a sequence of additions and shifts 
More precisely, a multiplier digit in q^ is sensed. There are two 
cases, (a) If q XQ = 1, the multiplicand is added to the partial 
product in A. A right shift in AQ follows which halves the number 
in A, moves a new multiplier digit into q^ and transfers 
a product digit from a to q^ (b ) If q ^ = 0, only the right 
shift occurs, which transfers a, n to q. and transfers the next most 
significant digit of the multiplier into q^. After thirty-ni ne 
right shifts have occurred, the sign digit of the multiplier is s ensed 
If the multiplier is positive, the multiplication is c omplete, with 
a double precision product (sign digit and seventy-eight non-sign 
digits) in AQ. For a negative multipli er x, the process of thirty- 
nine shifts and conditional additions constitutes a multiplication 
of the multiplicand y by 1 + x (Equation 2 . 3 ) to form a product 
y(l + x) = y + xy. In this case the Illiac automatically subtr acts 
the multiplicand y to produce the correct product xy. 
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A number of variations of the basic multiplication pro¬ 
cedure described in the previous paragraph are possible. The initial 
contents of A may be either 0, l/2, or some quantity previously 
calculated. If we designate t he initial contents of A by a, the 
final double precision product in AQ is xy + 2~^ a , where g. and y 
are multiplier and multiplicand , respectively. For a/= 0. the 
i nstruction specified is "clear multip ly" and the result in AQ is 

the exact 78 digit signed product xy . For a = l/2, the instruction 

* __ 

specified is "roundoff multiply" ancf ~bhe result, in A is a rounded 
/$9 digit signed produ ct. If a is arbitrary, the instruction 
specif ied is "hold multiply", and the result in AQ is xy + 2~^a: 
t ^e quantity a is thus added to the least significant part of the 
product in Q . I n. all multiplication instructions the s ign digit 
q n of the quotient register is set to zer o. tPVi- 

For each of the three types of multiplication instructions 



described in the previous paragraphs, four additional variations 
can be specified by the programmer. Either N(n), -N(n), |N(n)|, 
or - |N(n)| can be used as the multiplicand, where w(n) is the 
number transferred from memory location n to the number register 


2.11 DIVISION (Order Type 6 ). Initially, the doubl e 
precis ion dividend lies in AQ. The divisor is tran sferred at the 
beginning of the division instruction from a memory location 


specified by the address of the divide order to the number register 



r 3 . For .positive divisor and dividend, the process is analogous 



tp elementary long divis ion. The d ivisor is subtracted f rom a 

■ partial remainder in A and the sign of the differ ence (in the adder) 

* ” ' ~ ~ _ . . *7 

is^ sensed. >-If the difference is negative, 0 is inserted in * 

as,quotient digit and AQ i s doubled to form a new p a rtial remaind er. 
I f the difference is positive, 1 is inserted in q^ as quotie nt 
d igit and the difference in the adder is placed in A, and AQ is 
doubled to form a new partial remainder. At each doubling of AQ, 
is shifted into q^ as well as into a^ . Thus after 39 step s 
q n has the sign of the quotien t. The dividend is used as the initial 
partial remainder; after 39 quotient digits have been generated, the 
process is complete. A similar procedure is employed if divisor 
or dividend is negative or if both are negative. 

The Illiac division has the following properties: 

(1) A rounded quotient Q is always generated. 

The roundoff is achieved by setting q_ Q = 1 
—-.-- _Jlz _ 

in all cases. 

(2) The thirty-ninth partial remainder is left 
in A and is called the residue r . The true 
remainder R co rresponding to the rounded quotient 
is related to r approximately as f ollows: 

R - r 4- (2q 0 - l) y, , (2.4) 

w here y is the diviso r and q n is the sign 
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digit of the quotient - Thus if the quotient 


is positiv e 

„ .= L - y, ( 2 . 5 ) 

and if the quotient is ne gative 

r / . (2.6) 

/The equations given abov e for the true remainder 

,are valid if the absolute value of the quo tient is 

less than one, and they yield results in error by 
-■50 

j}ot more tha n 2 . 

The sign digit of the quotient replaces the least 
gj?gnifleant ( 78 th) digit of the double -precision 
. dividend. One effect is that the least significant 
digit of the residue (a^) is 'the same as the sign 
digit of the quotient q Q . 

If we have a priori knowledge of the true value q 
of the quotient (such as, for example, a division 
of qy by y), the relationship between the machine 
quotient Q and the true quotient q is 

Q = q + 2’ 39 (1 - q 39 ) (l - 2 y 0 ) ( 2 . 7 ) 

where q is the least significant digit of q and 
y Q is the sign digit of the divisor y.. Equation 
2.7 is essentially a description of the division 
roundoff. If q^ = 1, then Q = q and the machine 



qua ti ent is the true quotient . If q^ = _Q__ 

r .the nature of* the roundoff depends upo n the _ 

sign of the divisor. Suppose, for example, that \ 7 
__ \ / 

q is 0 . 101 ; then Q is either 0 . 101000...001 or j 
0 . 100111 .. .111 depending upon whether the divisor/ 
was positive or negative, respectively. ^J) 

2.12 PRECISE CALCULATION OF THE DIVISION REMAINDER. We 
define th e remainder R in relation to a . _ quQ . t i eht ... Q .i » a divisor 
and a dividend & hy the equation 

Qy + 2 -59 R = d, 
or 

R = 2 ^ 9 (d - Q-v)._ (2.8) 

However, the exact relationship involving the Illiac residue r and 
including the replacement of the least significant dividend digit 
d^g by the quotient digit q Q is 

Qy + 2 -59 [r + (2q Q - l)y] = d + 2 - <* 73 )* ( 2 *9) 

Solving for R, we have 

R = 2 59 (d - Qy) = r + (2q Q - l) y + 2 59 (d^ - q^)- (2.10) 

Equation 2.10 gives the exact expression for the remainder R which 

corresponds to the machine quotient Q. 

If the remainder R corresponding to the true quotient q 

is desired, then R is defined as 
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(2.11) 


R = 2^ (d - qy). 

The value of R is found by substituting equation 2.7 in equation 
2*9 and solving for R* which yields 

5 - r + [ 2qo - i}9 - (l . 4j9 ) 2yo ] y + 2-39^ . qg y (a.ja) 

2 -13 THE DIVISION HANGUP. SPECIAL CASES OF DIVISION. 
Cir cuits are incorporated in the Illiac for stopping the Ill iac if 
the quotient re sulting from a division exceeds one in absolut e value. 

Thensign digi t of the quotient is predictable from the s igns of divisor 
ang^ dividend. The sign digit of the quotient is calcu lated in the 
ty compa ring the dividend and divisor arithmet ically. Thus, 
by sensing the sign digits of quotient, divisor, and dividend, it 

% 

is possible to detect the fact that the quotient exceeds one in 
absolute value, and stop the Illiac. 

cy Th e J- equations derived in Sections 2.10 and 2.11 are valid 
22^- if th j L abs olute valu e o f the dividend is less than the ab solute 
valu e of the divisor . Wfreh absolute va lues of dividend and divisor 
.are equal, the Illiac generates a quo tient -1 + 2 ~^ if the dividend 
l-I s negative and the divisor is posit ive. The quotient is +1 - 2~ 39 
^ dividend and, divisor are both negative. I f the dividend is 
pos itive and equal to the absolute value of the divis or, the Illiac 
. wii; stop . If the _divisor is -1, the Illiac generates aj piotient 
w hich is the digitvise complement of the dividend, except for the 
quotient roundoff. 
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2.14 MEMORY TO Q (Order Type 5 ) AMD STORE (Order Type 4 ) 


INSTR UCTIONS . Instructions are provided in the Illiac for trans¬ 
ferring a number from a specified memory location to the Q, register, 
and for transferring a number in A to a specified memory location. 

The former i nstruction can be u sed to transfer the tttttI tipi i er - 

v| 

to a before a multiplicat ions the latter instruction is used to- 

transfer a result from A to the memor y. If a result li es in. Q(e.g., 
a division quotient) it can be transferred to the memory by using tw o 
orders. S5 and 40 which transfer ,Q to A a nd then A to the memory . 


2.15 ILLUSTRATIVE EXAMPLES. 

A. The Leapfrog I Division Test. In the Leapfrog 
I division test the product xy is formed and is then divided by y. 

The sum of the quotient in Q and residue in A is formed and stored.^ 
^he sum of quotient and residue is then calculated independently 
and compare d with the sum previously stored . 

The value of the machine quotient Q in the quotient regis- 
ter after division of xy by y is given by equation 2*7.> with q = x. 
Thus 

Q = x + 2“ 59 (1 - x ?9 ) (1 - 2y 0 ). 

The value of the residue r left in the accumulator is found from 
equation 2.12 by setting R = 0, x = q, and solving for r: 
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r = [x 39 + (1 - x^) 2y 0 - 2 x q ] y + 2" 39 (x Q - y^). 

The sum Q + r, after rearranging terms, is 

Q + r = x - 2~ 39 x ?9 y^ - x Q (y - 2~ 39 ) + (l - x Q ) y + 

(1 - * 39 ) (2y Q - l) (y - 2~ 39 ). 

The independent calculation of Q + r thus consists of forming 

-39 _zq 

x ” 2 x 39 y 39 and either adding y or (-y +2 ) depending upon 

whether x is positive or negative; and finally, if x,_ is 0. adding 

39 

or subtracting (y - 2 ^) depending upon the sign of y. 

B. The "Double Precision" Division. It is 
sometimes convenient to consider a single precision divisor y as 
exact and form a "double precision" quotient s + 2~ 39 t utilizing 
the double precision dividend d originally in AQ. The procedure 
used i6 as follows: 

1. Form d/y, yielding s + 2~ 39 in Q and r + 2" 39 s Q in A. 

2. Shift right, forming r/2 in A and s/2 in Q. 

. 7 

3« Form r/2 + y leaving ^T*in Q and a residue u in A. 

4. Assemble s + 2 ^ 9 t by setting t Q (the sign digit of T) 
to zero, inserting s/2 in A and shifting left once. 

-39 

The precision of s + 2 ■'■'t can then be calculated as follows: 

Step 1 yields s and r such that (by equation 2*9) 

(s + 2 -39 )y + 2 -39 [ r + 2 _39 s Q + (2s Q - 1 )y J = d + 2“ ?8 (s Q - d^). 
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or 


(2.13) 


sy + 2~ 59 (r + 2s Q y) = d - 2~ 78 d...g. 

Step 3 yields 

Ty + 2' 59 [ u + (2t 0 - l)y ] = r/2 +• 2" 78 t Q . # (2.14) 

The following substitutions are made; 

(a) s q “ since the.sign of r/2 is the same as the sign of 

d. 

(b) From step 4, t = 2(T + ) or T = l/2t - s^. 


Substituting in equation 2.14 and solving for r, we have 
0 ,J L ” ' s ~~0 


r = (t - 2s n )y + 2 -58 [ u + (2s n - l) y ] - 2 77 s Q . 


Substitution of the value for r in equation 2.15 yields 
sy + 2 -59 [ (t - 2s q )y + 2s^y ] + 2 -77 [ u + (2s - 1 )y ] 

- 2 " 116 s o ■ d - 2 ‘ 78 d 78 - 


We thus have 

(s + 2 " t)y t 2 ^ [ 2u + 2 (2s^ - l) y + - 2 s Q ] = d. 

It can be shown that u + (2s^ - l)y < 1 so that the quantity 

within square brackets is less than 5 in absolute value. We. conclude 

-59 -76 

that the quotient s + 2 t is in error by not more than 2 

The program for forming a "double precision n quotient 
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.from a double precision dividend and an exact single precisio n 
divisor is given in wor ds 57 to 66 of li brary routine A 4 entitled, 

J ’l.7 Precision Fl oating Decimal ". A major difficulty encountered 
is the formation of-r/2 from r (Step 2 above), for r may be as large 
as 2y and may therefore exceed range- It is therefore necessary 
to set the sign digit of r/2 to that of the original dividend d. 

2‘l6 INTEGER OPERATIONS• It is sometimes desirable to 

use integers for computations in the Illiac . Suppose we have an 

integer a stored in the memory or arithmetic unit. In terms of the 

formulation of previous sections, we would store 2~^a, where a 

lies in the range -2 59 < a < 2 59 . If we wish to add or subtract 

two integers a and b, no difficulty is encountered, for 2 ^a + 2~^b = 

2 (a + b ) indicating that the correct sum or difference lies in 

A after the instrucion is performed. Multiplication of two integers 

a, b yields (2 ^a) (2 ^b) = 2 ^ ab. The product ab lies in AQ 

TS 

and is in the range -2 < ab < 2 . If the programmer scales 

all quantities so that the product remains in the range -2'^ < ab < 2^, 
then the 40 digit signed product can be transferred to A by a left 
shift of 39 digital positions. It should be noted that the sign 
digit of Q is set to zero during the multiplication so that for 
positive products, N(Q) = ab, if 0 < ab < 2^. 

Division of integers presents certain difficulties . An 
example is given here of a method of dividing a positive dividend 
a by a positive divisor b to yield a quotient f and remainder g. 
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The steps are as follows; 

(1) Place the dividend 2 h in AQ. ) < a < 2^1 

. * - " ’ 

(2) Shift left one digital position, leaving 2a in AQ, with 



(3) Divide by 2~^ 9 b, leaving (2f + 1) 2 ^ in Q and (2g) 

2 -59 in A. 0 < b < 2 39 , 

-39 

(4) Shift right one digital position, leaving 2 f in Q 

“39 

and 2 g in A . 

It can he proved that bf + g = a by substitution of the 
appropriate quantities in equation 2.9* as follows; 

(2f + 1) 2 -59 (2 _59 b) + 2' 59 [ 2 _59 (2g) - 2“ 59 b ] = 2 _T8 (2a) 

which yields 

2 -78 (2bf + b + 2g - b) = 2 - ^ 8 (2a) or bf + g = a. 

~zQ ~z Q 

The ranges of f and g are 0 < f <2? and 0 < g < 2^ . 

2.17 SUMMARY. In the Illiac arithmetic unit are two 
registers, A and Q, which are directly accessible to the programmer. 

A single arithmetic order of a program utilizes the initial numerical 
operands in A, in Q, and in a specified memory location, and trans¬ 
forms these quantities to produce desired results which are left in 

~"i 

the registers. The programmer must know where the operands are 
initially located and where the results are to be found. The functions 
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of the registers for the operations of arithmetic are indicated in 
Table 2.1. 

The Illiac has a fixed point arithmetic unit; the binary 
point is fixed so that any number x used in computation must lie 
in the range -1 < x < 1. The programmer must insure that all 
quantities remain within this range during a computation. The 
sign of a numerical quantity is indicated by the leftmost of the 
forty binary digits stored in a register of a memory location. 

The sign digit is 0 for a positive number 1 for a negative number. 
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TYPE 

INSTRUCTION 

INITIAL CONDITIONS 

FINAL CON 

DITIONS 

NOTES 

Specified 

Memory 

Location 

A 

. Register 

Q 

Register 

A 

Register 

Q 

Register 

L 

Add 

Addend 

Augend 

Arbitrary 

Sum 

Unchanged 

Augend can be 
set to 0 or l/2 

Subtract 

Subtrahend 

Minuend 

Arbitrary 

Difference 

Unchanged 

Minuend can be 
set to 0 or l/2 

7 

Clear 

Multiply 

Multipli¬ 

cand 

Set to 0 at 
start of 
i ns+.-rnnti on 

Multiplier 

Double Preci 

sion Product 

q^ set to 0 

Rounded 

Multiply 

tl 

Set to l/2 at 
start of in- 

Q +- r* +1 i rvn 

tl 

Rounded 

Product 

Usually 

Unwanted 

m 

Hold i 

Multiply 

n 

y 

Accumulant 

a 

It 

X 

xy + 2' 

‘ 39 a 

Tt 

6 

Divide 

Divisor 

Double Precision 
Dividend 

Residue 

Rounded 

Quotient 


0 

Left Shift 
(n 1) 

See Note 

2 

Double Precision 

Number x 

2x 

1. q Q unchanged 

2. Address digits 

1 

Right Shift 
(n 1) 

tf 

It 

X 

x/2 

specify num¬ 
ber (n) of shifts 


Table 2,1 Use of Illiac Registers During Arithmetic Instructions 




CHAPTER 3 


THE ORDER CODE 


The Illiac is a binary computer in which the storage 
capacity of each register or memory location is 4-0 binary digits. 

The orders which the machine carries out are represented by numbers 
in the machine. The relation between the order types and the corre¬ 
sponding numbers is called a code, and the collection of all such 
numbers is called the order code of the machine. The order code 
is interpreted by the control circuits of the machine and completely 
determines what the machine does. The machine is designed so that 
any storage location in the memory may be used either for orders or 
for numbers, the only distinction being that the control must be in¬ 
structed properly so that orders and numbers will be treated 
appropriately. 

3-1 THE MAKEUP OF ORDERS . An order for a digital computer 
consists of an instruction to say what to do and one or more addresses 
to say where to get the quantities to be used in carrying out the in¬ 
struction. In contrast to some existing computers, the Illiac uses 
what is called a one-address code. For a one-address code, the order 
is of the kind that says "add the number in memory location 12 to a 
number already in the arithmetic unit, leaving the sum in the arith¬ 
metic unit." 

Since a one-address code is used, it is not necessary to use 
40 binary digits to describe an order. The Illiac uses 20 digits and 
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packs two orders (an order pair ) into one location. These are the 


left-hand order and the right-hand order. Since many orders must 
refer to locations in the memory, each order contains an address . 

It is still called an address in those orders which do not refer 
to the memory. 

The electrostatic memory of the Illiac has 1024 locations 
and because 1024 = 2 10 we require 10 binary digits for the address. 
These are the rightmost 10 digits of the 20 digits assigned to an 
order. 

Of the remaining 10 digits of each order, eight are used 
for the instruction or (function ) and the other 2 are unused. The 
digit makeup of an order pair is shown in Figure 3.1. 


8-DIGIT 

INSTRUCTION 


10-DIGIT 
ADDRESS 


8-DIGIT 

INSTRUCTION 


10-DIGIT 

ADDRESS 


LEFT-HAND ORDER 


DIGITS 0-19 



RIGHT-HAND ORDER 


DIGITS 20-39 



Figure 3-1 
Order Pair Makeup 

The instruction digits are 8 in number because of the 
convenience obtained by using a base 16 (sexadecimal) number system 
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in which 4 binary digits may be represented by one sexadecimal digit. 
Thus each instruction may be coded as 2 sexadecimal digits. The 
symbols K. S, N, J, F, L are used for 10, 11, 12, 13, 14, 15- 

As an example of an order pair let us consider the following 
40 binary digits. 

1111010100000001110101000000000000000110 


When divided into instruction and address digits, these digits look 


like this: 

1111 0101 00 
INSTRUCTION 


* 4- f 

}■ i. } 9? 

OOOOOl'llQl 
* %' 
L.H„ ADDRESS 


0100 0000 00 0000000110 
INSTRUCTION R.H. ADDRESS 


The left-hand instruction is made up of the two 4-digit numbers 
1111 0101 which are the sexadecimal digits L5. The left-hand address 
is interpreted as an integer which may go from 0 to 1023 if we use 
decimal notation or from 0 to 3LL if we use sexadecimal notation. 

In sexadecimal representation the left-hand address is 1J 
which corresponds to the decimal number 29. Thus the left-hand order 
is L501J where the 0 has been supplied so that all 20 binary digits 
(including the unused 2) are accounted for. (We could have set the 
unused digits to l's and used N rather than 0 if we had wished.) 

Similarly the right-hand order is 40006, and we have, in 
sexadecimal notation, t he order pair 

L501J 40006 
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This order pair says, "Transfer the contents of memory location 1J 
to the accumulator; store the accumulator contents at location 6 
of the memory." 

It is quite inconvenient to have to write addresses in 
sexadecimal form, and it is unnecessary if the machine contains a 
program which will take addresses written in decimal form and con¬ 
vert them to sexadecimal (i.e., binary) form for machine use. A 
program of this type is available for the Illiac. It is called 
the Decimal Order Input Routine (See Chapter 5). 

3-2 EXECUTION OF ORDERS . The Illiac operates by using 
orders which have been stored in the memory according to a plan deter¬ 
mined by the programmer. The program is begun with a particular 
order chosen by the programmer. Let us suppose that it is the left- 
hand order at location 10. (We shall refer to addresses in decimal 
notation.) The control will put into R^ the order pair from location 
10. Then, until something is said to the contrary the control will 
follow a fixed pattern in executing orders. It will do the left-hand 
order and then the right-hand order in R^. Then it will put the order 
pair at location 11 into R^ and again do the left and right-hand orders. 
It will continue to withdraw and execute order pairs from successive 
memory locations until one of two things occurs: 

(a) One of the orders brought out says "stop," 

(b) One of the orders brought out says "change the 
sequencing." 
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Control Transfer Orders . The second Kind of order, a con ¬ 
trol transfer order , permits the programmer to change the sequencing 
of orders and provides the flexibility required for iterative processes. 
It works in the following way. Let us suppose that after the machine 
has executed the left-hand order at location 17 the programmer wishes 
to move to a sequence of orders beginning, say, with the right-hand 
order at location 35. Then the right-hand order at location 17 will 
say "Transfer control to the right-hand order at location 25-" The 
execution of this order will consist of arranging that the next order 
pair is brought to R^ from location 35 and that the left-hand order 
is skipped. Having done the right-hand order in R^, the control brings 
out the order pair from location 36 and proceeds in the usual way. 

Conditional Transfer Orders . There are two kinds of con¬ 
trol transfer orders, conditional transfer order s and unconditional 
transfer orders . We have just described the unconditional transfer 
order. The conditional transfer order does the same thing as the 
unconditional transfer order provided that the number in the accumu¬ 
lator is not negative . Otherwise it does nothing. In the example 
we have just used, if the right-hand order at location 17 bas been 
conditional, then if the accumulator had held zero or a positive 
number the next order executed would have been the right-hand order 
at location 35. But if the accumulator had held a negative number 
the next order executed would have been the left-hand order at location 
18 . 
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3„3 STOP ORDERS . There are four contr ol transfer orders, 
the right-hand and left-hand unconditional and c onditional orders. 

These have been further combined with a stop or n ot-stop order, giving 
eight orders in all. If the programmer chooses one of the four "stop 
and transfer control" orders the machine will stop before transferring 
control and must he started again with a panel switch. The switch can 
he set so that these stop orders are ignored, which of course makes 
these stop orders into ordinary control transfer orders. This is often 
useful in programming, and when it is done the machine can still he 
stopped with another kind of stop order which cannot be ignored. 

3.4 ORDER TYPES . The number of orders which the Illiac 
can execute is quite large (more than 100), hut not all of them are 
of general interest. What we shall do here is discuss the various 
types of orders, giving the variants of each. Following this dis¬ 
cussion will he a list of those Illiac orders which are of most use. 
This list is adequate for the coding of any problem, and beginners 
are urged to confine themselves to it. Experienced coders will find 
uses for other variants. 

The two sexadecimal instruction digits of an order give 
the order type and the order variant . We shall refer to them as the 
T- and V- digits, respectively. In the example L5 cited earlier the 
T- digit is L and denotes addition. The V- digit is 5 and chooses 
one of the variants of the addition order. 
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The order types are given in Table 3.1. 




T-Digit 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

K 

S 

N 

J 

F 

L 


Order Type 
Left Shift 
Right Shift 

Unconditional Transfer 

Conditional Transfer 

Store from A register 

Memory to Q, register 

Divide 

Multiply 

Input or Output 

Special Input or Output 

Increment Add from Q 

Add from Q register 

Not Used 

Extract 

Increment Add 

Add 


Table 3-1 
Order Types 


The meaning of the address digits of the different order 
types is given in Table 3*2 
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ORDER-TYPE 


ADDRESS SIGNIFICANCE 


0 , 1 
2 , 3 


4 

5 

6 
7 

8o, 82 

92 


K, S 


J 

F 

L 


Number of shifts (cannot exceed 63 ) 

Memory location from which next order 
pair will come 

Memory location at which storage will 
occur 

Memory location from which word is 
brought to Q register 

Memory location of divisor 

Memory location of multiplicand 

Number of binary digits to be input 
or output 

Character to be punched and number of 
punchings 

Address not used 

Same as for 5 order 

Memory location of addend 

Memory location of addend 

Table 3.2 

Meaning of Address Digits 


It will be seen from Table 3«2 that except for the 0, 1, 
8 ^9, K and S orders the address of an ord er always refers to the 
memory. 

3*5 ORDER VARIANTS . Let us now consider the variants 
obtained by changing the V-digit. The sexadecimal V-digit is made 
up of the four binary digits V8, V4, V2 and VI. There are 16 possi¬ 
ble combinations, giving V-digits from 0 to L, but not all are used. 
The results obtained with different V-digits are as follows: 
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( a ) The Digit VI . If VI = 1, A will always be cleared to 
zero at the beginning of any order. If VI * 0, A will not be cleared 
ffhus an odd V-digit means that A will be cleared. 

(b) The Digit V8. If VI « 1 and V8 « 1, the quantity l/2 
will be put in A at the beginning of any order. This is how rounded 
multiplication is carried out. Jf V8 ■ 1 and VI = 0, the Illiac will 
hang ^p. 

(c) The Digit V2 . This digit affects all o rders ex cept 
those of types 0, 1, 5 and J. See Table 3 .3. 

(d) The Digit V*K .This digit affects all o rders except 

those of types 0, 1, 5 and J. See Table 3.3. 

«, — _ ...... — . . 

The following notation is convenient for a more detailed 

description of the orders: 

A = accumulator register 
a Q = sign digit of A 
Q = multiplier-quotient register 
q Q = sign digit of Q, 

AQ = the 79 binary digit doiible register formed 
from A and Q by omitting q Q . 

N(R) = contents of register R 

N(n) = contents of memory location n 

If no ambiguity is possible, the symbols A and Q will also 
be used to denote the contents of A and of Q. 
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On the following pages there are detailed descriptions of 
the results obtained by changing the V-digit of the 15 order types 
used in the Illiac. After the detailed descriptions is an abbrev¬ 
iated list of orders. The orders in the abbreviated list have been 
underlined in the detailed descriptions. 
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3-11 


TYPE 

v4 

EFFECT OF V4 

V2 

EFFECT OF V2 

0, 1 

— 

NONE 

— 

NONE 


0 

RIGHT-HAND ORDER 

0 

STOP IF SWITCH SET TO OBEY 

2, 3 

1 

LEFT-HAND ORDER 

1 

NONE 



NONE 

0 

STORE FULL WORD 

4 

0 

RIGHT-HAND ADDRESS 

1 

STORE ADDRESS 


1 

LEFT-HAND ADDRESS 



5, J 


NONE 

_ 

NONE 

6 

1 

V4 MUST BE 1 

1 

V2 MUST BE 1 

8 

0 

TAPE 

0 

INPUT 

1 

DRUM 

1 

OUTPUT 




0 

INPUT 

9 

0 

TAPE 

1 

OUTPUT 

7, K, S 

0 

SUBTRACT 

0 

NUMBER 

F, L 

1 

ADD 

1 

ABSOLUTE VALUE 


Table 3-3 Effect of Digits V4 and V2 
























OV n 


Left Shift (Double) 
Final Stop 


3 Orders 
1 Order 


If n = 0, the machine will stop. If not, repeat n times 
the operation which replaces the contents 

a 0 , B S , , a 3 g, a^j g y ... , q^, , 

of AQ by 


a x , a 2 , a 3 , , a^ q 2 , q 4 , . . . , q^, 0 


leaving q Q unchanged. 

The number n will be interpreted modulo 64. 


Variants 


0 , 2 , 4 , 6 

1, 3, 5, 7 

8, K, N 

9, 5, J, L 

F 


Order as described above. 

Clear A and then execute as described. 

Illiac will hang up. Avoid these. 

Clear A, insert l/2 in A, and then execute 
order as described above. 

Final Stop. 


Use only _00, 01, _0g for shifts 
Use OF with n = 0 for final stop 


39 
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IV n 


Right Shift (Halve) 


3 Orders 


If n = 0, the machine will stop. If not, repeat n times 
the operation which replaces the contents 

a 0' a i' a 2> * * ’ ' a 38> a 39 ; q l' q 2> q 3' ' * * ' q 38' q 39 

of AQ by 

a 0 / a 0' V • • • ' a 37' a 38’ a 39’ q l' q 2' • * * ' q 37' q 38' 

leaving q Q unchanged. 

The number n will be interpreted modulo 64, 

Varients 

Order as described above. 

Clear A and then execute as described. 

Illiac will hang up. Avoid these. 

Clear A, insert l/2 in A, and then 
execute order as described above. 

Use only 10, 11, 19 


0, 2, 4, 6 
1, 3, 5, 7 

8, K, N, F 

9, S, J, L 
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2V n 


Unconditional Control Transfer 


12 Orders 


Bring the next order pair from memory location n and 
choose the left or right hand order of this pair, stopping before¬ 
hand or not, depending upon V. 


Variants 


0 


2 

4 

6 

l, 3, 5, 7 

8, K, N. F 

9, S, J, L 


Stop. The first order after starting with 
the START switch will be the right-hand order 
at memory location n. The stop can be ignored 
by setting a panel switch. 

Transfer control to right-hand order at memory 
location. 

Same as 0 except take left-hand order. 

Same as 2 except take left-hand order. 

Same as 0, 2, 4, 6 except clear A first. 

Illiac will hang up. Avoid these. 

Same as 1, 3, 5, 7 except also make A = l/2 
after clearing. 


Use 20, 21, 22, 23, 24, 25, 26, 


27, 29, 2S, 2J, 2L. 


Starting After A Stop . When the Illiac has been stopped 
by one of the control transfer stop orders, it is usually started 
again by moving the black switch to START, from which position the 
switch automatically returns to OBEY. 

The Illiac can also be started again by moving the white 
switch through EXECUTE to FETCH and then back to RUN. If this is 
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done, the control transfer order which stopped the Illiac will be 
ignored. The normal sequencing will then follow unless the stop 
order is a right hand order transferring control to the right hand 
side of a word. In this case, the order first obeyed after starting 
will be the right hand instead of the left hand order of the new 
order pair brought out. 

For example consider the following orders;. 


p L5 F 


20 p+2 

p+1 40 IF 

24 p 

p+2 7J 2F 

L4 3F 


If we stop with the 20 order, the black switch will start with 
L4 3F and the white switch with 24 p. If we stop with the 24 
order, the black switch will start with L5 F and the white 
switch with 7J 2F. 
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3V n 


Conditional Control Transfer 


4 Orders 


If A > 0, bring the next order pair from memory location 
n and choose the left or right hand order, of this pair, stopping be¬ 
forehand or not, depending upon V. If A < 0, go on to the next 
order, 

Variants 

0, 2, 4, 6 If A > 0, do the same operation as for 

the corresponding 2V order. If A < 0, 
go on to the next order. 

1, 3, 5 s 7 Identical with corresponding 2V orders. 

8, K, N, F Illiac will hang up. Avoid these, 

9, 8, J, L Identical with corresponding 2V orders . 

Use only 30, 32, 34, 36. 

Starting After A Stop. The discussion given with the 
2V orders applies here to the corresponding 3V orders. 
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4v n 


Store 


9 Orders 


Copy into memory location n all of the contents of A, 
the contents corresponding to the address of a left-hand order, or 
the contents corresponding to the address of a right-hand order, 
depending upon V. 


0 , k 
l, 5 
2 

3 

6 

7 

8, K, N, F 

9, J 
s 


Variants 

Replace N(n) by A. 

Replace N(n) and A by 0. 

Replace address digits of right-hand order 
at memory location n by the corresponding 
digits of A. 

Same as 2 except clear A first. 

Same as 2 except take left-hand order. 

Same as 6 except clear A first. 

Illiac will hang up. Avoid these. 

Replace N(n) and A by l/2. 

Replace A by l/2 and address digits of 
right-hand order at memory location n 
by 0. 

Same as S except take left-hand order. 


L 



5V n 


Memory to Q 


3 Orders 


Transfer N(n) to Q, 

Variants 

0, 2, 4, 6 Transfer N(n) to Q 

3, 5, 7 Clear A and transfer N(n) to Q. 

8, K, N, F Illiac will hang up. Avoid these. 

9, S, J, L Put l/2 in A and transfer N(n) to Q. 

Use only 50♦ 51/ 59* 
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6V n 


Divide 


3 Orders 


Divide AQ by N(n), placing the rounded quotient in Q, 

(the least significant digit being 1 for the roundoff) and leaving 
a residue in A. If |a| > |N(n)| the Illiac will stop after 
dividing. If |A| = |N(n)| and if A > 0, the Illiac will stop 
after dividing; if |A| = |N(n)| and if A < 0, the Illiac will 
not stop after dividing. 


6 

7 

L, S 

8, K, N, F 
0, 1, 2 , 3 
b, 5, 9, J 


Variants 

As described above. 

Make A = 0, then proceed as above. 

Make A = l/2, then proceed as above. 
Illiac will hang up. Avoid these. 

These give incorrect results or results 
which are correct only under certain 
conditions. Avoid them. 


Use only 66, 67, 6l. 
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7V n 


Multiply 


12 Orders 


Put Q x P(n) + 2 ^ A into AQ, the least significant 
39 digits being in Q with = 0. 

Variants 


0 

1 

2 

3 

5 

6 

7 

8, K, N, F 

9 

s 

J 

L 


P(n) = - N(n) 

P(n) = - N(n); A = 0 "Unrounded Negative Multiply" 
P(n) = - |N(n)| 

P(n) = - |N(n)| j A = 0 
P(n) = N(n) "Hold Multiply" 

P(n) = N(n) ; A = 0 "Unrounded Multiply" 

P(n) = |N(n)| 

P(n) = |N(n)| ; A = 0 


Tlliac wi 1J—han g up . Avoid these. 

P(n) = -N(n) j A = l/2 "Rounded Negative Multiply 


It 


P(n) = - |N(n)| j A = 1/2 

P(n) = N(n) j A = l/2 "Rounded Multiply" 

P(n) = |N(n) | ; A = 1/2 


Use only 70, 71, 72, 73, 74, 75, 
76, 7J, 79, 7S, 7J, 7L 
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8V n 


Input-Output 


9 Orders 


Transfer words between A and the input tape, output punch, 
or magnetic drum. 

The address n must be a multiple of 4 for the tape and punch 
orders and must be 11 for drum orders. 

Variants 

0 Shift AQ four places left as in the 00 order 

and replace a^£> a 39 ^ ^ nai T 

digits corresponding to the sexadecimal 
character being read. This is done n/4 times. 

1 , 9 Clear A and then do as in 80 order. 

2 Punch the digits a Q , a^, a 2 , a^ as one sexa¬ 
decimal character and shift AQ four places 
left as in the 00 order. This is done n/4 
times. 

3 Clear A and do as in 82 order. 

5 This is a 4o-digit order of the form 85 11 TV p. 

We distinguish two cases. 

( 1 ) T is not 0, 1, 8 , 9" In this case, A 
and Q are shifted left eleven places as 
in the 00 order and the word at drum 
location p is placed in A. Then the T 
order is executed using address p. Com¬ 
plete freedom is not available in drum 
addresses because p may interfere with V, 
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6 


7 

s 

L 

8, K, N, F 
h, J 


( 2 ) T is 0, 1, 8 , 9 . In this case after 
the word at drum location p is placed 
in A the right-hand order is skipped. 

This permits use of any drum address 
for p. 

This is a 40-digit order of the form 85 11 TV p. 

We again have two cases: 

(1) T is not 0, 1 , 8 , 9- In this case A is 
transferred to drum location p, and A 
and Q are shifted 11 places left as in 
the 00 order. Then the T order is ex¬ 
ecuted using address p. Complete freedom 
is not available in drum addresses be¬ 
cause p may interfere with V. 

( 2 ) T is 0, 1, 8 , 9 . In this case the right- 
hand order is skipped after doing the 
left-hand order as in case (l). This 
permits use of any drum address p. 

Same as 86 except clear A first. 

Put l/2 in A and do as in 82 order. 

Put l/2 in A ana do as in 86 order. 

Illiac will h ang mx.— A - v o M- t he se, 

These are not useful. Avoid them. 

Use 80, 81 , 82, 83 , 85 , 86 , 87 , 8 S, 8 L. 
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9V n Special Input-Output 2 Orders 

Variants 

1 Five hole input. Clear A, shift AQ, four places right 
and replace aa^, a^g, a^ by the binary digits 
corresponding to the four least significant holes on 
the tape. Place the contents of the fifth hole in 
position a Q . The address part of this order must be 4. 

2 Letter output. Punch on the tape a character de¬ 
pending upon the address digits n. Three quantities 
are defined by the 10 binary address digits: 

(l) The leftmost 4 digits define the usual 4 digit 
positions in the output tape. 

(2 ) The rightmost digit defines the 5th hole in 
the output tape. 

(5) The rightmost 6 digits determine the number 
b of times that the above-defined character 
is punched and also the number of right 
shifts executed. The number of characters 
punched will be found by dividing the number 
in the rightmost 6 digits by 4 and rounding 
up to the next integer. 

The address n may always be found from the following 
formula: 

n = 64a + 4b + c -2 

where a is the character punched, a = 0, 1, 2, . . ., J, F, L 

b is the number of characters punched, 1 < b < l6 
c determines the fifth hole, c =0, 1. 
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The number of right shifts executed is 4b + c - 2. 

Example . Punch the character 7 thirteen times. 

n = 64x7 + 4xlj) + 0 - 2 = 498 

It will be found that the last 6 digits contain the 
number yj which when divided by 4 and rounded up gives 13 . There 
will be 50 right shifts. 

Figure 3.2 shows the relationship between the tape 
holes and the address digits of the 92 order. The address shown 
will print the character 7 thirteen times. 



Figure 3-2 

Address Digits of 92 Order 


3-24 




KV n 


Increment Add from Q 


12 Orders 


These orders are identical with the FV orders except 
that Q, is used instead of N(n). 


Variants 


0 

Add -Q - 2" 59 to A. 

(See note). 




1 

Put -Q - 2 -59 in A. 








2 

Same as 0 if Q > Oj 

same 

as 

4 

if 

Q 

< 

0. 

3 

Same as 1 if Q > 0; 

same 

as 

5 

if 

Q 

< 

0. 

4 

Add Q, + 2“ 59 to A. 








5 

-39 

Put Q + 2 ^ to A. 








6 

Same as 4 if Q > Oj 

same 

as 

0 

if 

Q 

< 

0. 

7 

Same as 5 if Q > 0» 

same 

as 

1 

if 

Q 

< 

0. 

8, K, N, F 

T11 i an wi 1 1 hang-up. 

-—-Aaim.ri 

.these. 



9 

Put -Q - 2 -59 + 1/2 

in A 







S 

Same as 9 if Q > 0; 

same 

as 

J 

if 

Q 

< 

0 . 

J 

Put Q + 2~ 59 + 1/2 : 

in A. 







L 

Same as J if Q > 0; 

same 

as 

9 

if 

Q 

< 

0. 


Use KO, Kl, K2, K3, K4, K5, K 6 , K7, K9, KS, KJ, KL 

NOTE: -Q - 2 ^ is the digitwise complement of Q. 
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SV n 


Add from Q 


12 Orders 


Q is added, subtracted, etc., to A. These orders 
identical with the LV orders with Q used instead of N(n). 


Variants 


0 

1 

2 

3 

4 

5 

6 

7 

8, K, N, F 

9 

s 

J 

L 


Subtract Q from A. 

Put -q in A. 

Subtract |q| from A 
Put " | Q | in A. 

Add Q to A. 

Put Q in A. 

Add |q) to A. 

Put |q| in A. 

^lliac will hang up. Avoid the se 
Put l/2 - Q, in A. 

Put 1/2 - |q| in A. 

Put 1/2 + Q in A. 

Put l/2 + |q| in A. 


Use SO, SI, S2, S3, S4, S5, S6, S7, S9, SS, SJ , SL 


are 
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JV n 


Extract 


3 Orders 


If two corresponding digits of N(n) and Q are both 
l's, put I in that place of Q. Otherwise put 0. This order gives 
the logical product of N(n) and Q. 

Variants 

As described above. 

Cl ear A and do JO order. 

Illiac will hang up. Avoid these . 

Put l/2 in A and do JO order. 

Use only JO, Jl, J9 


0 , 2 , 4, 6 
1, 3, 3, 7 

8, K, N, F 

9 , S, J, L 
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12 Orders 
1 Order 


FV n 


0 

1 

2 

3 

4 

5 

6 

7 

8, K, N 

9 

s 

J 

L 

F 


Increment Add 
Programmed Stop 


Variants 


Add -N(n) - 2~ 39 to A. 
Put -N(n) - 2~ 39 in A. 
Same as 0 if N(n) > 0; 
N(n) < 0. 

Same as 1 if N(n) > 0; 
N(n) < 0. 

Add N(n) + 2~ 39 in A. 
Put N(n) + 2~ 39 in A. 
Same as 4 if K(n) > 0; 
N(n) < 0. 

Same as 5 if N(n) > 0; 
N(n) < 0. 


(See note) 
same as 4 if 
same as 5 if 


same as 0 if 

same as 1 if 


Il liac will hang up. Avoid these. 
Put -N(n) - 2" 39 + l/2 in A. 

Same as 9 if N(n) > 0; same as J if 
N(n) < 0. 

Put N(n) + 2*" 39 + l/2 in A. 

Same as J if N(n) > 0; same as 9 if 
N(n) < 0. 


Stop. This stop is used to indicate 
failure of a programmed check. 


NOTE: 


Use F0, FI, F2, F3, F4, F 6, F'J, F9, FS, FJ, FF, FL. 

-N(n) - 2 39 is the digitwise complement of N(n). 
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LV n 


Add 


12 Orders 


N(n) is added, subtracted, etc., to A. 


0 

1 

2 

3 

k 

5 

6 

7 

8, K, N, F 

9 

s 

J 

L 


Variants 


Subtract N(n) from A. 
Put -N(n) in A. 
Subtract |N(n)| from/ 
Put - |N(n)| in A. 

Add N(n) to A. 

Put N(n) in A. 

Add |N(n)| to A. 

Put |N(n)| in A. 
Illiac will hang up. 
Put l/2 - N(n) in A. 
Put 1/2 - |N(n)| in A 
Put l/2 + N(n) in A. 
Put l/2 a JN(n)| in A 


Use LO, LI, L2, L3, iA, L5, L6, L7, L9 


Avoid these. 


LS, LJ, LL. 
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ABBREVIATED ORDER LIST 


ORDER 


DESCRIPTION 


00 n 
09 n 

OF 0 
10 n 

19 n 

20 n 


22 n 

2k n 
26 n 


Shift AQ left n places, 1 < n < 63 . 

Make A = l/2 and then shift AQ left n 
places, 1 < n < 63 . 

Final stop. 

Shift AQ right n places, 1 < n < 63 , 

Make A = l/2 and shift AQ right n places 
so that AQ contains 2“ n_1 , 1 < n < 63 . 

Stop. The first order after START will 
be the right hand order at location n. 
The stop can be ignored by setting the 
panel switch to IGNORE. 

Transfer control to the right hand order 
at location n. 

Same as 20 except take left hand order. 
Same as 22 except take left hand order. 


30 n 
32 n 

3k n 
36 n 

kO n 
41 n 


If A > 0 do as in the corresponding 2V 
order. 

If A < 0 go on to the next order. 


Replace N(n) by A. Do not change A. 
Replace N(n) and A by zero. 
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ORDER 


DESCRIPTION 


42 n 


46 n 


k -9 n 
50 n 
66 n 

70 n 

71 n 

72 n 
75 n 

74 n 

75 n 

76 n 

77 n 

79 n 
7 S n 
7J n 
7 L n 

80 n 

8 1 n 


Replace address digits of the right hand 
order at location n by the corresponding 
digits of A. Do not change A. 

Replace address digits of the left hand 
order at location n by the corresponding 
digits of A. Do not change A. 

Replace N(n) and A by l/ 2 . 

Replace Q by N(n). 

Divide A + 2 99 (Q + Iq) by N(n). 

Put -N(n) Q + 2 - " A into AQ. 

Multiply -N(n) by Q, putting result in AQ. 

.Put - |N(n) | Q + 2 _59 A into AQ. 

Multiply - |N(n)| by Q, putting result in AQ. 

Put N(n) Q + 2 ~ 59 A into AQ. 

Multiply N(n) by Q, putting result in AQ. 

Put |N(n)| Q + 2 -59 A into AQ. 

Multiply |N(n)| by Q, putting result in AQ. 

Put the rounded product -N(n) Q into A. 

Put the rounded product - |N(n)| Q into A. 

Put the rounded product N(n) Q into A. 

Put the rounded product |N(n) | Q into A. 

Input n /4 sexadecimal characters from the tape. 

Clear A and input n /4 sexadecimal characters 
from the tape. 
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ORDER 


DESCRIPTION 

82 n 

Punch n/4 sexadecimal characters on the tape 

85 n 

Repla.ce A by the word at drum location n. 

See page 3-21. 

86 n 

Replace the word at drum location n by A. 

See page 3-22. 

91 4 

Read one five-hole character from the tape. 

92 131 

Punch one carriage return and line feed 
character. 

92 64 2 

Punch one + sign (or K). 

92 706 

Punch one - sign (or S). 

92 769 

Advance film. 

92 963 

Punch one space character. 

KL n 

Put - Q - 2 ^ in A. (This is the digitwise 
complement of Q). 

K4 n 

Add 

Q + 2' 39 to A. 

K5 n 

Put 

Q + 2“ 59 in A. 

SO n 

Subtract Q from A. 

SI n 

Put 

- Q in A. 

S2 n 

Subtract |q| from A. 

S3 n 

Put ■ 

- |Q 1 in A. 

S4 n 

Add Q to A. 

S5 n 

Put Q in A. 

S6 n 

Add | 

[QJ to A. 

S7 n 

Put 

|Q| in A. 
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ORDER 

DESCRIPTION 


JO n 

Replace Q by the logical product 

N(n). 

of Q and 

FI n 

Put -N(n) - 2 -59 in A. (This is 
wise complement of N(n). 

the digit 

F^ n 

Add N(n) + 2 -59 to A. 


F5 n 

Put N(n) + 2 -39 in A. 


LO n 

Subtract N(n) from A. 


LI n 

Put -N(n) in A. 


L2 n 

Subtract |N(n)| from A. 


L3 n 

Put - |N(n) | in A. 


Lk n 

Add N(n) to A. 


L5 n 

Put N(n) in A. 


L 6 n 

Add |N(n)| to A. 


L7 n 

Put |N(n ) | in A. 
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CHAPTER k 


ROUTINES 


One of the first tasks in programming a computation is 
to break it down into a number of small or medium size operations, 
each of which forms a fairly distinct logical step in the whole 
calculation. Such a step may be the evaluation of a function such 
as a square root, logarithm or cosine. Each step must be carefully 
defined so that it fits correctly with all the other steps. Having 
decided exactly what is required of each step, it is possible to 
proceed with the coding of the individual steps. 

This procedure has three important advantages. Firstly, 
it enables the programmer to concentrate on one part of the job at 
a time. Secondly, the coding of each step can be tested separately 
before incorporating it in the program. Thirdly, certain steps 
are common to a very large number of different calculations; these 
have been coded and tested once for all, and may be used by anyone 
in the Laboratory. 

The coding for one step constitutes a routine . The 
collection of routines for performing common operations in the Illiac 
is called the library of routines. It includes routines for evaluating 
many simple functions like those mentioned above, also for printing 
numbers in various layouts, for integrating differential equations, 
for solving simultaneous equations, and for many other operations. 
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A routine, when used to perform part of the work of another 
routine, is called a subroutine . The simplest type of subroutine 
is merely a string of order pairs which can be inserted in the appro¬ 
priate place among the other orders of a program. This is called an 
open subroutine. However, for various reasons another type, called 
a dosed, subroutine, is more commonly used. 

4.1 CLOSED SUBROUTINES . A closed subroutine is also 
a string of consecutive words but these do not have to be placed 
among the other orders of the program. Instead, they may be put 
in any convenient place in the store. Each time the subroutine 
is to be executed, control must be transferred to it in a certain 
special way (the subroutine is said to be entered ). It is so arranged 
that when execution is complete, control is automatically returned zo 
the point from which the subroutine was entered, so that execution of 
the rest of the program may continue. 

In this way a program is seen to consist of several dis¬ 
tinct, self-contained blocks, namely the various subroutines and 
the part of the program (usually called the main program or master 
routine) which makes use of its subroutines by sending control to 
them. Internal rearrangement of a routine is usually difficult, 
but the routines comprising a program can be shuffled about very 
easily, and this fact makes the coding of large problems much simpler. 

It is not necessary to enter each subroutine directly from 
the master routine; there is nothing to prevent one subroutine being 
entered from another. A subroutine may itself have subroutines. 
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4.2 ENTERING A CLOSED SUBROUTINE. The following two 


orders must he inserted in the master routine at the point from 
which a subroutine is required to be entered. 

k Apy 
50 k'\ 

k+1 26 y Subroutine starts at m. 

Any 

These orders must be, as shown, in the right-hand half of one word 
and the left-hand half of the next (it may be necessary to waste an 
order to do this). The 50 order contains the address of the memory 
location in which it is itself contained. The 2 6 order contains the 
address of the memory location containing the first word of the sub¬ 
routine (it actually transfers control to the left-hand half of this 
word). 

The subroutine, after execution, automatically transfers 
control to the order following these, i.e., to the right-hand order 
in memory location (k+l), so that there is no break in the logical 
continuity of the master routine. There is no need for the user to 
know how the subroutine does this, but it is not difficult to under¬ 
stand and it helps to complete the picture. 

4.3 RETURNING CONTROL TO THE MASTER ROUTINE . It will 
be seen that the effect of the 50 order above is to place the order 
pair, of which it is the right-hand half, in the Q register. Thus 
whenever a closed subroutine is entered, the Q register contains a 
pair of orders of which the right-hand order contains the address 
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of the memory location from which the pair came. Control must be 
returned to the right-hand order in the memory location following 
this. 

The subroutine uses the information in the Q register to 
set the address in a transfer of control order (a 22 order) called 
the link , which ultimately causes the transfer back to the master 
routine. It does this by means of the orders shown in Table 4.1. 


m K5 0 

42 m+n 

m+1 


These orders "plant" 
the link 

a 

\ 


m+n 

22 (0» link 

Table 4.1 

Forming a Closed Subroutine Link 


The K5 order transfers the order pair from the Q register 
to the A register, increasing the right-hand address by 1 as it 
does so. This address is now the address to which control is 
eventually to be returned. The 42 order places this address in 
the link. The link is situated so that it will be encountered 
by control when the execution of the subroutine is complete. Note 
that the link should be a right-hand order so that its address may 
be Inserted from the right-hand side of the A register. 

It may sometimes be convenient to use a 32 order as a link 




this order may then perform a test within the subroutine and return 
control to the master routine only if some condition is satisfied. 
Also the above procedure for "planting" the link may be varied if 
desired, so long as the result is the same. In particular note the 
variants in sections 4.5 and 4.6 below. 


4.4 PLACING THE ARGUMENT . All subroutines operate on 
at least one number somewhere in the machine, and there must be 
agreement between the subroutine and the master routine on the 
placing of these numbers, and also on the placing of the results 

of the subroutine operation. 

If only one number is operated upon, it is convenient 
to use the A register to hold this number when the subroutine is 


entered(since the A register is not being used for anything else 
at this moment). Similarly, if only one number is produced by the 
subroutine it is convenient for the subroutine to leave this number 
in the A register. 

For example. Library Routine R 1 is a closed subroutine that 
finds the square root of the number given in A, and leaves this 
square root in A. Suppose we wish the Illiac to find the largest 

root of 


4 5 2 

x +2x +x - c = 0 

where c = N(lO) and satisifies 0 < c < l/2• It can be shown that 
the root is 


x = -l/2 + y/l/4 + 7c • 
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Suppose that the constants l/2 and l/4 are given in 20 and 21 
respectively, and that x is required in 11. Routine R 1 will be 
used to find the square roots; let its first word be in 100. Then 
the master routine might contain the orders of Table 4.2 


50 

L5 

10 

Put c in A 


50 

50 

Enter Code R 1 

51 

26 

100 

to form sJ~o~ 


L4 

21 

Add l/4 to ftr 

52 

22 

52 

Waste order 


50 

52 

Enter Code R 1 

55 

26 

100 

to form J 1/4 + 


L0 

20 

Subtract l/2 to form x 

54 

4o 

11 

Put x in 11 



Table 

4.2 


Master Routine Using A Closed Subroutine 

The subroutine itself must use a slightly different method 
oi link planting from that given in section 4.5, which would des¬ 
troy the argument in A before using it. It is here necessary to 
rescue the argument and hold it in a storage location until the 
link has been planted, as in Table 4.5 


m 

4o 

<1 

Put argument in q 


K5 

0 


m+1 

42 

m+n 

Plant link 

m+n 

# , 

• • • 

• # # 



22 

(o) 

Link 


Table 4.5 

Modified Link Planting 
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4-5 PROGRAM PARAMETERS. Sometimes a subroutine is 


made to carry out slightly different operations on different occasions. 
For example, one subroutine might be made to print numbers to any 
number of figures. Another might form the n ^ root of a number, 
where n is any positive integer. The value of n, or the number of 
figures to be printed, is called a parameter of the subroutine. 

The value of a parameter must always be specified if the 
subroutine is to operate correctly. There are two standard ways of 
specifying parameters} parameters specified in the way now to be 
described are called program parameters . (The others are called 
"preset" parameters - see below, section 4.7). 

A program parameter is specified in the first half of 
the word containing the 50 order which is required on entering any 
closed subroutine. Let us take for example Library Routine R 2, 
which is a closed subroutine for replacing N(A) by its p root. 

Here p is a program parameter, to be specified as follows: 

k 50 P Program parameter 

50 k 

k+1 26 m Subroutine starts at m 

Any 

To illustrate this, suppose we wish the machine to com¬ 
pute the positive real root of 

» ’ 

x^ - 3 x:^ + + 2 c = 0 , (l/ 2 <c<l) 

which is given by the formula 
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Suppose that c = N(lO) and that x is required in 11} suppose also 
that N(20) = -1, and that Routine R 2 starts at 100. Then the master 
routine might contain the program given in Table 4-.4 


50 

L9 

10 

N(A) = (1/2) - c 


00 

1 

N(A) = 1 - 2c 

51 

50 

3 

Program parameter 


50 

51 

Link 

52 

26 

100 

To Routine R 2 


L0 

20 

N(A) = 1 + 1 - 2c 

53 

50 

2 

Program parameter 


50 

53 

Link 

54 

26 

100 

To Routine R 2 


40 

11 

x to 11 


Table 4.4 

Use of Program Parameter 

It will be noticed in Table 4.4 that the order contain¬ 
ing the parameter (50 3 or 50 2 ) is actually obeyed by the machine 
before the subroutine is entered. This is just a waste of time 
which is not worth avoiding} it does no good and, provided the order 
containing the parameter is suitably chosen, no harm. It is usual 
to use a 50 order here; this puts rubbish in the Q register, but 
the latter is then immediately reset by the following 50 order and 
no damage has been done. 

A more elaborate example of the use of a program parameter 
is to be found in Library Routine P 1 which prints A as an integer 
or fraction. The entry is by means of the orders of Table 4.5- 
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k XY d 

50 k 
k+1 26 m 

Table 4.5 

Multiple Program Parameters 


In Table 4.5 x = 5 or J; if x = 5 and A is negative the printed 
number is preceded by a minus sign, otherwise it is preceded by 
a space. Also Y = 0 or 2; if Y = 0, A x 2 39 is printed as an 
integer; if Y = 2, A is printed as a fraction (correctly rounded 
off). The layout of digits is given by d = lOq + s; q digits are 
printed with a space after the first s. 

It will be noted that all combinations of X and Y form 

harmless orders. ^ 1 

4.6 INTERPRETIVE ROUTINES .' There is a type of sub¬ 


routine which, instead of executing a single distinctive operation, 
carries out a whole series of operations. Each operation requires 
a parameter for its specification, so that the master routine contains 
a string of parameters, one for each operation. The string may be 
of indefinite length. 


Such subroutines are called interpretive routines. Their 


use lies in programs that involve actual operations on elements which 
are not numbers stored in the usual form but may be numbers stored 
in some special form, or different mathematical entities altogether 
such as expansions in Boolean algebra. The commonest application 
is to numbers stored in the so-called "floating point' form (see 
Chapter 6). There are certain advantages in storing numbers in this 
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form; however, such numbers cannot be added or multiplied in a single 
Illiac operation. A routine is required to handle simple arithmetic, 
and for this purpose an interpretive routine is used (e.g. Routine Al). 
Each parameter of the interpretive routine corresponds to one arith¬ 
metical operation, just as each Illiac order corresponds to one opera¬ 
tion in the Illiac. Thus, for example, to place in register 6 the 
sum of the numbers in registers 4 and 5 we write the parameters 

85 4 
84 5 
8s 6, 

which act in a similar way to the Illiac orders 

L 5 4 
L4 5 
4o 6 . 

Owing to the close similarity between the parameters and 
the ordinary Illiac orders, the parameters are themselves often re¬ 
ferred to as interpretive orders, or merely as orders, "orders," or 
orders . The "order code" of the interpretive routine may be des¬ 
cribed in a similar way to the order code of the machine] it involves 
reference to an "accumulator" which behaves, for floating point num¬ 
bers, like the A register of the Illiac. 

An interpretive routine is entered in the same way as a 
closed subroutine, but the parameters (i.e., the interpretive "orders") 
follow the orders causing entry. Finally a special parameter may 
be used to cause control to the Illiac to be transferred out of the 
interpretive routine and back to the master routine. Thus if we 
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imagine that the above example is a complete set of operations to 
be carried out by the interpretive routine, the master routine 
would contain the program of Table 4.6. 


k 

Any 




50 

k 

Enter interpretive routine 

k+1 

26 

m 



85 

4 

Parameters (interpretive orders) 

k+2 

84 

5 



8s 

6 


k+3 

8j 

00 

k+4 

0 

Send control of Illiac to L.H. 
side of k+4 


etc 

9 



Table 4.6 

Program for Interpretive Routine 

4.7 PRESET PARAMETERS . By making use of a library routine 
a programmer avoids not only the necessity of writing out the orders, 
but also the labor of punching them, since a master copy of the tape 
is kept on file in the Teletype Room to be copied when required. 

If the routine has one or more parameters associated with 
it, it can meet a wide variety of requirements. However, the use 
of program parameters consumes both computing time and storage space. 

A program parameter has the property that it can be varied 
from one application to another within the same program. Frequently 
it happens that, although the ability to choose a value of a para¬ 
meter to suit any particular program is desired, the ability to change 
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the value during the execution of a program is not needed. The 
value can be set before execution begins; there is no need for the 
orders which, in the case of a program parameter, set the value 
afresh each time the routine is entered. A parameter whose value 
is set before execution begins is called a preset parameter. 

The master copy of a library routine must be valid for 
all values of any preset parameters involved. The fixing of a pre¬ 
set parameter for a particular program must therefore be done after 
(or during) the copying of the master tape; in fact it is done as 
the program is read into the Illiac. This operation is carried 
out by the Decimal Order Input Routine which is described in 
Chapter 5- 
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CHAPTER 5 


THE DECIMAL ORDER INPUT 

The Illiac is a binary machine intended for University 
research. Since it will be used by a large group of people with 
various scientific backgrounds, it is important that its use be 
made as simple as possible. To this end a library of subroutines 
has been organized and aids in using it have been devised. One of 
these aids is the Decimal Order Input, Library Code XI. This code 
has two principal purposes: 

(1) To make it possible for programmers 
to use decimal notation in coding, 

(2) To make use of the library easy. 

The use of subroutines is discussed more fully in Chapter 
4 , but one way to characterize a subroutine is to say that it is 
essentially an extension of the order code of the machine. It is 
a group of orders used to carry out one or more operations which 
may be as simple as a square root or as complicated as a complete 
program for executing the details of floating point arithmetic. 

In any case the interior structure of a subroutine will depend upon 
its location in the memory. For example, the left-hand address of 
the second word in the square root routine (Library Code Rl) refers 
to the eighth word of the routine. If the routine begins at loca¬ 
tion 10, this address must be 18 . But if the routine begins at 
location 97, this address must be 105 . The use of subroutines is 
very awkward unless problems of address changing can be easily 
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handled. 


Another principal nuisance in coding is the number system 
used by the Illiac. In the sexadecimal system using Illiac notation 
(see Chapter 2), the 1024 addresses of the memory are represented 
by numbers between 0 and 3LL* It is much more convenient if addresses 
can be expressed in decimal form with the Illiac doing the necessary 
converting to the binary system. 

With the Decimal Order Input every number corresponding 
to the address part of an order must be written in the decimal system. 
These decimal numbers are always converted to sexadecimal (i.e., binary) 
form by the Decimal Order Input. 

The Decimal Order Input overcomes the obstacles of changing 
addresses in subroutines and of converting from decimal to binary 
form. It also provides certain other useful services. 

5.1 RELATIVE AM? FIXED ADDRESSES . If we can devise some 
way of using addresses in decimal form, we never need the 6 characters 
K, S, N, J, P, L In an address and can use them for other purposes. 

Let us consider a program beginning in location 11 of the memory, as 
s-hown in Table 5«1* 


11 

4o 

1 


S5 

1 

12 

Lb 

13 


46 

20 

13 

51 

1 


10 

1 


Table 

5.1 


Program Beginning at Location 11 
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The same program, if begun at location 20, would read as in Table 

5.2. 


20 

4o 

1 


S5 

1 

21 

l4 

22 


46 

29 

22 

51 

1 


10 

1 


Table 5.2 

Program Beginning at Location 20 

The second word .has changed. But if we mark each address which de¬ 
pends upon the location of the routine with the symbol L and each 
address which is independent of the location of the routine with 
the symbol F, we have the program in Table 5*3- 


0 

40 

IF 


S5 

LF 

1 

L4 

2L 


46 

9L 

2 

51 

IF 


10 

IF 


Table 5.3 

Program With Relative and Fixed Addresses 

The addresses ending in L are relative to the location of the first 
word of the routine. 
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We can now write each subroutine as if it begins at location 
zero provided we add to each L-terminated address the location of the 
first word of the routine when we store the routine. We shall do this 
by introducing a new symbol K. 

5.2 DIRECTIVES. Orders which are written for the Illiac 
use two sexadecimal function digits followed by an address. We 
follow this convention even in the case of pseudo orders. A directive 
is such a pseudo order. When we want a group of orders to be placed 
in memory locations n, n+ 1 , n+ 2 , . . . , we punch on the tape the 
directive 

00 nK ^ . 

The Decimal Order Input will recognize this as a directive and will 
place the orders following the directive in pairs in the locations 
n, n+ 1 , n+ 2 , . „ . , starting with the left-hand side of location n. 

It will add n to the address of any order terminated with L before 
placing it in position. The directive 00 nK is not placed in the 
memory. Thus if we have 00 11K followed by the orders of Table 5-3 
we will place in the memory the code given in Table 5-l> while the 
code of Table 5-2 would be obtained by using the directive 00 20K. 

5.3 ASSEMBLING OF ORDERS . Since each address ends with 
an alphabet character and since there are always two function digits, 
the Decimal Order Input can distinguish between address digits and 
function digits. The Decimal Order Input uses the fixed storage 


1 


We shall use small letters to represent decimal quantities and 
capital letters to represent sexadecimal quantities. 



locations 0, 1, and 2 as temporary storage. After two function digits 
have been read and shifted to a right-hand order position in location 

-39 

1 the routine converts the decimal address n to binary form n x 2 
and adds it to the function digits. The order pair in location 1 is 
then stored. Next the contents of 1 are shifted left 20 places and 
the next order is formed as before. Again the order pair in location 1 
is stored in the same address as the previous time , but now the correctly 
assembled order pair has been stored. The address of the store order 
which is assembling the program is increased every other time . Thus 
the left-hand 20 digits of location 1 always contain the previous 
order while the next 8 digits of location 1 contain the function digits 
of the order being put in. 

When the address ends in K (a directive) it is stored in 
location 2 and is added to each L terminated address in location 1 
before the appropriate order is stored. 

5.4 DECIMAL ADDRESSES . The address is converted one 
digit at a time as it is read from the tape. This kind of con¬ 
version makes it unnecessary to write non-significant zeros, so 
orders may have the form L5 7L, 40 1021F, 26 L, etc. Indeed, 
we are not restricted to addresses smaller than 1024 and may use 
anything we please. If the address exceeds 4095 it will add to 
the function digits (remember that there are two unused digits 
between the function digits and the address), but even this may 
be made use of as we shall see below in Section 5*6. 


5-5 



5.5 STARTING THE PROGRAM . After the Decimal Order 
Input has placed a program in the memory we must start the pro¬ 
gram. This is done by using the terminating symbol N. The symbol 
N causes the order which it follows to be obeyed. It is used with 

a control transfer order which must appear on the tape as if it were 
a left-hand order . It will never be stored in the memory, other 
than as the right-hand order in location 1 . 

For example, the order 26 pN will be obeyed and will trans¬ 
fer control to the left-hand order at address p. Any unconditional 
transfer order may be used. The order 20 qN will stop the computer, 
and when the computer is started again control will go to the right- 
hand order at address q. The addresses p and q are fixed addresses. 

Control can also be transferred using a relative address, 
for the previous order with its adjusted address is always in the 
left-hand side of location 1. Therefore the two orders 22 rL 
26 IN will cause control to be transferred to address r relative 
to the last directive if 22 rL is the right-hand order of the last 
order pair, i.e., if the phase is correct. 

5.6 INPUT OF DECIMAL FRACTIONS . In section 5.4 it was 

pointed out that an address is formed and added to the function 

39 

digits. This address could be as large as 2 - 1 without getting 

into the sign digit, and any positive integer n could be input as 
-39 

n x 2 by letting the left-hand order and the right-hand function 
digits be zero. For example, the "order pair” OOF OO 2896 F would 
appear as 2896 x 2 J . Hence a 12 digit right-hand address smaller 
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than 2 ^ (about 5.5 x 10 11 ) could then be converted to a decimal 

on i p 

fraction by multiplying it by 2 jy /lO . 

This is done when the terminating symbol J is used. For 
example, the characters OOF 002969 0000 0000J would cause the 

quantity 0.2969 to be placed in the memory. Zeros at the end can¬ 
not be omitted, but preceding zeros can be omitted. Do not omit 
function digits . Remember that this is essentially an integer 
input. The range can be extended by using theleft-hand function 
digits, for they are simply a number to which the right-hand 
address is added after being multiplied by 2^/l0^. The left- 
hand function digits 40, 80, NO represent l/2, -1, -l/2. Thus 
numbers could be input as follows: 


-0.8888 

8888 

8888 as 

80F 

00 1111 

1111 1112J, 

0.7854 

3216 

0000 as 

40F 

00 2854 

3216 0000 J. 


This is not an efficient way to read decimal fractions into the 
machine, but it is convenient for occasional numbers scattered 
through the program. 

5.7 PRE-SET PARAMETERS . MODIFICATION OF ORDERS . The 
remaining unused alphabetic sexadecimal character is S. This is 
used to modify orders by using pre-set parameters. The symbol S 
differs from the other terminating symbols in that it is always 
followed by another single character which may be any of the 13 
sexadecimal characters 3 to L. The termination SD on an address 
causes the content of loqation D = 3* 4, . . . , F, L, to be 
added to the order while it is in a right-hand position as 
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described in Section 5*3 and before it is placed in the memory. 

This is for either right or left-hand orders of completely assembled 

-39 

order pairs. For example, if location J contains 7x2 , the 

order pair L5 20SJ 40 30SJ will be modified to read L5 27 F 40 

37F before being stored. This facility is very convenient for using 
parameters with a program because the program can be written with 
orders of the form L5 S3, L4 S4, and if locations 3 and 4 have 
been previously set the parameters will be added appropriately as 
the program is read into the machine. 

Many examples of pre-set parameters can be found in the 
library. In Illinois Code P6, Single Column Print, the parameter 
OOF OOmF is used to specify the printing of m decimal digits and 
must be placed in location 3 before Code P6 is read. Thus if we 
wish to place Code P6 in locations beginning with 50 and to print 
7 decimal places, the pertinent part of the program tape would 
read 


00 

3K 

Directive 

00 

F 

Place 7 x 2 in 

00 

7F 

location 3 

00 

5 OK 

Directive 

Code 

P 6 

Place Code P6 beginning 


with location 50 . 

5.8 EXAMPLE OF USE OF DECIMAL ORDER INPUT . The decimal 
Order Input has 25 words and occupies locations 999 to 1023 in the 
memory. It is placed in these locations with its own bootstrap 
input which then transfers control to the Decimal Order Input 


5-8 



so that it can take over the control of program input. Let us con' 


sider the following simple example: 

Compute to 10 decimal places the square roots of 
it/lO and e/lO using the Square Root Routine (Code 
Rl)and the Single Column Print routine (Code P6). 

We have the following data: 

(a) Code Rl - 9 words, closed, finds 
square root of argument placed in 
A and places answer in A. 

(b) Code P6 - l4 words, closed, prints 
words in A to m decimal places followed 
by carriage return and line-feed. Ccn- 

-39 

tents of address 3 must contain m x 2 
as Code P6 is input. 

In this program, the code for which is given in Table 
5.1+ on the following page, we have scattered the orders through 
the memory to indicate how directives are used. Notice that the 
arrangement on the tape is arbitrary after the Decimal Order Input 
except that the parameter in address 3 must be in place when Code 
P6 is being input. The parameter is used in word 6 of Code P6, 
this word being 

19 63 S 3 

50 F 

Thus the address of the 19 order was set to 9 (73 = 9 mod 64) for 
counting the number of digits to be printed. 

Notice also that the directive 00 56 OK can be changed 
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MEMORY LOCATION PROGRAM TAPE 


REMARKS 


999 - 1023 

Decimal Order Input 

Routine XI 


00 10K 

Directive 

10 - 18 

Square Root Routine 

Routine R1 


00 3K 

Directive 

3 

00 F 

Parameter 


00 10F 



00 30K 

Directive 

30 - h3 

Print Routine 

Routine P6 


00 50K 

Directive 

50 

00 F 

rt/lO 


00 314159265359J 


51 

00 F 

e/lO 


00 271828182846J 



00 56OK 

Directive 

560 

L5 50F 

jt/lO to A 


50 L 

Link 

561 

26 10F 

To Routine R1 


50 1L 

Link 

562 

26 30F 

To Routine P6 


L5 5LF 

e/lO to A 

563 

22 3L 

Waste Order 


50 3L 

Link 

564 

26 10F 

To Routine R1 


50 4L 

Link 

565 

26 30F 

To Routine P6 


OF F 

Stop 

566 

22 6L 

Waste 


26 L 

Start Program at 



relative location 


H 

=t 

CVJ 

0 sifter stop 


Table 5-4 



Use of Decimal Order Input 
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so that the words following it are placed differently bnt that no 
other change need he made to move these woras^because the program 
was started using a relative address. The waste order at location 
566 is required so that 24 IN will have a left-hand location 
5-9 USE WITH INTERLUDES . RETAINED DIRECTIVE . An 
interlude is a computation performed during the input of a program, 
the input being interrupted and then resumed. 

A tape bearing a library routine may begin with an inter¬ 
lude which is- placed in locations destined eventually to hold the 
routine itself. When the words of the interlude have been read 
control is directed to it, the interlude is executed, and then in¬ 
put is resumed. The next part of the tape carries the routine it¬ 
self which is written over the interlude. The purpose of the inter¬ 
lude is usually to prepare some orders or constants required for 
the routine. For example, a printing routine, where the number of 
digits printed is determined by a preset parameter, may use an 
interlude to compute the roundoff constant. 

Input is resumed after the interlude by transferring con¬ 
trol to the left side of location 999 (3F7 sexadecimal). Either 
the first word on the tape must contain the needed directive 
00 mK or Q must contain m x 2 

If, upon resuming input, it is desired to retain the last 
used directive , control should be transferred to the right side of 
1014 with m x 2 ^ in A. The next words on the tape will be placed 
in m, m+1, . . . , retaining the previous relative address. The 
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Interlude must not uae location 2. 


5-10 STOPPING THE TAPE , The order 20 1019N on the 

tape will stop the computer and will have no other effect. Upon 
being started the tape will continue being read from where it stopped. 

5-11 PLACING THE DECIMAL ORDER INPUT . BOOTSTRAPS . 

Up to this point we have not said how the Decimal Order Input is 
itself put into the Illiac. It occupies locations 999 to 1023, 
the last 25 positions of the electrostatic memory, and it is placed 
in the memory with a bootstrap input routine . 

With panel switches we can clear the control counter to zero and 
place the sexadecimal equivalent of the order pair 80 40F 40 F in the 

order register. The code in Table 5.4, which must be written with sexadeci- 
mal addresses, will then place the Decimal Order Input in the memory. 

80 028 

40 001 

80 028 

40 002 

19 026 

26 000 

80 028 V 

4o(ooo) d) 

L4 001 
40 001 

80 028 

4o(3F6) 

Table 5.4 

Tape for Bootstrap Input Routine 

s 

Parentheses are often placed around addresses which change during 
the course of a program. 
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This bootstrap actually places the 3-word routine (shown in both 
sexadecimal and decimal forms) of Table 5-5 in the memory and starts 
it 



l4 

001 


L4 

IF 

0 

4o 

001 

0 

4o 

IF 


8 o 

028 


80 

4 OF 

1 

40 

(3F7) 

1 

4o 

(999)F 


19 

026 


19 

38F 

2 

26 

000 

2 

26 

F 


Table 5-5 

Memory Contents for Bootstrap Input Routine 

at location 1. Clearly it will take the next words on the tape 
and start putting them at location 999F. To stop the input we 
place the order pair 22 3LS 00 001 (This is 22 1019F 00 IF) 

on the tape so that it gets put into location 0. The control 
will be transferred to the right-hand side of location 1019F and 
the Decimal Order Input will be started. 

This bootstrap input may be used with any code and it 
or something like it must be used to input programs whenever the 
Decimal Order Input has been overwritten. 

The term bootstrap start is often used for tapes which 
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are started by setting the order register to 80 40F 40 F and 

the control counter to zero. 
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CHAPTER 6 


SCALING 


With the convention adopted for the Illiac (See Chapter 2), 
only numbers which lie in the range -1 < x < 1 can be held in the 
registers. Since most problems require numbers outside this range, 
some scaling process is usually needed to fit a problem to the 
machine. It is necessary that each number at every stage of a 
calculation lie within the capacity of the machine. The organization 
required to assure this is sometimes trivial, but in many instances 

it is the very essence of the problem. 

6.1 SCALING BY SHIFTING . Although the number 2 lies out¬ 
side the range of Illiac numbers, we can multiply and divide numbers 
by powers of 2 by shifting. Thus the left shift order 00 IGF will 
cause AQ to be multiplied by 2 10 = 1024. Similarly, the right shift 
order 10 9F will divide AQ by 2 9 = 512. A knowledge of the use of 
the shift orders is essential to an understanding of scaling. 

6.2 NUMBERS WITH THE BINARY POINT SHIFTED . Let us con¬ 
sider the problem of computing with numbers in which the binary point 
has been moved 10 places to the right of its Illiac position. We 

-29 

are then dealing with numbers in the range - 1024 < y < 1024 - 2 

Let these numbers be designated by N^. Then we have N^q(iii) - 2 N(m). 

We can formulate rules for doing arithmetic with the numoers 
N Addition and subtraction are simple. If N(q) = N(m) + N(n), 
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then 


2 10 N(q) = 2 10 N(m) + 2 10 N(n) 

and 

N 10 (q) = N 10 (m) + N 1Q (n) . 

Thus the Illiac addition rules hold. 

Multiplication requires a shift to the left of 10 places. 

We want 

N 10 (q) = N 10 (m) x N 1Q (n) . 

Thus we require 

2 10 N(q) = 2 10 N(ra) x 2 10 N(n) 

and 

N(q) = 2 10 N(m) x N(n) . 

We may consider several simple routines to carry out the 
multiplication. The method given by (b) is probably the best. 

(a) The shortest method merely multiplies and shifts . 

Notice that a 75 order must be used rather than a 7^ order because 

the least significant digit of N 1Q (.p) is now 2 -29 and not 2 _39 . 

-30 

A bias of -2 is introduced by the absence of a roundoff. The 
program is given in Table 6.1. 

(b) This method rounds off by adding 2 ~^° to N, ,(p), 

10 

giving an unbiased result. The program is given in Table 6.2. 
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50 mF N(ra) to Q 

75 nF N(m) x N(n) to AQ 

00 10F 2 10 N(m) x N(n) 

40 pF to p. 

Table 6.1 

Multiplication with Binary Point Shifted 

19 10F 2” 11 to A 

50 mF N(m) to Q 

74 nF N(m) x N(n) + 2"to AQ 

00 10F 2 10 N(m) x N(n) + 2"^° to A 

40 pF to p. 

Table 6.2 

Unbiased Multiplication with 
Binary Point Shifted 

(c) A roundoff similar to that of division is obtained 
with the program given in Table 6.3. 


50 mF 

N(m) to Q 

75 nF 

N(m) x N(n) to AQ 

00 9F 

2^ N(m) x N(n) to AQ 
-1 „ 

50 tF 

2 Q = 1 

00 IF 

-79 

This makes 2 A = 

N(A) = 2 10 N(m) x N(: 

40 qF 

to q. 


Table 6.3 

Division-Type Roundoff in Multiplication 
with Binary Point Shifted 
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In Table 6,3 the order 50 tF is used simply because the word at 
location t has a 1 in the proper place. 

In division we also need an extra shift to restore the 
quotient to the proper range. But here the shift precedes the divide 
order and no special arrangements for roundoff are necessary. We 
want 

N 10 (q) = N 10 (m)/N l0 (n). 

Hence 

2 10 N(q) = N(m)/N(n), 

and 

N(q) = 2' 10 N(m)/N(n). 

The following program in Table 6.4 will carry out the 
required operations: 

L5 mF N(m) to A 

10 10F 2 -10 N(rn) 

66 nF 2 -10 N(m)/N(n) 

S5 F to A 

40 qF to q. 

Table 6.4 

Division with Binary Point Shifted 


It can be noted that by using such relations as 


10 


P 
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where 10 P < 2 n we can use decimal scaling although it will he slower 
and clumsier to handle because of the factors 10 5 / 2 n . 

6.3 SCALING A FULL PROBLEM . There are conceptually two 

ways in which we can approach the scaling of a problem; both give 
the same program. 

(a) We can alter the problem using such 
substitutions as ^ = lOOx or p = 32p, 
so that' the modified problem has all its 
variables less than one but retaining full 
significance. 

(b) We can use scaled numbers inside the machine 
to represent the variables. That is, we can 
use x/100 instead of x and p /32 instead of p. 

The final result must be such that the variables lie 
within machine range and retain sufficient accuracy. Constants 
greater than one can be represented by numbers less than one in 
conjunction with scaling factors. For example, multiplication by 
5.63 can be done by multiplying by 5-63/8 and shifting left 3 places. 

Let us consider the following simple problem: 

Example 1. Program the Illiac to compute the quantities 

x 

e 

y -- 

1 + X 2 

at intervals of 0.01 from x = 0 to x = 6. 

We note that we must scale x. Since the largest value of 
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x is 6 , let us use x /8 inside the computer because this will minimize 

the loss of significant figures . 

What scaling factor is needed for y? A rough estimate shows 

6 

that y does not exceed e / (l + 36 ) ~ 11 , so that we shall compute 
and store y/l 6 . 

Library Routine S 2, the exponential routine, will give 
valid results only if x is negative. We therefore write 

x /8 ^ „ -1 + x /8 

so that 

e X = (e X / 8 ) 8 . 


X 

Now y/l6 = --- X 

16 x 64 


1 

» 

l/64 + x^/64 


and 


e 

16 x 64 


|[(e 1 + x / 8 x e/h -) 2 x 2] 2 X 2 j 2 . 


We shall store the constants e/4 and l/64. Each part of the computa¬ 
tion is within machine capacity. We proceed as follows: 

(a) Square x/8 and add l/64, 

(b) Call the result R and store it, 

(c) Evaluate e * 1 + X ^ 8 with Library Routine S 2, 

(d ) Multiply by e/4 and square, 

(e) Double and square. 
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(f) Double and square again. Call the result P. 

(g) Form y/l6 by dividing P by P>. 

Accuracy in the Result . Let us now consider how much 
accuracy we obtain. If x/8 is near unity, numbers remain large 
during the calculation and we do not lose significance by sub¬ 
tracting nearly equal numbers or by other ill-conditioning. 

However, when x/8 is small we notice that we form e X / (l6 x 64) 
quite accurately and then divide by R ~ 1 /64 which is about equiva¬ 
lent to a left shift of 6 places and losed 6 binary digits on the 
right. We can prevent this if we perform the division while the 
double-length P is in AQ. 

Accuracy in the Argument . We should also consider the 
accuracy of the argument x/8. It can be formed either by successively 
adding the Increment 0.01/8 o 0.00125 or by counting and multiplying. 

The adding method is not so good because the quantity 0.00125 will 
not be stored exactly and the accumulated roundoff error obtained 
by the time x is 6 (which requires 600 additions ) may be troublesome. 

The counting method avoids this trouble. And x/8 is 
n x 0.01/8. Hence we have 

x/8 = n x 0.01/8 

= (n x 2 ^) x 0.64 x 2^° 

If we store 0.64 and count to get. n x 2 - ^, we can use the following 
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set of orders to get x /8 (where we have used arguments in place of 
addresses): 

50 n x 

75 0.64 

00 30 

4 o x/8 . 

The program for the entire calculation of y/l 6 is given 
in Table 6.3 on page 6-9 where we have again used arguments instead 
of addresses. 

Example 2: Solve the pair of simultaneous equations 


ax + by + c =0 
dx + ey + f = 0 

where the coefficients are in absolute values less than l /2 and 
where the answers are known to lie within machine range. Retain 
as much accuracy as is reasonably possible. 

We can distinguish two cases: 

(a) If |a| > |d| , then 


y = 


dc/a - f 
bd/a - e 


x 


c + by 
a 
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ORDER 

ARGUMENT 

ORDER 

ARGUMENT 

F5 

n X 2 *"^ 9 

50 

e -l + x /8 

p 

1+0 

(n + 1) X 2 ~ 59 

7J 

e X 2 

50 

(n + 1 ) X 2 -59 

40 

A 8 / 4 

75 

.64 

50 

e x / 8 / 4 

00 

50 

75 

00 

0 

l+o 

x /8 

00 

1 

50 

x /8 

40 

e x /“ / 8 

7J 

x /8 

50 

A 4 / 8 

Ll+ 

1/64 

75 

e x / 4 / 8 

1+0 

l/64 + x^/64 

00 

1 

L5 

x /8 

40 

e ^ 2 / 32 

l4 

-1 

50 

e x/2 / 32 

50 

(link address) 

75 

OJ 

N"\ 

OJ 

0 

26 

to S 2 

66 

l/64 + x^/64 

40 

-1 + x /8 
e ' 

S5 

(any address) 



4o 

y/l 6 


Table 6.5 

2 

Calculation of e X / (l + x ) 
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(b) If |d | > |a| , then 


y = _ c - fa / d 
b - ae/d 

x = . —L± e y 

d 

In this program we shall follow the conventions of 
the Decimal Order Input (See Chapter 5)- Notice that divisions are 
always made using a full Jd digit dividend to retain as much assuracy 
as possible. The program treats the two cases (a) and (b) separately, 
distinguishing with the 36 order at 1L. Locations 0 and 1 in the 
memory are used as temporaty storage. The program is given in Table 
6.6 on pages 6-12 and 6-13. 

6.4 ADJUSTABLE SCALING FACTORS . It is not always 
possible to arrange a program so that a single scaling factor can 
be used throughout the calculation. Then it is necessary to make 
tests at appropriate places to discover when variables are becoming 
too large or too small and to make proper adjustments in the scaling 
factors. For many problems it is advantageous to have the variable 
less than l/2. Then two numbers can be added or multiplied without 
exceeding capacity. By using the LL n order we place’l/2 + |N(n)| 
in A. Hence A is positive if |N(n)| < l/2 and A is negative if 
|N(n)|>l/2. 

6.5 CONTINUOUS SCALING . FLOATING POINT ROUTINES . For 
calculations in which continual tests are required to maintain accuracy 
floating point routines (See also Section 4.7) may be used. These 
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routines represent numbers as y = a x 10* 3 and store a and b. Thus 

they can represent accurately the numbers in some large range such 

-65 63 

as, for example, 10 < y < 10 where y has 30 significant binary 

digits. There are two such routines in the Illinois Program Library. 

The first. Library Routine A 1 is as described above. The second. 

Library Routine A 4 is a multiple precision floating point program 

-153 153 

in which numbers lies in the range 10 < y < 10 with y 

having 20 significant decimal places. 

Floating point routines are slow because numbers are 
scaled at each step of the calculation. Certain conveniences have 
been programmed in, however, and these to some extent compensate for 
the extra time required and also simplify the programming. 
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0 

L7 22L 

L2 25L 

|a| - |d| 

1 

36 13L 

50 27L 

|a | < |d | } f to Q 

2 

75 22L 

66 25 L 

fa/d 

3 

SI F 

L4 24L 


4 

40 F 

50 22L 

-fa/d + c to 0 

5 

75 26L 

66 25L 

a e/d 

6 

SI F 

L4 23L 


7 

40 IF 

L5 F 

-ae/d + b to 1 

8 

66 IF 

SI F 

(-fa/d + c) / (-ae/d 

9 

40 29L 

50 29L 

y to 2% 

10 

75 26l 

L4 27L 

(ey + f) 

11 

66 25L 

SI F 

-(ey + f) / d = x 

12 

40 28 l 

OF F 

x to 28L 

STOP 

13 

50 23L 

75 25L 

|a| > |d| } b to Q 

14 

66 22L 
•S5 F 

bd/a 


Table 6.6 

Solution of Two Simultaneous Equations 
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15 

L0 

26L 

bd/a - e to 0 

40 

F 

16 

50 

25L 

d to Q 


75 

24L 


17 

66 

22L 

dc/a 


si 

F 


18 

L4 

27 L 

(-dc/a + f) / (bd/i 


66 

F 

19 

S5 

F 



40 

29L 

y to 29L 

20 

75 

23 L 

by 


L4 

24L 


21 

66 

22L 

(by +c)/a = -x 


22 

11L 

Control to 11L 

22 



a 

23 



b 

24 



c 

25 



d 

26 



e 

27 



f 

28 



X 

29 



y 


Table 6.6 (Continued) 
Solution of Teo Simultaneous Equations 
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CHAPTER 7 


MACHINE METHODS AUTO CODING TRICKS 

There are usually a number of special techniques which 
can be used on any particular digital computer and which will simplify 
programming. Some of these techniques are applicable on many 
different computers but usually, as is the case in those which follow, 
they result from particular orders or combinations of orders which 
are peculiar to an individual machine. This chapter is concerned 
with a number of unrelated sections having to do with operations 
which frequently arise in programming. 

7 .I THE SUMMATION OF PRODUCTS . We often need to form 
sums of products, and on the Illiac this cannot be directly done in 
the accumulator. The accuracy can often be enhanced by performing 
a summation either exactly or with only°Dne round-off error. This 
is comparatively easy to do using the 7^ order. All that we need 
to do is to place the least significant half of the partially summed 
products into the accumulator before performing the 74 order. In 
fact, this can usually be done by an 85 order because the quotient 

register will usually hold the last single half of a summed product. 

-39 

Then, since 74 n gives N(n) Q + 2 A, we obtain the double-length 
product in AQ. Of course, the most significant part of the pre¬ 
viously s umm ed products needs to be added using an L4 order. 

Using similar schemes we can also arrange to add or 
subtract products with double-length accuracy in a program. As a 
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first example we shall sum 50 double-length products t assuming we 
do not exceed capacity. 


Example 1 : Place the rounded sum 

49 

Z N (100 + i) N(l50 + i) 
i=0 


in location 0. The program is given in Table 7.1 


m 

4l F 

26 1 L 

Clear location 0 for sum. 
Wasted order 

m+l 

50 8 L 

L5 11 L 

Put 1/2 in Q 

Set i » 0 

m +2 

40 3 L 

S5 F 

Round off (See Section 7*2) 

®T5 

50 ( )F 

74 ( )F 

N (100 + i) N (150 + i) 

m+4 

L4 F 

40 F 


m+5 

L5 9L 

L4 3L 

Increase i by 1 

m +6 

40 3L 

L0 10L 

Test for i > 150 

m+7 

32 2L 

OF F 

Re-enter loop 

Stop 

m +8 

40 F 

00 F 

Roundoff constant *» 1/2 

m+9 

00 IF 

00 IF 

Increment 

m +10 

JO 150F 

74 200F 

End constant 

m +11 

50 100F 

74 150F 

Starting constant 


Table 7.1 

Program for Example 1 
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7.2 REVERSING THE CONTROL TRANSFER . FURTHER DISCUSSION OF 
EXAMPLE 1 . There is a second coding trick in Example 1 . The end 
constant, instead of being 50 150F 7^ 200F has had -1 added to it, 

making the first order JO 150F. The effect is to reverse the sense 
of the following 32 order (in location m + 7 ) so that we transfer 
control to re-enter the repetitive loop. If this had not been done 
a 22 order following the 32 order would have been required and the 
technique is equivalent to having a conditional transfer order act 
when the accumulator is negative. 

Another coding trick might have been used to save a full 
word. The left-hand order at m + 1 puts l/2 in Q, so that it may be 
used to round off on the first step of the sum, this being the sole 
roundoff. Instead of storing l/2 in m + 8 we could have used the 
order pair in m + 2 as the roundoff constant. This order pair is 
l/ 2 , plus at most 2 ~^ and would serve quite well. 

As a second example we consider a summation of two products 
with a single roundoff. 

Example 2: Given 

cos 0 in 10 , 
sin © in 11 , 
x in 12 , 
y in 13> 

place the rounded quantity (x cos 9 + y sin ©) in 
location 20 and the rounded quantity (-x sin © + y cos ©) 
in location 21. The program is given in Table 7*2. 
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m 

50 

10F 


7J 

12F 

m +1 

4o 

F 


S5 

F 

m +2 

50 

11 F 


74 

13F 

m+3 

L4 

F 


4o 

20F 

m+4 

50 

11F 


79 

12F 

m+5 

4o 

F 


S5 

F 

m +6 

50 

10F 


74 

13F 

m+7 

L4 

F 


4o 

21F 


-4o 

x cos 6 + 2 

Most significant half to location 0 
Least significant half to A 

y sin 0 + 2 (l.s. half of x cos 0 + 2 ”^) 

-40 

Add most significant half of. x cos 0 + 2 
Store x cos © + y sin © + 2"^ 

-40 

-x sin © + 2 


y cos © + 2 ~^ (l.s. half of -x cos © + 2 -HU ) 


-40, 


Store -x sin 0 + y cos © + 2 


-40 


Table 7.2 

Program for Example 2. 


7.3 BINARY SWITCHES . It is sometimes necessary to do two 
different operations alternately. This can be done by changing the 
sign of a number each time we pass it so that it will be alternately 
positive and negative. Usually it is not necessary to use a special 
number for this because some number or order pair in the rest of the 
program may be used. To accomplish the switchman LI order followed 
by a 40 order is used to change the sign of the number and put it 
back in its location with sign changed. A conditional transfer order 
may then be used to decide which of two sequences will be performed. 

A variation is the requirement that an order (or order pair) 
take on two different values alternately. This can be accomplished 
by using the identities 
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and 


b a (b + a) - a, 
a = (b + a) - b. 

Thus, if the current value of an order (or order pair) is subtracted 
from the sum,the other value is obtained. 


Example 3: Arrange a program to alter the address 
of the left-hand order at (m+ 2 ) so that it takes on 
the values 0 and 5 alternately. A program for this 
is given in Table 7*3* 


m any order 

L5 pF 


Put sum of orders in A 


m+1 LO 2L 
40 2 L 


Form alternate address 

Store alternate address at (m+2) 


ra +2 L5 ( 0 )F 

L4 12F Normal program.order 

p FK 6 F Sum of L5 OF L4 12F 

F 8 24F and L 5 5F L4 12F 


Table 7-3 
Binary Switch 


In Table 7.3 only a single address is taking on alternate 
values and it is also possible to carry out the switch using the 
program of Table 7.4. In Table 7-^ the 50 5F order at location 
m is provided for the switch. If some order which needed address 
5F could be used here, we should be very well off indeed. 
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m 50 5 F 

L5 mF 

m+1 LO (m+2)F 

46 (m+2)F 

ra+2 L5 (O)F 

L4 (12)F 

Table 7*4 
Binary Switch 

7.4 TESTS FOR 0 AND - 1 . In order to test for a particular 
number value held inside the machine it is generally necessary to 

use two tests. However, the numbers 0 and -1 can be tested for using 

absolute value orders and a single test. In machine l anguage 0 is 
the onl y number whose, .ne gative abseil ute value is positive , and -1 

is the only number whose positive absolute value is negative. 

/—'----- 

Thus we can test for 0 using an L3 order followed by a conditional 

transfer order, and we can test for -1 using an L 7 order. Similarly, 

-39 -39 

we can test for -2 J and 1 - 2 by using F 3 and F 7 orders, res¬ 
pectively. 

Example 4: Transfer control to location 200 if A 

is zero but transfer control to location 300 if A 

is non-zero. Two ways to do this are given in Tables 
7-5 and 7 . 6 . The program of Table 7.5 has only two 
words. The program of Table 7.6 has four words but 
will be faster than the other if A is usually negative. 
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Moreover, the program of Table 7.6 can be used 
to transfer control to any of three l ocaLbn s de¬ 
pending upon whether A is positive, negative or zero. 


m 

40 

F 

A to 0 



L3 

F 

- |n(o) | to 

A 

m+1 

36 

2 OOF 

To 200 if A 

> 0 


26 

300F 

To 300 if A 

i 0 


Table 7-5 
Testing for Zero 


m 

36 

26 

1L 

300F 

To 300 if A < 0 


m+1 

LO 

pF 

A " 2 " 39 2TQ 



36 

300F 

To 300 if A - 2'^ > 0, 

i.e., if 

m+2 

26 

2 OOF 

To 200 otherwise, i.e.. 

if A = 0 

P 

00 

F 

Constant 2 - ^ 



00 

IF 




Table J.6 
Testing for Sign 


7-5 USE OF ORDERS AMD ADDRESSES AS CONSTANTS . SV and 
KV orders do not use their addresses, so these addresses can often 
be used for other purposes. For instance, they may be used to store 
a starting address taken by a cycling order or an increment which 
is used to change an address. In such cases we naturally use 42 or 
46 orders to make certain that the function digits do not become 
altered. 
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In the following example the address of a K5 order is 


used as a counter. 

Example 5s Given the positive number a in A, 
write a closed subroutine which will furnish 
the positive integer m such that l/2 < 2 m a < 1. 
The program is given in Table 7*7* 


m 

40 

F 

Store a at location 3 


K5 

F 

Form link 

m+1 

42 

4l 

Plant link 


45 

L 

Clear counter 

m+2 

L5 

F 



00 

IF 

Shift a 

m+5 

40 

F 

■Q 


36 

5L 

Test to see if 2a > 1 

m+4 

L4 

L 

Counter to A 


22 

( )F 

Link 

m+5 

F5 

L 

Count 


4o 

L 


m+6 

26 

2L 

Re-enter loop 


00 

F 

Waste order 


Table 7-7 

Address Use in K5 Order 


7.6 RESETTING AND STARTING OF CYCLES OF ORDERS- 


In 


many cases we have cycles of orders, some of which are being modified 
by the same increment. In such cases the variable addresses can all 



be modified by modifying one order and then deriving the other orders 
(or addresses) by adding the constant difference between the variable 
orders. When this is done it is economical to use the same orders 
to set these addresses when the cycle is begun. The following example 
with the program given in Table 7.8 illustrates this: 


Example 6: 

For 

i = 0, 1, . . . ,99 place in 

location 200 + i 

the sum N(lO) + N(l00 + i) + N(200 + i). 

m 

1*5 

7L 



26 

4l 

Set i = 0 

m+1 

1*5 

10F 

Form N(lO) + N(l00 + i) + N(200 + i) 


l4 

F 

and put in 200 + i 

m+2 

L4 

F 



4o 

F 


m+3 

L5 

8L 



L4 

2L 

Increase i in m + 2 

itW-4 

4o 

2L 

Increase i in m + 1 and test for end 


L0 

9L 


m+5 

42 

1L 



36 

1L 


m+6 

OF 

F 

Stop 


00 

F 

Waste order 

m+7 

l4 

200F 

Starting constant 


4o 

2 OOF 


m+8 

00 

IF 

Increment 


00 

IF 


m+9 

74 

300F 

End constant 


00 

100F 

Constant to change i in 2L 


Table 7.8 

Resetting of Addresses 
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In Example 6 the addresses in location m + 2 have been 
changed and then that in m + 1 has been obtained from one of them. 
Notice that the end test has been combined with the second address 
change and that the end test constant is not LU but 7^ so that the 
36 order at location m + 5 will cause re-entry to the repetitive loop. 

When two orders have to be varied in a single cycle, it is 
advantageous to let these form a single order pair as in Example 6 
so that they can both be modified simultaneously by the same orders. 
This arrangement is not always possible and the second best arrange¬ 
ment is to place the variable orders on the same side of their res¬ 
pective order pairs so that the orders required to modify them will 
be as simple as possible. 

For very simple operations it is sometimes advantageous to 
do three or four operations in a single cycle. This saves time 
although the advantage is bought at the expense of more orders. This 
is illustrated in Example 7 • 

Example 7 : Add the numbers in memory locations 
10 to l 4 , putting the sum in location 15* It is 
simpler and faster, both in coding and in machine 
operation, to use the program given at the top of 
page 7-11 than to write a repetitive code which counts. 
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1*5 

L4 


10F 

11F 


1.4 12F 

L4 13F 

L4 i4f 
4o 15F 

7-7 USE OF THE QUOTIENT REGISTER FOR INTERCHANGES . In 
many programs we wish to replace the number in a certain storage 
location and yet use the value that is there to continue with the cal¬ 
culation. In such cases the old value can be placed in the quotient 
register before the new value replaces it in the memory. Thus, the 
old value is available in the quotient register for further computation. 

Example 8: Store A in location 10, but use the old 
N(l0) as a dividend to form N(l0) / N(ll). The program 
is then 


50 10F 

40 10F 

S5 F 
66 11F 

7 *8 TESTJJlg IF NUMBERS -ABE GREATER THAN ORE-HALF . When 
scaling numbers it is very often necessary to test when numbers are 
larger in magnitude than one-half. This can easily be ddne with the 
appropriate L or S order. For example, the order LL n will cause the 
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accumulator to be negative if the magnitude of N(n) is greater than 
or equal to l/2. 

Example 9: If JN(10) | > l/2 replace it by half 

its value. The program is given in Table 7-9* 

m LL 10F 1/2 + |N(lO)| to A 

32 (m+2 )F 

m+1 L5 10F 

10 IF 

m+2 40 10F 

Table 7-9 

Scaling by Testing for One-half 

7-9 CONVERGENCE CRITERIA . When iterations or repetitive 

calculations are carried out we frequently want to stop when we have 
achieved the maximum accuracy. In some cases it is difficult to 
specify in advance the tolerances which can be used as end criteria 
because we have to compromise between achieving the greatest accuracy 
and yet assuring that we terminate the processes (i.e., don't loop). 

In such cases it is worthwhile to use more complicated 
criteria which will give us maximum accuracy but which will not loop. 
One such criterion is to terminate the process if e n > e n-1 where e^^ 
is some positive number which tends to zero as the process converges 
with increasing n. This criterion will terminate the process only 
when either e ± is the same for two successive iterations or when the 
roundoff error has actually caused it to increase. 
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7.10 MARKING . It is often possible to use marking 
techniques instead of the more usual counting processes. The 
simplest illustration of such a technique is Library Routine N3 
where a sequence of numbers read from the tape is automatically 
terminated by the mark N. In this code instead of counting up to 
some predetermined number we test each character as it is read from 
the tape until it is N. 

A binary digit is sometimes shifted as a marker; this is 
illustrated in the next example. 

Example 10: Using a print routine stored at locations 
beginning with 50 print the 7 numbers in locations 10 -l 6 . 
The program is given in Table 7 .10. 


0 19 

7F 

0 

2~ to IF. This is the marker 

4o 

IP 


1 L5 

IP 


00 

IF 

Advance marker by shifting 

2 40 

IF 


32 

3L 

Test for end 

3 OP 

F 

Stop 

L5 

10F 


4 22 

4l 

Waste order 

50 

4l 


5 26 

50F 

Enter print rojutine 

F5 

31 * 

6 40 

3L 

Increase address of number to be printed 

26 

1L 

re-enter loop 


Table 7-10 


Use 

of a Marker 
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In the program of Table 7-10 the marker ia shifted into the sign 
digit to indicate the end of the repetitive process. 

This technique can be used in a similar way for completely 
internal programs. For example, when we are dealing with a group 
of numbers in the memory we may arrange that the storage location 
following the group contains some unique number such as 0 or -1. 

Then the code merely has to test for the presence of one of these 
numbers rather than for a predetermined count. 

7.11 REMAINDER IN INTEGER DIVISION . In the general case 
it is difficult to compute the remainder from the residue that is 
left in the accumulator after a division. However, if we are dealing 
with positive integers less than 2 ^ in magnitude we can do this 
quite readily. We place twice the dividend integer in the quotient 
register, clear the accumulator and divide by the divisor integer. 

The accumulator then contains twice the integer remainder and the 

quotient contains twice the integer quotient. We store an integer 

-39 

m as m x 2 

Example 11: Divide the positive integer in location 
20 by 10. Place the quotient in location 11 and the 
remainder in location 12. The program is given in 
Table 7.11. 
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m 51 fiOF m x 2 -59 to Q 

-38 

00 IF m x 2 J to AQ 

m+1 66 4l (in x 2 - ^) / (10 x2^ 9 )x2^ = ra/lO 

10 IF 

m+2 40 12F Store remainder 

S5 F 

m+3 40 11F Store quotient 

26 pF 

m+4 00 F 

00 10F 

Table 7-U 

Remainder in Integer Division 

7.12 BINARY CHOPPING . This is the method of repeated 

subdivision of an interval. It is easy to code although it may be 

wasteful of memory space and it is slow because it will usually take 

the full 59 steps to go from an interval of length unity down to 
-39 

one of length 2 

If we use binary chopping to find the zero of a function, 

we proceed as follows. We choose bounds for the zero, perhaps -1 as 

a lower bound and +1 as an upper bound. Then we bisect the interval 

and compute the function at the midpoint. Depending upon whether the 

sign of the function at the midpoint agrees with the sign of the lower 

or upper bound, we substitute the midpoint for the appropriate bound. 

-39 

After 39 steps the difference between the two bounds will be 2 
and the zero will be determined. 
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The code will be simpler if the signs of the upper and 
lower bounds are known so that a comparison with the sign of the 
midpoint is not needed at each step. This is the case in Example 
12 . 

Example 12; Find the square root of a = O. 269 L 3 

by using a binary chopping technique on the function 
2 

a - a. 
n 

Here we choose initial upper and lower bounds of 

1 - 2 "^ and 0. Instead of counting 39 steps, 

2 -37 

a test has been included so that if a - a < 2 

n 

the code stops. The program is given in Table 7 *12. 

7.13 EVALUATION OF POLYNOMIALS . Polynomials are best 
evaluated by use of a recurrence solution. Given the polynomial 

f( x ) - Z a x - V s + a x x + • • • + Vi x + a n ’ 
i =0 


we can express it in the form 



0 

L5 

12L 



10 

IF 


1 

4o 

F 

a/2 to locatii 


L5 

13 L 


2 

10 

IF 

(a + a)/2 = a. 


L4 

F 


3 

40 

i4l 



50 

14L 


4 

7J 

14l 



L0 

11L 


5 

40 

15L 

a ^ - a to 15 ' 

y-> 


32 

7L 

11 

6 

L5 

14L 

Change a 


4o 

13L 


7 

22 

8l 



L5 

14L 


8 

4o 

12L 

Change a 


L7 

15L 


9 

L0 

i6l 

2 

Test for a 


36 

L 

n 

10 

OF 

F 

Stop 


00 

F 


11 

00 

F 

a = 


00 

269430000000J 


12 

7L 

4095F 

a = 1 - 2“ 59 


LL 

4095F 


13 

00 

F 

a = 0 


00 

F 


14 

00 

F 

a 


00 

F 

n 

15 

00 

F 

2 

a - a 


00 

F 

n 

l6 

00 

F 

-*37 

2 


00 

4f 



- a < 2 


-37 


Table 7.12 

Square Root by Binary Chopping 
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If, as is often the case, the coefficients a i are the 
quantities N(m+i) the recurrence relation becomes 


= x + N(m+i) , 



S 


n +1 


f(x) 


Example 13 shows how a polynomial may be evaluated. In practice 
the summation of products should be done with a 74 order as in 
Example 1, but we do not wish to obscure the general idea here. 


Example 13t Given x in location 49 and coefficients 
in locations 50 through 67 evaluate the polynomial 

£ x 7-1 N(50 + i) . 
i =0 

The program is given in Table 7.13. 

7 *14 SHIFTING A LEFT WITHOUT SHIFTING Q . The contents 
of A may be shifted to the left (doubled) by adding A to itself. 

This is slower than using a shift order but it has the advantage 
that the Q register is not altered (as it would be by a shift order). 
The order pair 40 nF L4 nF will shift the contents of A one 
place to the left for each use. 
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p 

4l 

F 


50 

49 F 

P+1 

7J 

F 


L4 

50F 

p +2 

40 

F 


F5 

1L 

p+3 

40 

1 L 


L0 

5L 

p+4 

32 

L 


OF 

F 

P+5 

LJ 

F 


L4 

68 f 


Table 
Evaluation of 



x 


x S i + N(50 + i) 


Increase i 

Test for end 


End test constant 


7-13 

a Polynomial 


7-19 



CHAPTER 8 


CHECKING METHODS 

To obtain a correctly working program from a written one 
the coder must find and remove all of the coding errors. To the 
person unfamiliar with coding it might seem that a careful inspection 
of a program before it is run on the Illiac would remove all of the 
coding blunders but unfortunately this is not true. It is necessary 
to check the program on the machine to remove the remaining blunders 
from a code, and it is unusual to find all blunders in one run on 
the machine. As a general rule it will be necessary to have several 
checking runs on the machine before a code is correct. 

The number of coding errors in a program depends upon a 
number of factors, including complexity, length of untested program, 
programmer, and the care with which the code was inspected. It is 
possible to remove most coding blunders from a program by carefully 
inspecting it and it is best to obtain a person other than the coder 
to inspect the code. In many cases this cannot be done, and then the 
code should be set aside for as long as possible before being inspected. 
This to some extent prevents a grooved mind from missing the flaws in 
the program, for the coder remembers the general nature of the program 
but not the details where most errors occur. 

8.1 COMMON BLUNDERS . A great many of the coding mistakes 
made by programmers are familiar ones which are made over and over again. 
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The list given in Table 8*1 is one of common errors. It should be 
pointed out, however, that this list is not complete, and that 
programmers develop their own idiosyncrasies and should be on guard 
for their own pet blunders. 

8 ‘ 2 THE DISCOVERY OF ERRORS„ Errors are found by running 
the program and comparing the actual performance with the designed 
performance. That is, to find an error it is necessary to obtain 
information about the way a program runs. It is almost useless to 
have a program run through a long calculation and then print out the 
result. If the result is wrong, no information is usually available 
to explain why. Therefore, when checking a program we have to print 
out more information about the intermediate results than is required 
in the actual running of the program. There are many ways in which 

this can be done and the remainder of this chapter describes some of 
them. 

Because we have to compare the intermediate stages of the 
calculation with those estimated by other means, it behooves us to 
choose parameters and variables so that the initial calculation is 
as simple as possible. Further trials will probably have to be run 
with more complicated calculations and with values designed to test 
the special cases and boundary limitations. Simple blunders should 
be detected and removed before an attempt is made to look for more 
subtle errors. Localization of mysterious errors can be effected 
hy continually printing more and more information from smaller and 
smaller parts of the program until the error is found. 
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1. L5 orders used instead of L4 orders 

2. S5 orders omitted after divisions. 

3. Orders terminated by L instead of F and vice-versa. 

4. The renumbering of a code not complete after a 
modification has been made. 

5. Rounded multiplication used when dealing with 
integers. 

6. Control transfers to the wrong address or wrong 
order of an order pair. 

7 . Accumulating storage registers not cleared before 
a cycle of orders is entered. 

8. The end condition for a cycle of orders not correct. 

9. Allowing the temporary storage of a subroutine to 
erase useful data. 

10. Using a 46 order instead of a 42 order and vice-versa. 

11., Omitting directives and starting orders on the 
program tape. 

12.. Incorrectly remembering the specifications of a 
subroutine. 

13» Forgetting to reset addresses when coming back to 
a cycle of orders. 

14. Making corrections incorrectly. 

15. Using the same relative addresses on correction 
words although the preceding directive is different 
from that of the program. 

16 . Overlooking the digits shifted from the quotient 
register to the accumulator on a left shift. 

17 . Attempting to convert fractions greater than one- 
half by using the J terminating symbol. 

Table 8.1 
Typical Blunders 
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8.3 SOURCES OF INFORMATION . 

The Punch . The punch is the most effective way by which 
data can be extracted from the Illiac. However, there are other ways 
to obtain information, particularly when a program does not run far 
enough to punch any data at all. 

The Reader . If the input tape stops before the entire 
program has been read into the memory, then an examination of the 
characters punched on the tape just ahead of the place where it stopped 
will often provide an explanation. 

The Order Register and Control Counter . If the program 
is read in correctly but comes to a sudden unexpected stop, then the 
order register R^ will exhibit the order on which the program stopped 
while the control counter will contain a number one greater than the 
storage location from which the order pair came. With this information 
the progr amm er can often discover the cause of failure of his program. 

If the Illiac loops, observation of the slave tube will give 
an indication of how extensive the loop is. Then, if the machine is 
stopped and the contents of the order register and the control counter 
noted, we usually have enough data to identify the loop in our program. 

8.4 MODIFICATIONS TO PROGRAMS . Temporary modifications 
often need to be made to a program while it is being checked. These 
may be for the purpose of correcting the program, for temporarily 
arranging the punching of extra data, or for some other purpose. In 
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such cases it is undesirable to repunch the entire program tape. The 
modification can be carried out by "overwriting" the original program, 
that is by replacing some of the orders already placed in the memory 
by other orders. To do this a short extra tape is punched with suit¬ 
able directives and words so that, when it is read into the memory 
by the Decimal Order Input, the appropriate words of the original 
program are replaced by those on the tape. For example, 00 523K 
L5 4 f 4o 2F will cause the word at 523 to be replaced by L5 4 f 40 2F. 

In order to be able to use this technique of modifying a 
program we must: 

(a) Prevent the original program from being started 
before the modification is made. 

(b) Transfer control to the Decimal Order Input so 
that the modification tape can be read into the 
memory. 

(c) Start the program after the modification has been 
made. 

This can be done efficiently if the original program ends 
with a stop transfer of control to the Decimal Order Input followed 
by a transfer of control to the program. An example might be 24 999N 

26 93 N if it were necessary to transfer control to location 93 to 

start the program. 

If we wish to modify the program before starting it, all 
that is necessary is to place the modification tape in the reader 
before moving the switch to START. This will cause the modification 
tape to be read instead of the 2 6 93N. The modification tape will 
naturally have to end with 26 93N in order to start the program. 
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If the original program ends with a stop transfer of control 
to itself instead of to the Decimal Order Input a slightly more compli¬ 
cated modification tape is necessary. An example might be 2b 93N- 
Then the sexadecimal order pair 2 6 3^7 00 000 is punched at the 

head of the modification tape. When this tape is placed in the reader 
the machine is completely restarted by setting the order register to 
82 40F 40 F and the control counter to zero without, however, clearing 

the memory. This results in the first order pair on the tape (namely 
26 3F7 00 000) being transferred to location zero and being obeyed, 

so that control is transferred to 3 F 7 sexadecimal or 999 decimal (the 
Decimal Order Input), and the rest of the tape is read in the usual 
way. This is a bootstrap start. See Section 5«H« 

Corrections . As each coding blunder is found, a modification 
tape should be repunched to include all of the corrections. It is not 
worthwhile repunching the entire program tape until all or a large 
number of coding errors have been found. 

8.5 BLOC KING ORDERS . This is the name given to control 
transfer orders which are used to replace normal program orders, so 
that some printing or checking can be done at the point of replace¬ 
ment. Before control is restored to the program the replaced order 
of the program is executed and the contents of the arithmetic registers 
are restored. Thus the original program is unaffected, but the 
extra orders that are obeyed can be utilized to do printing of desired 
data. 
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An example may illustrate this: Let us suppose that 

(a) we have a program in which we wish to 

print N(19) after the left-hand order of 
the order pair 40 9F L5 29 F at storage 
location 100, has been obeyed, 

0>) we wish to preserve the quotient register 
but not the temporary storage of the print 
subroutine, 

(c) locations 800-804 are unused by the program, 

(d) the program print routine starts at 200, 

(e) the original program ends with 24 999N 26 nN. 

We then prepare the following modification tape: 


MODIFICATION TAPE COMMENTS 


00 

100K 

Directive 

4o 

9F 

Plants blocking order in progr 

26 

800F 


00 

800K 

Directive 

S5 

F 


4o 

4l 

Stores Q at 8 o4f 

L5 

19F 

Enters print subroutine 

50 

1L 


26 

2 OOF 


L5 

29F 

Does omitted order 

50 

4L 

Restores Q 

26 

101F 

Control back to program 

26 

nN 

Starts program 


8.6 TYPES OF CHECKING ROUTINES . There are several types 
of checking routines. One prints out the contents of certain 
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memory locations after a program has stopped. We call this a post 
mortem routine. 

A second type takes a given program and allows it to be 
obeyed order-by-order while printing out information about the course 
of the program. It is often called a sequence checking code. 

A third type arranges for information to be printed Out at 
specified points of a program. It is usually known as a check point 
or blocking order routine. 

8.7 POST MORTEM ROUTINES . Library Routines C3, C4, C5 . 

These are printing routines and have been arranged so that they can be 
used with little or no preparation on the part of the programmer. They 
are read into the memory with bootstrap input routines and are located 
in storage locations at the end of the memory, the longest occupying 
locations from 986 to 1023 and using locations 0, 1, and 2, as temporary 
storage. The end of each of the post mortem tapes contains 100 two- 
decimal-digit numbers. These numbers are used to specify the locations 
from which printing will occur. Suppose, for example, that we wish to 
know the order pairs in memory locations 540 to 549. Then we read in 
Routine C5. When it stops we find the number 54 on the end of the tape 
and place it in the reader. When the START switch is moved, the order 
pairs in locations 540 to 549 will be printed and then the Illiac will 
stop. If we start again we will get the order pairs in locations 550 
to 559. Codes C3 and C4 perform similar functions for decimal fractions 
and decimal integers. 
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The programmer should keep in mind the storage locations 
used by the post mortem routines so that he will be able to take 
full advantage of these important checking routines. 

8.8 POST MORTEM VERSION OF THE DECIMAL ORDER INPUT . Library 
Routine Cl . This is a very important checking routine and is usually 
the first one used after a program failure. It compares the contents 
of the memory with the contents of the input tape. Only discrepancies 
are printed out, enabling programmers to discover which orders of a 
program have become altered while the program was in the memory. This 
is an aid in making sure that orders which should have been modified 
have been modified correctly and that no order has become modified 
accidentally. 

This routine is used in the following manner: The Post 
Mortem Version of the Decimal Order Input is read into the memory 
in the usual way with a bootstrap input routine. It occupies storage 
locations 962 to 1023. The program tape is then placed in the reader 
after its copy of the Decimal Order Input, so that the Decimal Order 
Input is not read into the Illiac. Then when the machine is started 
the program tape will be read and the words created from it compared 
with the corresponding ones in the memory. When a discrepancy is 
found it will be printed on a line giving, first, the location at 
which the descrepancy was found, second, the word read from the tape 
printed as an order pair and, third, the word found in the memory, again 
printed as an order pair. Thus a typical line of printing might be 
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3^5 L5 000 4o 354 L5 546 4o 354 

indicating that the left-hand address of the order pair in location 345 
had become modified, taking the value 546. 

It should be pointed out that closed subroutines which have 
been used will usually have their links printed. This gives an indi¬ 
cation of the part of the program from which they were last called in. 
Interludes cause a large amount of printing because the contents of 
the interlude locations are changed twice in the course of input and 
are printed out both times. 

The Post Mortem Version of the Decimal Order Input causes 
the memory to be changed to the original state as it is being executed, 
so that if the program is started the original program will be performed 
again. 

If desired, only selected parts of the tape need be compared 
with the contents of thememory. However, when doing this it is 
necessary that all the preset parameters pertaining to that part 
of the program be input and that the selected part of the program 
begin with a directive. 

All post mortem codes occupy locations at the end of the 
memory and all of them use storage locations 0, 1, 2 as temporary 
storage, so that it is desirable when coding not to use these storage 
locations for constants or numbers which may need to be printed out. 

It is worthwhile to note that if addresses which are to be changed 
are initially read in with their final values then no printing will 
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take place on the post mortem if they have been modified correctly. 

8.9 THE ADDRESS SEARCH ROUTINE . Library Routine C2 . 

Programs sometimes fail because of a transfer of control to an order 
which causes the machine to stop. The usual order causing the stop 
is a zero left shift order because the memory is normally cleared to 
zeros before a new program is read in. 

Under such circumstances the offending transfer of control 
order may be searched for with the aid of Library Routine C2. The 
search routine is read into the memory, occupying storage locations 
normally occupied by the Decimal Order Input. Next the address to 
be searched for is read into the machine as a three character sexa¬ 
decimal address. The routine then searches the memory (exclusive 
of itself) for order pairs containing this address. When found they 
are printed out (in sexadecimal form) together with their location 
(in sexadecimal form). 

The routine naturally has other similar uses. For example, 
if it is known that some number becomes modified but it is not known 
why, then the store order which does the damage can be sought for in 
the above manner. 

8.10 SEQUENCE CHECKING CODES . Library Routines D2 and D3 . 
These routines control a program order by order and print out suitable 
information about the execution of each order. This enables the action 
of a program to be traced order by order. They use a blocking order 
technique which enables a selected part of the program to be checked 
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in this manner. These codes are very slow on account of the printing 
involved and should not be used blindly . 

Routine D3 prints out the function digits of the orders which 
are actually obeyed, starting a new line of printing whenever a control 
transfer order is obeyed. This enables the sequence of orders obeyed 
by a program to be traced. 

Routine D2 prints in full each order that is obeyed, having 
one order pair per line of printing. After each store order or address 
order, the number transferred to the memory is also printed. This 
enables every step of a program to be completely checked. However, 
the amount of printing is such that this routine is very slow. It 
should be used only in the final stages of tracking down an elusive 
fault. 

8.11 CONTROL TRANSFER CHECK . Library Routine Pk „ This 
routine takes charge of a program and allows it to be obeyed order 
by order. Bach transfer of control that is obeyed Is placed in a 
list kept in a specified place. The list is cyclic, that is to say, 
the later entires overwrite th« earlier ones In a cyclic fashion. At 
the end of a program, the list can be printed, so that it can be dis¬ 
covered how the program reached its final end. There is no printing 
during the execution of the program, so that this routine allows the 
program to be obeyed at speeds much greater than those of the routines 
which print. 
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8.12 THE CHECK POINT ROUTINE . Library Routine D l. This 
routine is designed to print out intermediate information about some 
other program. It uses the blocking order principle, and the programmer 
prepares a specification tape to describe the kinds of information he 
wishes to obtain. It is possible to go through iterative loops and 
print results on various passages through the loops. Data can be ob¬ 
tained as an order pair, a right-hand address, a left-hand address, 
a 10 character sexadecimal word, a signed integer, a signed 12 decimal 

place fraction or a signed 5 decimal place fraction. 

Library Routine Dl is a very powerful checking routine be¬ 
cause of the great latitude given-the programmer in choosing where 
and how he will obtain information and because it utilizes the pro¬ 
grammer's own knowledge of his code. 

8.13 CONCLUSION . The fact that it will be necessary to 

check a program should be kept In mind when writing It. Its storage 
locations. Including those of the temporary storage, should not Inter¬ 
fere with routines likely to be used In diagnosing its faults. The 
program tape should end with a stop control transfer to the decimal 

order input as explained in Section 8.4. 

Occasionally a blunder turns up in the blind spot of the 
programmer and it appears to be Impossible to find it. It is not much 
comfort to him to point out that the blunder can be found by continually 
narrowing down the section of the code in which it is known to be. This 
is tedious but eventually all blunders yield to this technique. A final 
„ord to this chapter might be that all programmers make mistakes but good 

programmers find theirs first. 
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CHAPTER 9 


TAPE PREPARATION 

g.l THE ILLIAC INPUT . The input unit of the Illiac 
is a photoelectric tape reader that transfers binary digits from a 
punched paper tape to the A register. The tape preparation equip¬ 
ment is used to translate instructions and data from the programmer's 
manuscript into a binary-coded punched tape acceptable to the Illiac. 

9.2 THE ILLIAC OUTPUT . Output from the Illiac is 
usually in the form of punched tape. This tape may be printed on 
any of the Teletype page printers that are a part of the tape .pre¬ 
paration equipment. The page printer performs a conversion from 
the binary-coded representation on the tape to the sexadecimal or 
decimal characters on the printed page. 

9.3 THE PERFORATED TAPE . SEXADECIMAL TAPE CODE CHARACTERS . 
The tape preparation .equipment presently in use consists of Teletype 
equipment which has been modified to correspond to the Illiac tape 
code. This tape code is shown in Figure 9.1. The paper tape can be 
punched in any one of 6 positions across its width. One of these 
positions is always punched. This is the feel hole or sprocket hole 
position. Feed holes are smaller than the other holes. 

Four of the remaining positions are used to represent the 
2* sexadecimal characters. There are 1 6 keys on the Teletype tape 
punch labeled 0 through 9 and K, 3, N, J, F, L. When one of these 
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oi23^56789KSNJFL 



Figure 9*1 

Sexadecimal Tape Codes 


Delay Space 



I 


Tabulate Carriage Return 

and Line Feed 


Figure 9-2 
Format Tape Codes 
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keys is pressed a corresponding pattern of holes is punched across 
the width of the tape. Since a hole or absence of a hole is a binary 
affair, we speak of the character punched in the 4 positions across 
the tape as a binary tape code for the corresponding sexadecimal digit. 
If we speak of only the 10 sexadecimal characters 0-9 and their 
corresponding binary codes, we call the characters punched in a tape 
binary-coded decimals . 

The Illiac tape code for sexadecimal characters is very 
easy to memorize since the hole positions across the tapte simply 
correspond to powers of 2 from 2 to 2 . 

Page Printer Format Characters . Punched tapes are con¬ 
verted to a printed page by means of a page printer. Four addi¬ 
tional tape codes shown in figure 9.2 are used for format control, 
causing the printer to space, return the carriage and advance to 
the next line, tabulate, and delay, i.e., perform one printing 
cycle without printing anything. The functions are similar to 
those of a typewriter, but the tabulate feature occurs on only 
one printer. Successive tabulate characters set the carriage to 
10, 25, 40 and 55 spaces from the left margin. The carriage 
return-line feed and tabulate characters must always be followed 
by a delay character to allow sufficient time for the carriage to 
move. 
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Notice that these codes all have a hole in a position that 


is never punched for a sexadecimal. This is called the 5th hole 
position . When a punched tape is placed on the Illiac tape reader 
the circuits cause any character with a 5 th hole punched to be skipped 
over, i.e., not read into the computer, when reading in the normal 
way. (However, a special input order will read such a character; 
see the order code.) This means that in the preparation of an 
instruction tape, these format codes can be interspersed with sexa¬ 
decimal order digits in any desired way. Then when an instruction 
tape is printed these tape codes will control the printer. The usual 
method in preparing instruction tapes is to follow each order pair with 
a carriage-return and line feed code. This produces a single column 
print of words. 

Tape Codes for Letters . The printer will also print all 
the letters of the alphabet for identifying headings and the like. 

Tape codes are provided for these and are shown in the complete tape 
code list. Figure 9-3- To accomodate all of these symbols, there 
are also numbers shift and letters shift tape codes, analogous to the 
shift-lock on an ordinary typewriter. 

A list of instructions to be used with the Illiac and the 
Decimal Order Input (See Chapter 5) to produce any of the tape codes 
shown in Figure 9.3 is given in Table 9.1. 








PRIMED CHARACTERS 


MACHINE ORDER 
(DECIMAL ADDRESS) 


Space 92 963F 
Carriage Return and Line Peed 92 131F 
Delay 92 515F 
Letters Shift 92 259F 
Numbers Shift ' 92 707F 


AFTER 

AFTER 


LETTERS SHIFT 

NUMBERS SHIFT 


A 

) 

92 387F 

B 

( 

92 195F 

C 


92 835F 

D 

$ or tab* 

92 67F 

E 

3 

92 194F 

F 

F 

92 898F 

G 

= 

92 579F 

H 

apos trophe 

92 771 F 

I 

8 

92 5UF 

J 

J 

92 834F 

K 

+ 

92 642F 

L 

L 

92 962F 

M 


92 643F 

N 

N 

92 770F 

0 

9 

92 578F 

P 

0 

92 2F 

Q 

1 

92 66 F 

R 

4 

92 258F 

S 

- 

92 706 F 

T 

5 

92 322F 

U 

7 

92 450F 

V 

> 

92 323F 

w 

2 

92 130F 

X 

/ 

92 451F 

Y 

6 

92 386F 

z 

X 

92 899F 


Table 9-1 

Coded Tape Operation of the Teletype Printer 


*Only one printer has a tabulation mechanism. 
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9.4 CLASSIFICATION OF OPERATIONS . The operations in 
preparing tapes may be listed as follows: 

1. Punching a tape from a manuscript by 
means of a keyboard tape punch. 

2. Punching a tape from one or more previously 
punched tapes. This operation is called 
reperforating and is used for duplicating 
tapes, joining short tapes, and making 
corrections. 

3. Printing a tape on a page printer. 

4. Comparing two tapes to see if they are 
identical. 

5 . Comparing the manuscript with a printed 
copy (proof-reading). 


It is necessary to punch the main part of the code from 
a manuscript by means of a keyboard. It Is not necessary to punch 
routines that are in the routine library. 

When all the parts of the code that have been punched from 
manuscript are correct, obtain any necessary routine tapes from the 
routine library. Using the reperforator, these tapes can be assembled 
into one long problem tape. 

Every time a tape is reperforated it is wise to check the 

copy by comparing it with the original tape. 

9.5 DESCRIPTION OF EQUIPMENT . Several Teletype units 
are provided for making instruction tapes. These include a keyboard 
perforator unit, a page printer unit, a reperforator unit, a combination 
unit, and a tape comparer unit. 
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F i g u r-e 9.4 

KEYBOARD PERFORATOR UNIT 
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Keyboard Perforator Unit . This unit is used to punch tape 
from a manuscript; see Figure 9.4. Turn on the power switch under the 
table by pulling it forward. The power supply takes about a minute to 
warm up. 

The sixteen sexadecimal keys that will be used most are 
in the center of the keyboard and arranged to be operated with one 
hand. It is recommended that a touch system be learned both for maxi¬ 
mum accuracy and maximum speed. 

If a wrong key is pressed, operate the backspace lever and 

then the space bar. This will punch the space symbol (all five holes) 
over the incorrect symbol and so be skipped on the Illiac reader if 

the usual 80 instructions are used. 

Page Printer Unit . This unit is used to make a printed 
page corresponding to a punched tape; see Figure 9*5* Turn on the 
power switch. The power supply takes about a minute to warm up. 

In order to operate the printer mechanism the page printer 

motor switch should he turned on* 

The function of the transmitter-distributor is to read paper 
tape and translate the binary tape codes into electrical signals. 

These signals are sent to the page printer. The binary tape code 
is sensed by 5 pins which are periodically pressed against the tape. 

A hinged lid holds the tape down against the pins. The tape is moved 
by a sprocket wheel which engages the feed holes punched in the tape. 
Place a tape in the reader and close the lid. Then turn the start- 
stop switch to start. The tape will now be transmitted to the page 

printer. 
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POWER SWITCH 



TIGHT TAPE 
STOP LEVER 


Figure 9.5 
PAGE PRINTER UNIT 
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There is an end of' tape stop device that stops the transmitter- 
distributor when the end of the tape passes under the lid. A tight tape 
stop bar will also stop the unit and keep the tape from being torn should 
the tape become snarled or tight and raise the bar. 

Reperforator Unit . This unit is used to duplicate tapes; 

see Figure 9.6. 

Turn on the power switch by moving it towards the back of the 
table. After about a minute, the power supply will warm up and the 
reperforator motor switch should be turned on. 

The transmitter-distributor has been described in the pre¬ 
ceding paragraph. Place a tape in the reader and turn the start-stop 
switch to start. A copy of this tape will be made by the reperforator. 

When a few corrections need to be inserted in a tape, punch 
a second tape containing only the corrections. Then use the reperforator 
unit to duplicate the correct portion; stop; insert the second tape in 
the reader and reperforate the correction; then replace the first tape 
and continue with the correct portion. 

Combination Teletype Unit . This unit, shown in Figure 9-7, 
combines the functions of those units already described. 

Underneath the top of the table are two switches. Switch 
10 turns on the dc supply which must be on before any of the other 
units will operate. On some of the tables it takes about a minute 
for the dc supply to warm up. 

Switch 9 turns on the motor of the transmitter-distributor. 
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Figure 9.6 
REPERFORATOR UNIT 
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ON 



© DOOR TO P0NCHIN6S drawer COMBINATION TELETYPE UNIT 

(g) REPERFORATOR MOTOR SWITCH 






The transmitter-distributor has already been described. It sends 
signals to both the page printer and the reperforator and will operate 
either or both of these units if the corresponding motor switch is on. 

Turn on the reperforator motor by means of switch 13. The 
reperforator will now duplicate any tape that is run through the 
transmitter-distributor. 

Pressing lever 11 will cause the reperforator to punch 
space tape codes. Pressing button 8 will cause the reperforator to 
punch zeros. However, one unpredictable character will be punched 
at the end of the zeros. Zeros are useful at the beginning of tapes 
to provide space for writing identifying information. They should be 
followed by spaces to make it easy to set the tape in the Illiac 
reader. 

To operate the page printer turn on motor switch 7 B. The 
page printer will now print any tape that is run through the trans¬ 
mitter-distributor. The carriage may be returned manually by pressing 
lever 2. Crank 1 will turn the platen. 

It is possible to operate the page printer and the reperforator 
at the same time. However, it is possible for the reperforator mec han ism 
to make an error that would not appear on the printed copy. If the 
printed copy that is to be proof-read against the manuscript is pre¬ 
pared from the final tape, it is quite certain to show any errors in 
the tape as well as any possible errors in the printing mechanism. 



Set switch 7A to "tape". Tape can now be perforated from the 

keyboard. 

It is also possible to type on the keyboard, and perforate 
type in the reperforator at the same time. To do this turn 7A to 
keyboard, 7B on, and 15 on. There is no way to back-space tape that 
comes from the reperforator. It is also more difficult to read the 
last character that has been punched on the tape than when the keyboard 
perforator is used for punching tape. The transmitter-distributor 
may be used to reperforate correct portions of a tape. The tape should 
be stopped before the error by means of switch 5* TyP e the corrections 
from the keyboard. Then move the tape to the next correct section 
and continue reperforating. 

It is also possible to print with the page printer at the 
same time a tape is being punched from the keyboard. To do this turn 
on the printer motor by means of switch 7B and set switch 7A to 
"keyboard and tape". A disadvantage of this method is that the keys 
cannot now be operated as fast as when punching tape alone. There is 
no way to change errors on the printed copy that are corrected on the 
tape with the backspace lever. Since the printed copy is made by a 
separate mechanism at the same time that the tape is punched, there 
is also a small possibility that it may not correspond to the holes 
in the tape. This copy should not be used for proof-reading. 

Switch 7B can be turned on and off without affecting anything 
being reperforated if switch 7A is placed in the "tape" position. 
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FIGURE 9.8 
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before 7B is operated, and kept there while 7B is off. 

Comparing Tapes. A tape comparer consists of two tape readers 
connected to checking circuitry. Two kinds of tape comparers are avail¬ 
able in the Laboratory, one kind using mechanical readers and the other 
using photoelectric readers. When tapes are compared in the mechanical 
readers, it is necessary only to place the two tapes in the readers and 
to turn the power on. The tapes will advance in synchronism as long as 
the binary codes in the two tapes agree and will halt if a disagreement 
is detected. On of the machanical units will halt one row of holes 
beyond the disagreement and the other will halt on the holes which dis¬ 
agree. When the tapes stop, the power should be turned off and the tapes 
inspected visually. They may then be marked if in error, replaced in 
the readers and the comparison continued. One of the two mechanical 
units has push-buttons for advancing either tape one step at a time 
while they are stopped on a disagreement. 

The tape comparer using photoelectric readers is more versatile 
than those with the mechanical readers. By means of a switch it is 
possible to choose three modes of operation in which the reader will 
"skip no characters", "skip all 5th hole characters" or "skip all spaces. 
A set of neons is associated with each tape reader and when the tapes 
stop on an error the neons indicate the character read from each tape. 

As in one of the mechanical readers, it is possible to advance either 
tape one step at a time by means of switches. A full description of 
the operation of the photoelectric tape comparer is provided in Internal 
Report 62, and a view of the control panel is shown In Figure 9-8. 
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Short tapes may be checked more rapidly visually by super¬ 
imposing them and holding them up to a strong light. 

9-6 GENERAL REMIRKS ■ The misreading of even a single 
hole in a long tape can cause enormous changes in the behavior of 
Illiac. It is important to keep tapes clean. Dirt from tapes 
accumulates in the transmitter-distributor tape reading pins and 
causes errors. Worn or torn tapes will also cause errors. 

Do not mark tapes with a waxed pencil: use a lead pencil 

or ink. 

Library tapes and other tapes which need to be handled 
frequently are on a grey parchment stock which is heavier than normal 
tape stock. One of the reperforators is mechanically adjusted to 
punch this heavy stock. 

It is generally more convenient to make corrections to be 
made on the tape itself, as well as on the printed copy. 

If corrections are needed at the end of a long tape, it is 
sometimes more convenient to reperforate the tape tail-end first, 
so that the corrections can be made at the beginning, and then the 
reperforator left running unattended. 

A hand punch is available and can often be used to change 
one or two characters and so avoid reperforating a long tape. Holes 
punched with it should be inspected carefully to see that they line 
up with the other holes in the tape. If in doubt, it is better to 
have the extra hole slightly oversized so that it is certain to be 
read by the Illiac reader. 
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9.7 HIGH SPEED PRINTER . This printer performs the 
same function as the slower printers except that it operates at 
10 characters per second. The turn-on switch is on the lower 
front side. Carriage return and line feed buttons are provided 
on the front. After printing is complete the paper should be 
advanced before shutting the printer off. Otherwise unwanted 
characters will be printed on the page as the unit stops. 

9.8 HIGH SPEED REPERFORATOR . The high speed reper¬ 
forator uses two photoelectric readers to provide a flexible 
means of preparing tapes. Four switches and a push button con¬ 
trol the operation as follows. 

Switch S4, Stop-Run . This switch has only two positions 
which determine whether the unit is in a standby position or 
operating in accordance with settings of the other switches. 

It should be at STOP to change any switch except S2. 

Switch SI , Reader Selection . This switch has three posi¬ 
tions, allowing a tape to be read from the front reader or the 
back reader or allowing carriage return-line feed characters to 
be punched on the new tape while neither reader operates. 

Switch S3, Rode Control . This switch, in conjunction with 
32, determinesthe mode of reading a tape and punching another. 

It may be set as follows: 

a. "Single Character Punching" to punch one character and stop. 

b. "Stop Punching on CR & LF" to punch up to and including the 
next CR & LF character and then stop. 

c. "Single Character Bypassing" to skip one character and stop. 

d. "Stop Bypassing on CR & LF" to skip all characters up to and 
including the next CR & LF and then stop. 
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Switch S2, Start , Obey , Ignore , This switch obeys the 
instruction of S3 when set to OBEY. If moved to START after a 
stop it will cause the instruction of S3 to be executed once. 

If set to IGNORE it will cause all of the stops specified by 
S3 to be ignored. 

Set Button . This is a reset control which should be 
pushed when the unit is first started and whenever* a tape is 
manually shifted in either reader. 

Neon Indicators . A row of neons on the panel indicates 
the last character read and either punched or bypassed. 
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CHAPTER 10 


CALCULATION OF RUNNING TIME 


Each programmer should estimate as well as he can the 
length of time his program will run. This information is needed 
for efficient scheduling of machine time and to enable the computer 
operator to decide whether a program may have failed to stop when 
it should have. 

The estimation of running time consists of summing up the 
times for executing all of the orders of the program, taking into 
account that many orders are executed more than once. Thus, the 
time is calculated for one passage through an iterated loop and then 
multiplied by the number of trips through the loop. The number of 
passages through a loop is not always known (as in the square root 
code where the number of passages depends upon the quantity whose 
square root is being calculated) and the programmer may have to 
make rough estimates for such cases. He can make use of the code 
checking periods to help him here. 

10.1 ORDER TIMER. The programmer must know how much 
time the Illiac requires for each order. The order times are given 
in Table 10.1, the time required for getting the order from the 
memory being included. 

The values given in Table 10.1 are in some cases not 
exact. The 70# microsecond value for multiplication is a maximum, 
the t im e depending upon the multiplier. It could be as low as 625 
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microseconds . 


ORDER TYPE TIME 


On, In 

16 

n microseconds 

s ^i 5j J 

55 

microseconds 

3 not executed 

18 

microseconds 

6 

800 

microseconds 

7 

700 

microseconds 

8o, 81 

4 

milliseconds/character 

82, 92 

17 

milliseconds/character on punch 


1 millisecond/character on display 

K, S, P, L 

90 

microseconds 


Table 10.1 
Order Times 

The time required to execute an 80 or 8l order so as to 
read one character from the tape is 4 milliseconds. Since most 
programs use conversion routines for input it will be found that 
the time is closer to 4 l/2 milliseconds for decimal input. 

10 * 2 EXAMPLE OF RUNNING TIME CALCULATION. Let us con¬ 
sider a simple example. The program given in Table 10.2 will trans¬ 
fer the contents of memory locations 100 to 199 into memory locations 
200 to 299- How long will it take? 

The program in Table 10.2 consists of 7 words. There is 
a loop consisting of the 9 orders beginning with the left-hand order 
at location 0 and ending with the left-hand order at location 4. 
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This loop is executed 100 times. There are no other orders executed 
except the stop order. Thus we have the following orders: 

00 50K 
0 L5 (100)F 
40 (200)F 

1 L5 L 
L0 5L 

2 32 4L 
L5 L 

3 l4 6 l 
40 L 

4 26 L 
OF F 

5 L5 199F 
40 299F 

6 00 IF 
00 IF 
24 50N 

Table 10.2 
Repetitive Program 

Five L orders at 75 microseconds. 

Two 4 orders at 55 microseconds. 

One 2 order at 55 microseconds. 

One _2 order at 18 microseconds. 

The time T is then given as 
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T - (5 x 75 + 2 x 55 + 55 + 18 ) 100 microseconds > 

= 53,800 microseconds, 

= .054 seconds. 

10.3 A SIMPLE RUNNING TIME FORMULA . A rule which is 
accurate enough for most calculations is the following one: 

Let Nq = number of orders obeyed, 

= number of multiplication orders obeyed, 

Nj - number of division orders obeyed, 

N = number of shifts. 
s 

Then the running time T, exclusive of input and output orders, is 

N N + N 
m _ 0 s m d 

1 “ Y£r + "ScT + ~~2 milliseconds. 

In the example given above KL is 900 while N = N ■= w ~ 

0 s m iV d “ 

We thus have T = 9 OO /16 milliseconds or about ,056 seconds. 
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CHAPTER 11 


PREPARATION OF A COMPLETE PROGRAM 

The problems solved on the Illiac vary in complexity from 
programs requiring only the punching of parameters and data for use 
with a library routine to programs involving hundreds of words re¬ 
quiring the entire memory capacity of the machine. There is, of 
course, no "typical" program; however, the problem described in the 
following sections illustrates features common to many problems. 

The problem described is one that arose in the electron 
tube research laboratory of the University of Illinois.* The problem 
is first described qualitatively, and a general mathematical formu¬ 
lation is given. Approximations are then made to simplify the 
equations. The resultant equations are reformulated so that a 
library routine can be used for the solution, and reformulated again 
in terms of scaled variables so that machine range will not be ex¬ 
ceeded. The general organization of the program is then described, 
and the details of the coding are given. Those interested primarily 
in the programming may omit the introductory description. 

11.1 GENERAL STATEMENT OF THE PROBLEM . The purpose 
of the problem is to solve for the trajectories of electrons in an 

* The staff of the computer laboratory is indebted to Mr. Irving 
Kaufman of the electron tube research laboratory for permission 
to describe his problem and for his assistance in the preparation 
of this chapter. 
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electron cyclotron. An electron cyclotron is a device that accelerates 
electrons to energies of several million electron volts. The accelera¬ 
ting mechanism is a radio-frequency electric field in a microwave 
cavity. According to the well-known Lorentz force equation, the force 
on a charged particle is given by: 

f = q! + qv x B , (11.1) 

Here 

f = vector force (newtons), 

q = charge (coulombs), 

E = vector electric field (volts per meter), 

v = vector velocity (meters per second), 

B = vector magnetic flux density (webers per -square meter). 

It is seen that the component of force produced by the 
electric field E is parallel to E $ that produced by the magnetic 
fisld B is perpendicular to both v and B. The former causes accelera¬ 
tion or increase in energy? the latter merely changes the direction 
of motion. This combination gives rise to the mechanism of the 
electron cyclotron shown in figure 11.1. 

Here a microwave cavity, externally supplied with electri¬ 
cal energy, is immersed in a uniform steady magnetic field (shown 
by the dots) whose direction is perpendicular to the paper. The 
cavity has been constructed so that near the axis it has a time- 
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varying electric field E that is almost; errtirely parallel "to "the 

axis X-X' and whose magnitude is given by E = E ffl sineiOt. The raagni- 

7 

tude of E is on the order of 5 x 10 volts per meter, so that it 
m 

is possible for an electron passing through the 1 cm. gap of the 

5 

microwave cavity to acquire an energy of 5 x 10 electron volts. 

(This corresponds roughly to the relativistic rest energy of the 
electron). 


X 



Electron Cyclotron 
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The^action of this cyclotron is as follows; 

1* ver y high value of E forces some electrons out of 

the cavity wall during some parts of the cycle. (Let us consider 
here the trajectory of electrons emitted from surface R-R). 

2. The E-field accelerates the electron toward surface 

S-S. 


3* A fraction of the electrons emitted pass.through the 
orifice in S-S. The magnetic field as well as space charge pro¬ 
duces a slight bending of the trajectory. 

^• Ths magnetic field external to the cavity causes 
the electron to move in circular orbits . The orbit radius is pro¬ 
portional to the product of electron speed and total (relativistic ) 
energy. Consequently, a fraction of the electrons ejected are 

returned to the orifice in cavity surface R-Rj the rest collide with 
the walls. 


5• Those particles that re-enter the cavity when the field 
is in its accelerating phase (recall, E - E ffi sincOt) acquire 
additional kinetic energy during transity through the gap. They conse 
quently execute larger orbits externally. 

6. For those electrons that maintain the favorable phase 
conditions of (4) and (5), above, an indefinite number of orbits 
is possible. For a 1 cm. cavity gap, each orbit corresponds to a 
step in energy of approximately 1/2 million electron volts. 
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11.2 TRAJECT ORY EQ U ATIONS . The combination of the 
relativistic equation of motion for a particle and the Lorentz force 
equation yields, for electron motion in the cavity, the equation 

1/2 = e [ I (E m s in to t - Bv ) + J (v^B) ] (ll 

Here: 

m^ = rest mass of electron, 

V x = dX//dt = com P one nt of velocity in X-direction, 
v y = dy/d.t = component of velocity in Y-direction, 
c = velocity of light, 
i = 'unit vector in X-direction, 
j = unit vector in Y-direction, 
t = time, 

e = magnitude of electronic charge, 

E m = peak electric field intensity, 
cu = angular frequency of cavity field, 

B = magnetic flux density. 

For motion exterior to the cavity, the term E sinwt 

m 

should be deleted from equation 11.2. The two resulting vector 
equations are then to be solved for electron trajectories. 

To reduce the complexity for a first order solution, 
three approximations are made: 
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1. The effect of the magnetic field inside the cavity 
is neglected. 

2. The effect of v in the gap is neglected. 

3- Particles exterior to the cavity are taken to move 
in circular arcs. The arc length is taken as the difference be¬ 
tween a full circle and the cavity gap length, (in other words, 
the chord length of the gap is set equal to its arc length.) As 
before, the radius of the circle is proportional to the product 
of velocity and total energy. These approximations leave only one 

differential equation to be solved; for now, after some manipulations, 
in the gap 

h [ m^ v (l - ] / dt = e sintut , (q 

while outside of the gap 

9 i,n + l - ^ t 1 - K, (1 - p 2 ) V2 / (5 ] (1 . f) -1/2 


+ © e ^ n - (2n) (n + 1) 


(11.4) 


The symbols used in 11*3 and 11.4 and not previously 
mentioned are: 

v = electron speed (meters per second), 

P = v/c, 

n = number of the orbit in which moved; i.e., for 

first cavity crossing and exterior orbit n = 1, etc. 
Q e, n = the P hase angle (radians) in the microwave RF cycle 
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at which the electron leaves the gap in orbit n. 


(See Figure 11.2). 

= the phase angle (radians) in the RF cycle at which 
the electron re-enters the gap after completing 
orbit n. (See Figure 11.2). 

K x = 6.51599 

k 2 = 0.092910 




11 •5 PROCESSING OF EQUATI O NS FOR SOLUTION BY THE ILLI AC. 

Ele c tr on Traje ctory in the Gap . The equation of motion of 
the elctron in the gap is given by equation 11. 3 , which is reduced 
to a set of first order ordinary differential equations so that the 
solution can be obtained using library routine F 1. 

Using as variables 0 = urt, u = v/u), we can reduce equation 
H-3 to the pair of first order equations: 
dx/d© = u, 

du/d© = (eE m / m 0 o) 2 ) [ 1, - (a? u 2 ) / c 2 ] ^ 2 s in 0 . 

The sine function may be evaluated by using two additional 
equations with appropriate initial, conditions: 

dfsin 0) / d0 = cos 0 or df/d© = g 
d(cos 0) / d0 - -sin 0 or dg/d© = -f 

The equation of morion of the electron in the gap has thus been reduced 
to the following set of four first order differential equations: 
dx/d© = u 

du/d© = (eE m / m 0 u£) [1 .. (J u 2 / c 2 ) ] 3/2 f 

df/d© = g 
dg/d© = -f 

In itial Conditions for Gap Trajectory Equations . An 
electron is introduced into the gap at x = 0 with negligible velocity 
v =u*i at a certain phase © Q =u)t. The initial conditions are then. 
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g = cos 9 ^ • 


for arbitrary 6^, 

x = 0, u = 0, f - sin © Q , 

The effect of variation of the parameter 9^ will he studied. 

Motion of the Electron E xt erior to the Cav ity. The elec¬ 
tron is accelerated while in the gap and leaves the gap at x = 1 cm. 
with a velocity ulu and phase © . It then follows a circular 

C jil f *■* 

orbit at constant velocity and re-enters the gap at x = 0. The 
phase at the time of re-entry can be calculated by equation 11.4, 

in which ft is equal totOu^ n / c. 

Re-entry Conditions for Gap Trajectory Equations . The 

re-entry conditions for the gap trajectory equations are: 

x = °, u = u e ^, f = sin 9 i#n+1 > S = cos Q i, n +1' 9 = Q i,n+1' 

where © , is calculated from 9 as indicated above. 

i,n+l 

Information to be Printed . For a given value of the 
initial phase 9 Q , the following information is desired for each of 
ten orbits: 

(a) The orbit index n, 

(b) The phase 9 of the electron leaving the gap, 

(c) The phase 9 ± n+1 of the electron re-entering the gap, 
(d ) The relative velocity j3 of the electron in the orbit, 

(e) The quantity [ 1 - p J > 

(f) The energy of the electron in the orbit. 
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The calculation is to be repeated with different values of the para¬ 
meter 0 ^. The parameter is printed for identification of results. 

11^ SCALING OF THE QUANTITIES FOR ILLIAC SOLUTION . 

Magnitude of Quantities Involved . For the rnage of in¬ 
tegration (electron trajectory in the gap), the phase angle 1 was 
known to be in the range 1 < G < 2.5 from a previous desk calculator 

solution. The design constants and physical constants are E , u) 

nr 

e, m Q , and c, which are, respectively, the peak electric field 

intensity in the gap, the angular frequency of the field, the charge 

of the electron, the rest mass of the electron, and the velocity of 

light. The gap length of 1 cm. is the maximum value of the variable 

x. The maximum value of u is easily caluclated, sincea>u/c I s 

the velocity of the electron relative to the velocity of light and 

is less than one. The variables f, g, and their derivatives are 

the well-known and well-behaved sine and cosine functions. The 

maximum value of the derivative du/dO is simply the constant 

eE /mu). 
m J 

The ranges of the unsealed variables and values of 
the constants are indicated in Tables 11.1 and 11.2. The scaling 
of the variables is indicated in Table 11. 3 . 
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SYMBOL 

0 

x 

u 

f 

g 

dx/d9 

du/d© 

df/d© 

dg/d© 


DESCRIPTION 

NUMERICAL 

Min. 

VALUES 

Max. 

SOURCE 

Phase 

1 radian 

2 .5 radians 

Previous desk 

calculator 

solution 

Distance 

0 

0.01 meter 

Design constant 


0 

0.017 

me ter/radian 

aJu < c 

f = cos 9 

-1 

1 


g = sin 9 

-1 

1 


dx/d© = u 

0 

0.017 m/rad. 



-0.029704 

0.029794 

sin 9 

(1 - |3 2 ) 5//2 < 

df/d© = - sin © 

-1 

1 


dg/d© = cos © 

-1 

1 



Table 11.1 

Table of Variables and Derivatives 
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SYMBOL 


DESCRIPTION 


NUMERICAL VALUE 


SOURCE 


e / m Q 

Ratio of charge to 
rest mass of an 
electron 

1.7592 x 10 11 
coulombs/kg. 

Physical 
constant 

G 

Velocity of light 

2.99776 x 10® 
m/sec. 

Pnysical 

constant 

E 

m 

Peak electric 
field intensity 

5.2266 X 10 7 
volts/meter 

Design 

constant 

u> 

Angular frequency of 
the field 

( 2 jt)( 2 . 8 )(l 0 9 ) 
rad/sec. 

Design 

constant 

X 

max 

Gap length 

0.01 meter 

Design 

constant 

e E„ 
id 

'2 

m 0 (b 


0.029704 

Derived 

constant 

K 1 


6 •31599 


K 2 


0.092910 



Table 11.2 

Table of Constants (mks units) 
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SYMBOL 

NUMERICAL 

Min. 

VALUE 

Max • 

RELATION TO UNSCALED VARIABLE 

0 

0.1 

0.25 

0 = 0.1 0 

y 

0 

0.301 

y = 0.1 x 

p 

0 

0.017 

p = u 

r 

0 

0.005 

r = 0.005 f 

s 

0 

0.005 

s = 0.005 g 

dy/d0 

0 

0.017 

dy/dfl = dx/d9 

dp/d0 

- 0.29704 

0.29704 

dp/d0 = 10 du/d9 

ds/d0 

- 0.05 

0.05 

ds/d0 = 0.05 df/d5 

dr/d0 

- 0.05 

0.05 

dr/d0 = 0.05 dg/d9 


Table of Scaled 

Table 11 .5 

Variables and 

Derivatives 

The 

Scaled Equations 

. The scaled equations become 


dy/d^ = p , 

dp/d0 = E [ l/2 - u 2 /G ] 5//2 (2s/0.0l) , 

dr/d0 = -10s , 
ds/d0 = lOr . 

The quantity in brackets in the equation for dp/d0 has 
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been scaled by l/2. The constants G and E are: 

G = 2c 2 /(J = 0.00058078 j 

p 

E = 20 f2 eE /mu) = 0.84016 . 
m u 

The scaling for dp/d/ and for the bracketed expression 
have been absorbed by the constant E. 

Scaling for External Orbit Recycling Equation . The time 
phase of an electron re-entering the gap is calculated as follows: 

°i,n+l = Kj_ [ 1 - % (1 -o) 2 u 2 /c 2 ) l/2 / (urn/c) ] [ 1 - (6>u/c f ] " l/2 

+ " ( n + l) (2rt) , 

where = 6.31599 and = 0.092910. Since / = 9/l0 and n can be 
as great as 10, 


^i,n+l = ( K i/ 10) ^.l) [ 1 - K 2 (1 - u)p 2 /c 2 ) l/2 / (I0u>p/c) ] 


[ 1 - (wp/c r i 


-1/2 


- 2n(n + 1) / IOC + / 


e,n 


where / and p are values at the time of exit after the nth tra- 

G y n 

jectory of the gap. After performing the indicated subtraction, the 
quantity within the brackets J is on the order of 0.01. 

Additional Calculation for Printed Results . The additional 
calculations below are to be performed. 


p =<JL>u/c 


(relative velocity) 
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2 

u n /u 0 = l/l - (ratio of total energy to rest energy) 

T = m 0 c 2 [ (1 - (? 2 ) - 1 / 2 - 1 ] = (511.24) (lO 3 ) [U n /U 0 - 1 ] 
(electron kinetic energy) 


After scaling we have 


p =<*jp/c 


0.01 u n /u 0 - 


s-8„ 


o.oi vfi" 


[ 1/2 - p 2 /G ] 


1/2 


10 T = 0.51124 [ 0.01 U n /U Q - 0.01 ] . 


11.5 CODING OF THE PROBLEM . 

Organization of the Program . The following steps are to be 

done: 

(1) Read parameter / Q from tape and print. 

(2) Set initial conditiofts for integration in gap. 

(3) Perform one stop of integration. 

(4) Test to see if last step has been performed; if not, 
repeat step 3- 

(5) Perform external orbit recycling calculation. 

(6) Calculate and print data desired for each orbit. 

(7) Test for last orbit; if not, return to step 2. 

(8) Stop, then return to stop 1. 

Use of Library Routines . The following library routines 

were used: 
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(1) Decimal Order Input X 1 

(2) Constant Listing Auxiliary X 3 

(3) Differential Equations F 1 

(*0 Sine-Cosine q q 

(5) Decimal Number Input N 3 

(6) Print p q 

(7) Square Root p 1 


Use of Parameters . In order that the coding could be 
independent of the allocation of memory space, parameters in 
addition to those required by the differential equations F 1 were 
assigned for memory locations of the subroutines. The assignment 
of parameters is as follows: 

location of first variable, 

location of scaled derivative of the first variables, 

location of first word of temporary storage for 
differential equations routine, 

number of differential equations to be solved, 

location of first word of auxiliary subroutine, 

temporary storage, 

square root subroutine R 1, 

initial conditions setting program, 

main program, 

print routine PI, 

decimal number input routine N 3, 


S3 

SU 

55 

56 

57 

58 

59 
SK 
SS 
SN 
SJ 
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SF sine-cosine routine T 1, 

SL differential equations routine F 1. 

The parameters S3-S7 are those required for the differential 
equations routine. 

Use of Differential Equations Routine. In using the 
differential equations routine, a number of choices must be made, 
as indicated in the specification sheets. These include: 

(a) selection of increment length h and the parameter m, 

(b ) method for handling the independent variable 0. 

It has been noted that two first order equations are required for 
solution of the second order differential equation of motion of the 
electron in the gap. Two additional equations were included to 
evaluate the sine function required. The range of integration is 
from y = 0 to y = 0.001 where y is one of the dependent variables. 
Since it is not possible to predict the relationship between y and 0, 
it is not possible to determine an increment h of 0 in such a way 
that integrating over an integral number of equal increments h will 
result in y assuming its final vplue. The technique used is one 
of integrating beyond the range desired, then integrating in the 
opposite direction with a decreased value of h. The integration 
then oscillates around the end of the range, with decreasing values 
of h, until y lies within some predetermined interval containing its 
final value 0.001. The initial length of increment h is determined 
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from considerations of accuracy and time. Ihe value of h chosen 
was 0.101. For this value of h, the estimated number of steps of in¬ 
tegration is 150 and the accuracy is on the order of 10' 1J . Th e ,, alue 
of m selected is the largest possible such that none of the scaled 
derivatives exceeds range, namely m = 9, or 2“h = O.512. 


For homing on the final value y 0 of the dependent variable 

y, the value of the scaling factor 2»h is changed in the following 
manner. 



Adjustment of Interval Length 


y i+l " y e dlffers in sign from y ± 

value 




2“h 


is changed to a new 


2 V 


2m h(y 


i+l 


- y e> 


y i 


y i + l 
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The integration then proceeds with the smaller interval h' (of opposite 

sign from h) until two successive values of y again lie on opposite 

i -9 

sides of y . If, at any time \y ± - y g | < 10 the integration is 
complete. 

In the differential equations of motion of the electron 
in the gap, the independent variable 0 does not appear. The value 
of 0 is calculated by using a counter in the main routine, rather than 
integrating the equation 0' = 1. The number of equations used is thus 

4. 

Details of the Coding . The Initial Conditions Setting 
Program. Words 0 to 9L of this part of the code are devoted to 
reading a value of the parameter 0 Q from tape (using the decimal number 
input) and arranging the format and printing 0^ in both degrees and 
radians. (Print routine P 1 is used for printing). Since the scaling 
of 0 Q is by a power of 10, printing is arranged so that the decimal 
point appears in the printed results to yield the correct values for 

the unsealed quantity 9 Q . 

The remaining words are used for the following; 

(1) Set initial value of orbit index n, 

(2 ) Set initial values of p and y to zero, 

(3) Set temporary storage locations for differential 
equations routine (S5 ff•) to zero, 

(4) Calculate and set initial values of r and s. 

(The sine-cosine routine is used for the calculation), 

(5) Set initial value of 2%, 


11-19 



(6) Set initial value of y. , - y , 

J l+l •'e* 

(7) Transfer to main program. 

Details of the Coding . The Auxiliary Subroutine . The 
auxiliary subroutine is a closed routine whose function is the calcu¬ 
lation of the scaled derivatives from the values of the variables in 
accordance with the equations of motion of the electron. The variables 
V> P * and r are in locations S3, 1S3, 2S3, and 3S3, respectively. 

The scaled derivatives 2% y', etc., are to be placed in locations 
S4 to 3S4. Constants necessary, such as G, E, and negative powers 
of ten, are obtained by using the constant listing auxiliary (Library 
Routine X 3). Results of intermediate calcualtions are placed in 
temporary storage locations 1S8 through l4s8, for reference in checking 
the program for errors. The calculations are arranged so that none of 
the intermediate results exceeds capacity; otherwise the program itself 
is self-explanatory. The auxiliary subroutine uses the square root 
routine (entry from 6 l), the first word of which is in location S9- 

Details of the Coding . The Main Program . The main program 
can be subdivided into three parts; the supervisory program for inte¬ 
gration in the gap, (words OL to l8L), the external orbit recycling 
calculations (words l8L to 30L), and the orbit counting and printing 
routine (words 30L to 6lL ). 

The first part of the main program performs the following 

functions: 
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The independent variable f) is increased by h. 


( 1 ) 

( 2 ) 

(3) 

(M 

(5) 


Values for the ith step of the variables y, p, 

s, r and of (y^ - y ) are stored for future reference. 

The i+1 step of the integration is performed by 
entering the differential equations routine; i.e., 
new values of y, p, s, and r are calculated. 


y i+ ^ - y g ia calculated and its sign is compared 
with the sign of y^ - y . If the signs differ, 2 m h 
is adjusted; if the signs agree, 2 m h is left 
unchanged. 


If 


l y i + l " y e 


-9 

>10 , another step of the inte¬ 


gration is performed. If 


i y i + i - y e i 


< 10 


-9 


the 


external orbit recycling calculation is begun. 


The second part of the main program performs the calcula¬ 
tions indicated by the external orbit recycling equation in such a 

way that none of the intermediate results exceeds capacity. It should 

2 l/2 

be noted that the quantity [1 - p /g] ' is necessary for this calcu¬ 
lation. To obtain this quantity, the auxiliary subroutine is entered 
directly from the main program. The auxiliary subroutine leaves the 
desired quantity in location JS8, where it is then available for 
further calculation by the main program. It might be thought that the 
quantity could be obtained directly without re-entry from the main 
program; it should be noted, however, that the auxiliary routine was 
last used with intermediate values of the variables. For accurate 
results, it is therefore necessary to re-enter the auxiliary subroutine 
with the final values of y, p, s, and r in locations S3 through 3S3. 
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The last part of the main program, beginning with the right- 
hand order of word 30L, is devoted to orbit counting, and to calculations 
for the format of the printed results. The print routine P 1 is entered 
in such a way that decimal points are placed to correspond to the ion- 
scaled values of the results. The format used is indicated in Table 
11.4 
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degrees radians 

050.0000 O.872665 


Q 

e,n 



n 

radians 

degrees 


01. 

2.07513 

118.896 


02. 

2. 00608 

11 4. 940 


03. 

1.98559 

113.766 


04. 

2.00941 

115.131 


05. 

2.05493 

117.739 


06. 

2.07354 

118.805 

H 

07. 

2.04100 

116. 941 

H 

< 

N> 

U> 

08. 

1.99582 

114.352 

09. 

1.98125 

113.517 


10. 

2.00854 

115.081 


®i,nfl 


radians 

degrees 

$ 

V u o 

T 

1.36207 

078. 041 

0.863718 

01.98 428 

00.5^3205 

1.37242 

078.634 

0.942114 

02. 98271 

01.013640 

1.40726 

080.630 

0. 967971 

03.98364 

01.525355 

1.45802 

083.538 

0. 979635 

04.98157 

02.035536 

1.47955 

084.77 2 

0. 985871 

05.97197 

02.541871 

1.44881 

083.011 

0. 989611 

06.95861 

03.046279 

1.40483 

080.491 

0.992050 

07.95137 

03.553819 

1.39111 

079.705 

0. 99373 0 

08.95100 

04.064869 

1.41901 

081.303 

0.994930 

09.95242 

04. 576833 

1.46613 

084.003 

0.995812 

10.95028 

05. 08698 0 


Table 11.4 


Printed Results 



LOCATION 


ORDER 


NOTES 


00 5K 

3 00 F 

00 320F 

4 00 F 

00 324F 

5 00 F 

00 328 f 

6 00 F 

00 4f 

7 00 F 

00 200F 

8 00 F 

00 340F 

9 00 F 

00 64f 

K 00 F 

00 220F 

S 00 F 

00 250F 

N 00 F 

00 125F 

J 00 F 

00 7lF 

F 00 F 

00 95F 

L 00 F 

00 153F 

Auxiliary Subroutine 

00 200K 


number of equations 


11-24 



LOCATION 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


ORDER 


NOTES 


K5 F 

42 18L 


Set link 


50 IS3 

7J 1S8 

P 

2% 

2%y' = 2%? 


4o S4 

2 m hy' 



50 IS 3 

P 



75 IS 3 

66n oof 005 

8078 0000 J 

S9 F 

P 

G 


40 2S8 


P* 

!-1 

II 

0 


32 5L 


Waste 


50 5L 




26 S9 

40 3S8 


[ 1/2 - p 2 /G ] 1//2 at 3S8 

50 2S8 

7J 3S8 

40 4s8 


[ 1/2 - p 2 /G ] 3,/2 at 

4s8 

50 4s8 

75 2S3 

66n OOF 00 


s 


100 0000 0000 

J 

0.01 


7JN 40F OO 5401 



6000 0000 J 


E 1 


40 F 

50 F 

75 1S8 

00 IF 

40 1S4 

50 S8 


2%) ’ = 2 m h (20C S ) E' 

C\J 

[a 

L5 3S3 


r 
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LOCATION 

ORDER 

NOTES 

14 

66 n oof 00 1000 

0000 0000 J 

0.1 


75 1S8 


15 

4 o 2 S 4 

50 s8 * 

2^hs 1 = 2 m h (lOr) 

16 

Ll 2S3 

66 n oof 00 1000 

0000 0000 J 


17 

75 1S8 



4 o 3 S 4 

2V' = 2 I "h (-10S) 

18 

32 18 l 

22 ( )L By 0’ 

Waste 


00 220K 

Set initial conditions 

0 

4 l s8 

92 149F 


1 

92 513F 

92 961 F 


2 

4 o 5S8 

50 2L 

0 

3 

26 SJ 

L5 5S8 

-* D.N.I. Input 0 

4 

52 73F 

50 4l 


5 

26 SN 

L5 5S8 

-* Print routine Print 

6 

50 s8 

66 n 4of 00 729 



5779 5131 J 

1 -8/jt 


(degrees ) 
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LOCATION 


ORDER 


NOTES 


7 

8 

9 

10 

n 

12 

13 

14 

15 

16 

IT 

18 

19 

20 


S5 F 
40 5 S 8 
52 71F 

50 8 l 

26 SN 
92 129F 
92 513F 
19 5F 

40 6 s 8 

41 IS 3 
4l S3 
4l S5 
4l 1S5 
41 2S5 
41 3S5 
L5 5 S 8 
50 S 8 
00 2F 

L4 5S8 5 0 in A 

50 i 6 l 

26 SF 
40 7S8 

7JN OOF 00 100 
0000 0000 J 
40 2S3 
50 7S8 
7JN OOF 00 
100 0000 0000 J 
40 3S3 

L5N OOF 00 5120 
0000 0000 J 


Print 0 in radians 

Set n+l=2x2^ 

Set p to 0 
Set y to 0 

0 

Sine cosine routine 

0.01 

s = 1/200 sin 10 0 

0.01 

r = 1/200 cos 10 0 
2 “h 


Set 

starting 

values 


11-27 



LOCATION 


ORDER 


NOTES 


21 40 1S8 Set 2™h 

LIN OOF 00 10 

0000 0000 J 0.001 

22 ^0 8 S 8 Set (y 0 - y g ) 

26 SS 


MAIN ROUTINE 

00 250K 

0 L5 1S8 

10 9F 

1 L4 5S8 

4o 5S8 

2 L5 S3 

40 4S5 

3 L5 1S3 
4o 5 S 5 

4 L5 2S3 
4o 6S5 

5 L5 3S3 
40 7S5 

6 L5 8s8 
4o 9S8 

7 00 9F 
50 7L 

8 26 SL 
L5 S3 

9 LON OOF 00 10 
0000 0000 J 
40 8s8 


2 ^ 

h in A 

Replace 0 by 0 + h 

Store old values of y, p, s, r 

Store (y. - y g ) 

Differential equations Routine F 1 

y g = 0.001 


11-28 



LOCATION 

ORDER 


NOTES 

10 

36 12L 

L5 9S8 



11 

56 15L 


If y 1 - y e and y i+1 - y g agree in 


26 17L 


sign, integrate again 

12 

L5 9S8 

56 17L 



15 

L5 4S5 

y i 

if y t - y e and y i+l " y e disa e ree 


LO S3 

y i + l 

in sign, form 2 m h' 

14 

15 

40 10S8 

50 8s8 

75 1S8 

y i - y i+l 

( y i + l - y e ^ 

2% 


66 10S8 



16 

S5 F 

4o 1S8 


2 %' = 2% (y. +1 - y e ) / ( yi - y i+1 ) 

17 

L7 8s8 

LON OOF 00 

1000 J 


Is |y ± - y e l < i° 9? 

18 

36 L 

50 18L 



19 

26 S7 

L5N OOF 00 446 


Auxiliary subroutine 


6080 0000 J 


M 

20 

50 S8 

66 3s8 



21 

S5 F 

40 11S8 



22 

LIN OOF 00 

10000 0000 J 

50 S8 


10' 5 

25 

66 is 3 




S5 F 
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LOCATION 

ORDER 

NOTES 

24 

40 12S8 

50 6 S 8 (n+l)2 ~ 7 

-10 5 / P 

25 

7JN OOF 00 1000 



0000 0000 J 

40 F 

0.1 

2 6 

50 F 

75N 40F 00 1283 



1800 0000 J 

2 n/l 0 

27 

00 7F 

10 11S8 


28 

LO 12S 8 

66 n OOF 00 1000 



0000 0000 J 

0.1 

29 

SI F 



L4 5 S 8 


30 

40 13S8 

92 12 9F 

^i,n+l 

31 

92 513F 

L5 6s8 


32 

LON OIF OOF 

10 32F 

-7 

n x 2 

33 

JO 22F 

50 33L 


34 

26 SN 

L5 5S8 

Print n 

35 

52 6 if 

50 35L 


36 

26 SN 

50 N 40F,00 729 

5779 5131 J 

Print 0 in radians 

e,n 
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LOCATION 

ORDER 

NOTES 

37 

7J 5S8 



22 38L 

Waste 

38 

52 63 F 

50 38L 

Print 0 in degrees 

r e y n 

39 

26 SN 

L5 13S8 


4o 

40 5S8 

22 4lL Waste 

Set 0 for next integration 

4l 

52 6 if 

50 4lL 


42 

26 SN 

50N 40F 00729 

Print 0 i n+1 in radians 


5779 5131 J 

1.8 /n 

43 

7J 13 S 8 

22 44l Waste 


44 

52 63 F 

50 44L 


45 

26 SN 

36 46L Waste 

Print 0 i n+1 in degrees 

46 

L5 IS 3 

50 S 8 


47 

66 n oof 00 170 



4100 0000 J 

7JN OOF 00 1000 

0000 0000 J 

Calculate and print 

48 

52 71F 

50 48l 


49 

26 SN 

L5N OOF 00 70 

7106 7812 J 
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LOCATION 

ORDER 

NOTES 

50 

50 S8 



66 3S8 


51 

S5 F 

Calculate and print 


4o l4s8 

[ 1 - B 2 ] -02 

52 

52 72F 

1 


50 52 L 


55 

26 SN 



L5 14S8 


54 

LON OOF 00 100 



0000 0000 J 



40 F 

Calcualte and print energy (mev 

55 

50N 40F 30 112 



4000 0000 J 



7J F 


56 

52 82F 



50 56L 


57 

26 SN 



L5 6S8 


58 

L4n OIF OOF 

Increase n 


4o 6s8 


59 

LON 0NF OOF 



34 SK 

If n - 11. stop, read in uev. ; 

6o 

26 12SK 

If n < 11, integrate again 


00 F 


Temporary 

Storage Starting at 

Memory Location 340 

0 


0 

1 


2% 
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DOTATION 


ORDER 


NOTES 


2 

4 

J 

6 

7 

8 

9 

J.o 

11 

1:2 

1.3 

14 


1/2 - p7g 

[ 1/2 - p 2 /G ] 1/2 

[ 1/2 - p 2 /G ] 5//2 

0 

(n+1) 2' 7 
l/2 cos 10 0 

< y i + l - y e } 

< y i - y e ) 

y l - y i + l 

M [ 1/2 - p 2 /G 3 " 1 / 2 
- 10' 3 / p 

[ 1 - p 2 i - 1/2 


Al.locatlon of Memory Space . Memory space was allocated 

as follows: 

3-15 parameters 

1 . 6-63 available for list of constants for 

constant-listing auxiliary 

64-75 square root routine R 1 

74-94 decimal number input routine N 3 

95-124 sine-cosine routine T 1 

125 - 1.52 print routine P 1 
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153-193 

194-199 

200-218 

219 

220-242 

24J-249 

250-311 

320-323 

324-327 

328-331 

332-335 

336-339 

340-354 

355-940 

941-961 

962-998 

999-1023 

0-2 


differential equations routine F 1 
unused 

auxiliary routine 
unused 

initial conditions setting program 
unused 

main program 

variables for integration routine 

derivatives for integration routine 

temporary storage for integration routine 

previous values of variables for integration 
routine 

unused 

temporary storage 
unused 

constant-listing auxiliary X 3 

available for post mortem D.0.1. (Librarv 
Routine Cl) J 

decimal order input X 1 

temporary storage 


—? e Preparation. The tape was prepared in two steps . The 
non-standard parts of the program - namely the list of parameters, 
the auxiliary subroutine, the initial conditions setting program, and 
the main program —- were punched, printed, and visually checked for 
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errors. Using library routines, the tape was prepared as follows 


Decimal order input X 1 
00 9^1K 

Constant listing auxiliary X 3 

00 3K 

List of parameters 

00 64K 

Square root subroutine R 1 
00 74K 

Decimal number input routine N 3 

00 95K 

Sine-cosine routine T 1 
00 125K 

Print routine P 1 
00 153K 

Differential equations routine F 1 

00 200K 

Auxiliary routine 

00 22OK 

Initial conditions setting program 

00 250K 


Main program 
24 999N 
26 22ON 
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ko6n 


K 055 N 

K 05 N 

The tape was then checked for reperforating errors, using a tape 
comparer. The directive 24 999N was placed on the tape so that a 
correction tape could be placed in the reader before the program was 
started. Initially, however, the stop of the 24 999N directive was 
by-passed, and control was transferred to the initial conditions 
setting program by the directive 2 6 220 N. The quantities K 06 n, 

K 055 N, and K 05 N at the end of the tape are values of the parameter 
0 ^ to be read from tape by the decimal number input. 

11-6 CHECKING THE PROGRAM. The electron cyclotron pro¬ 
gram was typical in that a number of blunders in coding were mads. 

A description of the sequence of events in checking the program follows. 
When the program was first read into the Illiac, the 

quantity 

07680. 0000 00528 59832 06685 86693 4983 

was punched continually. The temporary storage was printed (P.M. 
routine C 3 ) indicating that 0 Q had been read from tape and stored 
(at 5S8). 0^ was to be initially expressed in degrees scaled by 
100 and was later to be converted to radians scaled by 10. Since the 
conversion had not occurred, the difficulty was isolated to the 
region 3 L- 7 L of the initial conditions setting program. A subsequent 
printing of the orders of the print routine (with P.M. routine C 5) 
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indicated an incorrect link and the difficulty was isolated to the 
print routine entry (word 4 l ) which was 

52 73F 50 2L rather than 52 75F 50 4L. 

When the program was read into Illiac again, noises in¬ 
dicative of integration were heard, and after an appropriate interval 
cf time punching occurred. The results were, however, not entirely 
correct. Values of 3 and were incorrect, although some small 

e, x 1 j 

comfort was gained from the fact that Q and n were printed correctly. 

The behavior of the variables f, y, p, s and r for the first ten 
steps of integration were then observed by using check point routine 
D 1 with a blocking order placed at location 256 (6 l in the main routine ) 
ahead of the entry to the Runge-Kutta routine. The initial values of 
f> V, P* s, and r were correct, exonerating for the moment the initial 
conditions setting program, but f behaved peculiarly on successive 
steps. Rather than increasing uniformly by increments of 0.011, 
l increased by 0.001, then by 0.0246, by 0.0006, and eventually de¬ 
creased slightly. It was finally found that the location 5S8 used 
for storage of ( was also being used for temporary storage by the 
auxiliary routine. This blunder was cured by replacing the orders 
40 5S8 50 5S8 in words 10 and 11L of the auxiliary routine by 

40F 50F. 

With a somewhat longer correction tape, the program was 
again checked on the Illiac . The results of the first integration and 
orbit cycle were correct, but the results for the second orbit were 
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all incorrect with the exception of the orbit index n. The re-entry 
into the initial conditions resetting program was checked and it was 
found that p was incorrectly being reset to zero; i.e., the program 
was robbing the electron of all the velocity it had gained during its 
first orbit. The correction tape increased in length; word 6lL of the 
main program became 26 12SK 00 F rather than 22 11SK 00 F. The 
orders for clearing ^ and p were interchanged in words 11L and 14L 
of the initial conditions setting program. 

After the following code check, the results for the first 
orbit were again incorrect. After mutterings of "something is wrong 
with the computer", it was discovered that a terminating symbol had 
been omitted on the correction tape. The correction process then 
converged and correct results were obtained. 

11 *7 CONCLUSION. The preparation of a problem for 

solution on a digital computer is by no means a completely objective 
process. No two programmers would prepare the same problem for solu¬ 
tion in an identical way. Furthermore, a second coding of the problem 
by the same programmer would differ from the first preparation. 

The personal preferences of the programmer affected the 
electron cyclotron program in a number of ways. Memory space other 
than location 0, 1, and 2 was assigned for temporary storage; it 
was felt that such an arrangement might aid in code checking. The 
S terminating symbol was used for designation of the locations of 
first words of subroutines, so that the detailed coding could be 
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completed before memory locations were assigned to the subroutines. 
Thus, the location of a subroutine affects only one of the S para¬ 
meters rather than the addresses of a number of orders in the program. 

Among the changes which might be made if the program were 
rewritten is a change in the handling of the variable 9. If 9 were 
expressed in revolutions, the overflow properties of the computer 
would simplify the recycling computation. 
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CHAPTER 12 


THE CATHODE RAY TUBE DISPLAY 

12.1 GENERAL DESCRIPTION . Results of certain types of 
calculations may He most conveniently obtained from the computer by 
use of the cathode ray tube display. Whenever results of a cal¬ 
culation can be represented pictorially as a graph or diagram, the 
cathode ray tube output from the computer provides a rapid and 
elegant method for obtaining these results. If results would other¬ 
wise be graphed manually the cathode ray tube saves the human time 
required for this additional processing and also reduces the machine 
time required to present the results to the user. 

Another use of the cathode ray tube is to provide inter¬ 
mediate results in such a form that they may be analyzed subjectively 
during the course of a program. In this way the programmer may 
discover errors in his program at an early stage in the calculation. 
He may also obtain immediate results which may be used to determine 
the way the remainder of the program is to be run. Even if the 
final results must be presented with greater accuracy than is 
possible using a graphical display, some advantage may be obtained 
by supplementing the digital results with a graph or diagram. The 
cathode ray tube output is called into use by turning the output 
switch on the input-output rack to the cathode ray tube position. 

When the switch is in this position the cathode ray tube is the only 
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form of output available to the user. The display on the cathode 
ray tube is limited to a 3" x 3" square area centered on the face 
of the tube. Within this square area, chosen points may be brightened 
on a 256 x 256 regular square raster. A finer raster is unnecessary 
because of the limited resolution of the cathode ray tube. This 
limited resolution makes the display of digital information such 
as numbers or figures somewhat unhandy. It has been found that no 
more than about 200 fully legible numbers or letters may be dis¬ 
played in a single frame. In general the requirements of the format 
will restrict it to somewhat less than this. This disadvantage, 
however, is somewhat offset by the possibility of using characters 
having any shape which the programmer desires and by the greater 
speed of output. A subroutine has been written which will display 
characters at the rate of about 4-5 per second. 

In order to retain a permanent record of displayed re¬ 
sults it is necessary to photograph the face of the cathode ray 
tube. Two cathode ray tubes are installed in the rack. They 
are driven so they operate in parallel and hence display the same 
results. One is used for visual observation and the other is 
equipped with a semi-automatic camera to permit photographing 
the display. The film advance mechanism is controlled by the 

program while insertion, removal, and development of film must 
be done by hand. 
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12.2 ORDERS CONTROLLING THE CATHODE RAY TUBE DISPLAY . 

Orders which would normally cause characters to be punched or printed 
will affect the cathode ray tube when the output switch is in the 
cathode ray tube position. The conventional output order to use for 
cathode ray tube display is 82 l6P. When this order is executed 
it will cause one spot to be brightened on the face of the cathode 
ray tube. The position of the brightened spot depends upon the 
contents of the l6 leftmost binary digits 2° - 2 of the accumulator. 
Let these digits be designated by the symbols a^, a^, • • • > a i5' 
Assume an origin or coordinates in the lower left-hand corner of 
the square raster and let d represent the length of one side of 
the square raster (d = 3")- Then the coordinates of the brightened 
spot will be: 

—8 

y = d (a Q 2 1 + a x 2 + . . . + a^ 2 ) 

8 

x = d (a 8 2" 1 + a g 2 + • * * + a 15 2 ) 

These formulae say, in effect, that the ordinate is obtained from 
the first eight binary digits of A, regarded as a number, and that 
the abscissa is obtained in the same way from the second eight digits. 

During the execution of the 82 l6F order the AQ register 
will suffer a left shift of l6 places. This output order takes 800 
microseconds. A variation of this order permits a spot to be 
brightened in only 400 microseconds. When an 82 8F order is executed 
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the ordinate is determined in the same manner as it is for an 82 l6F 
order, y = d(a^ 2 +a^2^+. . .+a^2^), but the abcissa will 

be the same as that of the spot produced by the last 82 l6F order 
which was executed. In this way time may be saved when several points 
having the same abscissa are displayed. 

A letter output order is used for operating the semi-automatic 
camera. The following order is used for this purpose: 

92 769 Advance the film one frame. 

The camera has no shutter, so that programs having long 
periods of computation (say 15 minutes) between displays of results 
should be arranged to advance the film before and after each display 
in order that the results will not be fogged from filament light in 
the cathode ray tube. 

Exposure of a single frame of the film takes place during 
the period of time that the spots that make up the picture are 
brightened. Amount of exposure is therefore controlled by the 
intensity setting of the cathode ray tube. Exposure of a spot may 
also be increased by multiple brightenings. If calculations are 
made while the frame is being exposed the resulting photograph will 
in no way be affected although the image on the visual cathode ray 
tube may fade partially before the completion of the frame. Ad¬ 
vancing the film one frame takes one second and the Illiac will 
wait until this operation is completed before executing other out- 
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put. orders • If the orders following the film advance do not in¬ 
volve output, their execution will be begun after only 200 microseconds. 

No output orders other than those mentioned are of practical 
use in the operation of the cathode ray tube display. A 92 order 
having an odd address will affect the camera mechanism, while any 
other output order will brighten a spot. 

12.3 PROGRAMMING FOR THE CATHODE RAY TUBE DISPLAY. As 
with typewritten output, the cathode ray tube display is usually 
programmed by means of subroutines. If digital results are desired 
the n umb er of orders concerned with output will be greater when the 
cathode ray tube display is used than when the punch is used. In¬ 
formation needed to form a character (such as a letter or decimal 
digit) may be stored in a single word in the memory, but the 
mechanism required to decipher this information rapidly may take 
as many as 35 words. Such a device is used in library program 
CO 1. This program uses 63 words to display information similar 
to that printed by library program C 3 which has only 31 words. 

Graphical display is also most conveniently accomplished 
by means of subroutines. Program 0 1 is an example of a subroutine 
which enables one to plot axes and points on the cathode ray tube. 

Any method of graphical display will almost inevitably require 
scaling of the coordinates of the points to be displayed. Since 
it is desirable to utilize the full screen area of the cathode 
ray tube, up scaling may be necessary in some cases instead of down 
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scaling. The full screen is regarded as haring a range of 1 for 
both coordinates when Library Routine 0 1 is used. The full range 
of all variables in this case should be scaled so as to correspond 
to 1. If more than one graph is to be displayed on a single frame, 
however, it is sometimes better to avoid entering the point plotting 
subroutines more than once. Consider the following examples The 
points (x^, y^) and (x^, y^) are both to be displayed relative to 
the origin. Let us assume that axes have already been displayed 
by program 0 1. The quantities x Q and y Q have been supplied to 
program 0 1 when the axes were plotted so as to locate the origin. 
They represent the x and y coordinates respectively of the center 
of the screen in the coordinate system used. Let the necessary 
parameters be those given in Table 12.1: 


12-6 



LOCATION 


PARAMETER 


10 

11 

12 

13 

lb 

15 

16 

17 


v Assume these coordinates have been 

x properly scaled before being stored. 

y l 

y 2 

scaling factor for y^ 
scaling factor for y 2 
1/2 

Table 12.1 

Parameters for Cathode Ray Tube Display 


A program to plot the two points might be the one given in Table 
12 . 2 : 


p 

50 13 

7J 15 

scale y^ 

P+1 

JO 12 

5 o p+l 

plot (x x , y x ) by use of subroutine 

p+2 

26 (to 

program 0 1) 


50 lb 

scale y 2 

p+3 

7J 16 

prepare 2(y 2 + l/2 - y^) 


Lb 17 

p+4 

LO 11 

00 1 

to be displayed 

P+5 

82 8 

plot (x x , y g ) without use of a 
subroutine 

Table 12.2 


Program for Cathode Ray Tube Display 
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In the program of Table 12.2 the second ordinate y 2 is first scaled 
and then translated so as to be present in the accumulator in the 
proper form for display purposes, it is then displayed by use of the 
82 8 order, thus saving machine time and also saving the programmer 
the trouble of re-introducing the abscissa x . 
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CHAPTER 13 


PROGRAMMED CHECKS 

13.1 GENERAL PRINCIPLES OF CHECKING . Any person using 
the Illiac as an aid to his research is likely to question the 
reliability of results of calculations carried out by the machine. 
Although the Illiac is quite reliable in its operation, one is 
almost never justified in using Illiac results directly for publi¬ 
cation or for further research without assuring oneself in some 
way of their correctness. If the calculations carried out by the 
computer are to be treated with the same care and thoroughness 
that is usually required of results obtained in other ways some 
method of checking must be found. Naturally the Computer Laboratory 
cannot guarantee the correct operation of the computer and it is 
the responsibility of the person preparing the problem to provide 

adequate checking of his results. 

The subject of checking is by its nature, a difficult 
one to treat in generality since the most suitable checking method 
usually depends upon the nature of the calculations. An obvious 
checking procedure is available in some cases and not in others. 

For example, one may check the solution of a set of algebraic 
equations by substituting the results into the original equations. 
Whenever a checking process of this sort exists it is preferable 
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to checking methods which only check some of the operations Involved 
in the calculation since it provides an overall check on the operation 
of the computer from the beginning to the end of the problem. If the 
results of a calculation may be plotted as a graph one may often verify 
the correctness of the results by noting whether the points follow a 
smooth curve. Such a simple check as this may be used only when the 
various points are calculated Independently and should then be used In 
conjunction with other checking processes. Some types of calculation 
may be closely approximated by much simpler formulae which may be 
evaluated independently. When such approximations exist they are 
Often adequate for checking purposes. A checking process must thus 

be tailored to fit the problem and its effectiveness will often depend 
upon the ingenuity of the programmer. 

Let us assume that a program is divided into distinct runs 
of length t x and that a checking period of length t £ follows each 
run. If a malfunction or machine error occurs either during the run 
or during the checking period it will be detected at the end of the 
checking period and the entire run will have to be repeated. The 
average time wasted due to errors and checking and the optimum length 
of run may be calculated using these assumptions. Take the time T as 
the mean time between machine errors, i.e. the reciprocal of the error 
frequency. It is not strictly correct to speak of a fixed time T 
since the error frequency depends upon the nature of the calculation 
and other time dependent factors. Nevertheless, we shall assume a 
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fixed T and show later that our conclusions do not depend upon an 
accurate knowledge of T. A reasonable choice might be T = 8 hours. 

The average useful time per run may be taken as the prob¬ 
ability that no error will occur during the run multiplied by the 
length of run. Assuming a Poisson distribution of errors one obtains 


Average useful time per run = t exp[-(t + t )/T] 

X X 2 

Since the total time taken by the run is t ± + t we find that the 
fraction of the total time which is useful is 


F 



exp [-(t x + t 2 )/T] . 


( 2 ) 


The first factor differs from 1 because of the time taken by the 
checking process while the second factor differs from 1 because 
of the likelihood of errors occurring during the run. If we 
maximize- (2 ) with respect to t^ we obtain 

t l + t l t 2 - t 2 T = 0 » 


(3) 


giving 


t ] _/t 2 V1/4 + T/t 2 - 1/2 . 


( 4 ) 


If we assume that T » t 2 we obtain the approximation 

t ± =Vv? . 


(5) 
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If the same assumption is made in formula (2 ) we obtain for the 
fraction of the total time which is useful 

F = 1 - v\ - V T • 

Substituting (5) in (6) we obtain 
F = 1 - 

as the largest fraction obtainable. It is important to notice 
that if T is much larger than t formula (6) will remain nearly 
1 over a wide range of t . In effect, this means that our selec¬ 
tion of t 1 is not critical provided it is much more than t and 
much less than T. For example, if the time for checking, tg, is 
5 minutes we obtain for the optimum length of run 

= "\/5 min. x 480 min. = 4-9 minutes . 

Yet as seen from formula (6), if the length of run were doubled 
or halved the wasted time 1 - F would not be much greater. 

In formula (l) if t^ is allowed to increase we see that 
the function reaches a maximum when t^ = T. Therefore if a run 
is permitted to exceed T the useful time achieved by the run (on 
the average) begins to decrease. Under no circumstances will any 
advantage be gained, therefore, by making t^ greater than T. From 
(4) we see that this limit is approached as tg increases. 



An alternative system of checking may be used in which 
the checking process is incorporated in the calculation and the 
program is arranged to stop whenever an error occurs. The run must 
be restarted from the beginning after the error. Each successful 
run (of length t 1 ) will be followed by a period of printing or storing 
intermediate results which will allow the following run to be re¬ 
started at that point. We shall assume this period to be of length tg. 
The difference between this method and the preceding one is that an 
unsuccessful run is cut short in this case. Although formula (l) is 
unaltered the ones which follow are changed somewhat. This is be¬ 
cause the total time taken by a run (on the average) is no longer 


t^ + t but is rather 

(t x + t 2 ) exp[-(t 1 + tg )/T] + 




U exp[- (u/T ) ] dt 



exp[ - (t 



The first term on the left is due to successful runs while the second 
term is due to unsuccessful runs. The formulae corresponding to the 
previous analysis are 
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F = 1 -V 2t ? /T . " / 

Here we see that the optimum length of run is increased by a factor 
ofV2~^ and that the wasted time is decreased by the same factor. 
Although no closed expression can be obtained corresponding to (4) 
we note that the remarks concerning in the neighborhood of T 
still apply. Qualitatively the behavior of F as a function of t^ 
is still the same so F remains nearly unity when % « t « T. 

Thus a careful choice of t^ and an exact knowledge of T are still 
not important. 

■*-3-2 INPUT CHECKING . Often, an overall check on the 
entire program is not possible and individual checks must be made 
on separate operations. The operations of reading and punching 
are more subject to error than other operations and are also easier 
to check in most cases. Standard techniques are also available 
for checking input and output operations so that relatively little 
burden is placed on the programmer. 

The sum check is commonly used to check input-output 
operations. It derives its name from the fact that a primitive 
form of the sum check consists of merely adding together all the 
words in the memory before the calculation begins, but after the 
program has been read into the computer. If this sum agrees with 
a previously obtained sum, one concludes that the program has been 
placed in the memory in exactly the same way as it was when the 
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sum was originally obtained. Agreement between the sums is taken 
to mean that no error was committed during the reading of the pro¬ 
gram. 

A sum check is particularly useful for checking the reading 
of a master tape which Is to be used many times. Such tapes eventually 
wear out and a sum check provides an indication of this as well as a 
check on the operation of the computer. Library Program X 7 is an 
interlude which computes the sum of the memory and compares the 
result with a sum which is especially prepared on tape the first 
time the tape is used. A somewhat more reliable check than would 
be obtained by simply summing is obtained in this program by per¬ 
forming a circular right shift of the accumulated sum during the 
summing process. By this technique one may avoid the possibility 
that two errors occurring in the same digit position, in different 
words may compensate, thus yielding a correct sum. Misplacement 
of words is also detected by this technique. A typical interlude 
to sum the memory might be the one in Table 13*1* 
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0 

1 

2 

3 

4 

5 

6 

7 

8 

9 


L 5 

F — 

10 

IF 

s 4 

F 

L4 

(i)f~ 

4o 

F 

F 5 

1 L - 

4o 

1L 

L0 

8 l = 

36 

L 

81 

40F 

L 0 

F 

40 

F 

L3 

F _ 

36 

999F 

FF 

F 

26 

999F 

3^ 

F -1 

L4 

1024fJ 

00 

F -1 

26 

L J 


26 in 


Place accumulated sum in A and 
perform circular right shift. 

Add next word to sum and restore 
in zero. 


Step address to select next word. 


Test for the end of the summing 
process. 


Read previous memory sum from tape 
and compare with memory sum obtained 
this time. 

If sums do not agree stop on an FF 
order; otherwise return to the D.O.I. 


Test constant. 


Enter interlude. 


•Table 13.1 
Sum Memory Interlude 

After computing the ''memory sum" this interlude will read a previously 
prepared word from the tape which it compares with the memory sum. If 
the two disagree the routine will stop on an FF order. This order has 
been specially chosen to indicate the failure of nrosramn^ 
should be used whenever a programmed check 

stop so the computer operator will be able 
the stop. 


is made instead of the OF 
to identify the nature of 
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So far we have discussed a method for insuring the correct 

reading of a program tape or master tape. A data tape which is to "be 
used many times may be treated in the same way. Tapes which are pre¬ 
pared by the computer constitute a special problem and may be best 
discussed in connection with output checking. There still remains the 
problem of checking the reading of data tapes which are prepared by 
hand and which need to be read only once. Wo entirely satisfactory 
method for checking such tapes has been developed which does not 
involve reading the tape twice. If arrangements are made for the 
master program to read a tape twice one may store the numbers the 
first time they are read and then compare them directly on the second 
reading. If memory space does not permit storing the entire set of 
numbers which are to be read one may effect an equally valid check 
by computing a sum of the numbers as they are read the first time 
and retaining only this sum to compare with a similar sum computed 
on the second reading. 

Other ways of checking data tapes involve special preparation 
of the tape. One such technique is to follow the punching of each 
number by the sum of its digits modulo 9- These sums must be com¬ 
puted by hand as the tape is prepared thus complicating the preparation 
of the tape and increasing the likelihood of human errors. As the 
tape is read the corresponding sums may be formed by the computer (a 
special input routine is required) and the reading of the tape checked. 
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A somewhat simpler technique requires the use of a 
code in the punching of characters. Punching in code may be 
simplified if one covers the Teletype keys with caps which 
indicate the coded characters while the tape is being punched. 

A convenient code uses the 5th hole of a character as a parity 
check on the other four; for example, the first four holes might 
be used to indicate a sexadecimal character as they are at present 
and the 5th hole would be punched or not punched so as to make the 
number of holes per character odd. A special input routine would 
be used to interpret the coded digits and check for errors. An 
odd number of errors occurring in the reading of a single charac¬ 
ter would be detected by means of such a check since the parity 
would then be altered. 

In order to obtain a typewritten copy of the coded tape 
It is necessary to clamp the 5th pin on the transmitter-distributor 
unit of the Teletype in the down position so as to suppress the effect 
of the 5th hole. Thus all carriage returns must be inserted manually 
as the typewritten copy is formed. Another drawback to this method 
of checking is that a character may be ignored by the reader or read 
as two identical characters without an indication of error. In 
order to check that the correct number of characters has been read 
one must include some additional restriction on the typing of the 
data. Such a restriction might be that all numbers which are typed 
should have an even number of digits. Non-significant digits would 
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then have to he included in some numbers thus increasing the time 


spent typing as well as the time required to pass the tape through 
the reader. Still the time would be less than that required to read 
the tape twice. 

13»3 OUTPUT CHECKING . The checking of punched results 
is usually a simple operation compared to the checking of the reading 
of hand-prepared tapes. A standard technique is to arrange to read 
the output tape back into the computer at the completion of the problem. 
If the punched results are retained in the memory they may be compared 
directly at this time. If it is not possible to store the punched 
results one may compute a sum of all results as they are punched and 
then recompute this sum as the output tape is read back into the 
computer for checking purposes. 

If punched output from the computer is intended to be read 
back into the computer at a later time it is convenient to include a 
sum of the punched results on the output tape. This sum may be com¬ 
puted during the punching process and finally punched at the end of 
the tape. When this tape is read back again it may be checked for 
punch or reader errors by recomputing the sum as it is read and com¬ 
paring this sum with the sum punched at the end of the tape. 

No convenient method exists for checking the Teletypewriter 
or the cathode ray tube output. This is because no input device exists 
which allows one to read printed or photographed results into the 
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computer. Three methods exist for checking such results but all 
three involve considerable human labor which must be done without 
error if a check is to be made. These methods are. 

1. Print the results twice from a checked tape 
and compare the two printed copies by eye. 

2. Copy the printed results by hand onto Tele¬ 
type tape and then use the sum check method 
of the previous paragraph to check this tape, 
treating it as though it were computer output. 

3- Include some redundancy in the printed results 
so that the output may be checked by hand 
computation. One might, for example, follow 
the printing of each number by the sum of its 
digits modulo 9. This computation could then 
be made by hand using the printed results and 
could be compared. 

13.4 WMORY CHECimo. As we have seen in the previous 
sections of this chapter, the most satisfactory check of a calculation 
is an overall check which checks all operations performed during the 
calculation simultaneously. When such checks are Impossible, however, 
one must resort to the use of a combination of partial checks. One 
then wishes to check the operation of the input-output equipment, 
the memory, and the arithmetic unit by means of independent checks. 

The checking of the high-speed memory is also usually treated separately 
from the checking of the drum memory. 
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During a long calculation it sometimes happens that words 
in the high-speed memory become altered for no apparent reason. If 
the altered word happens to be an order pair in the program the result 
will usually be a malfunction of the program which in many cases will 
not be difficult to detect. Yet there are many memory failures which 
can produce subtle changes in the results, and for this reason a memory 
check must be made. Some methods of arithmetic checking are complete 
enough to include the memory and no special check of the memory is 
needed,but if this is not the case some special method must be devised. 

Library program X 10 is a general subroutine for checking 
all parts of the memory which are not altered after the first few 
minutes of calculation. Provision is also made for removing the 
program being checked from the computer each time such a check is 
made. Complete input checking is also incorporated in this program 
as an additional feature. This versatile subroutine may either be 
arranged so that the removal of the program from the machine is under 
the control of the operator or it may occur at intervals determined by 
the program. When removal is to occur the master program must enter 
X 10 as a closed subroutine. It causes a tape to be punched con¬ 
taining information concerning the subroutine link, those words in 
the memory which may suffer modification, and finally the memory sum. 

A special tape is read at this time to specify the words which may 
be modified. After punching is completed the program stops on an OF 
order. To resume the program at a later time one merely reads the 
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initial tapes as if the program were being started from the 
beginning and then arrests its progress with the white switch. 
Using a bootstrap start the tape which was punched just prior 
to removal may be read. All words will be restored to their 
condition at the time of entry to X 10 and a sum check performed. 
If the sum check agrees with that punched previously one is 
assured that words have been read correctly either at this 
reading or the previous one, that the specially punched tape 
was punched correctly, and that no word in the memory suffered 
a change during the previous run with the exception of those 
words specially designated. These latter words are usually 
checked by the arithmetic check. 

Checking of the drum memory is achieved more easily 
than checking of the high-speed memory. In most applications of 
the drum memory words are transferred to and from the drum in 
blocks of perhaps ten or more. The standard block transfer 
routine effects a check on the drum by recording a sum check 
following each block which is recorded on the drum. Upon play¬ 
back from the drum a similar sum is computed to be compared 
with the one which was recorded. In this way the operation of 
the drum is completely checked without special attention on the 
part of the programmer. Only when individual words are to be 
transferred to and from the drum must special arrangements be 
made for checking. 
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13.5 ARITHMETIC CHECKING. 


In most cases a characteristic 


identity of the problem being solved is used for arithmetic checking 
purposes . Two examples are presented here to show how arithmetic 
checking is carried out. 

a) Assume that n values of cos ” k with k = 0, 1* ..., n 

are to be used in a calculation and that — is non-integral. Checking 

n -1 2jfm 

may be accomplished by means of the identity k Z Q cos — k = 0. As 

the cosine functions are computed they may be summed and the final 

result tested to check that it is sufficiently close to zero. 

b ) Assume that the n products ax i are to be formed where 

a is a given number and x ± is given for i = 0, 1, ..., n - 1. These 

n -1 n -1 

products may be checked by computing i Z Q x ± and ax. ± so as 

n -1 n -1 


- 1 


to allow a check by the identity 


ax. 


In case a convenient identity cannot be found one may 
sometimes avoid duplicating all calculations by using an approximate 
calculation as a check. An integration process for example, may be 
checked by carrying out the same process using a larger step (or 
mesh size) during the checking run. The checking run will thus take 
considerably less time than the original, yet it will detect any 
malfunction which produces a gross error. 

13.6 SUMMARY- Although the operation of the Illiac is 
reliable by computer standards one is rarely justified in using Illiac 
results for scientific purposes without checking them in some way. 
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Sometimes it is possible to check results visually or by hand, but 
when this is not possible one must make use of a programmed checking 
technique. As has been seen the simplest and most effective pro¬ 
grammed checks are those which check the solution of the problem 
as a whole rather than the individual operations. Such overall 
checks are possible when the solution may be substituted back into 
an equation or when some fundamental redundancy exists in the results 
Even when an overall check is not possible one may obtain 
a valid check by combining checks on input-output, memory and arith¬ 
metic units. Although such individual checks are often difficult to 
perform they are always preferable to running the problem twice and 
thereby doubling the time taken. A summary of methods of checking 
individual operations is given in Table I 3 . 2 . 
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Operation Checked 

1. Reading of program tape or 
reading of data tape which 
is to be used several times. 

2. Reading of computer output. 

3. Reading of data tape to be 
used once. 


4. Punching of output tape. 


5 . Printing of results. 


6 . Storage in high-speed memory. 

7 . Transfer to and from the drum 
and storage on the drum. 

8 . Arithmetic operations. 


Checking Methods 

Sum check such as Library 

Routine X 7 or X 10. 

Sum check to be punched on 

computer output tape. 

a. Read tape twice and arrange 
that one reading shall check 
the other. 

b. Use a parity check code on 
characters which are read. 

c. Use a modulo check on the 
digits of each number. 

a. Compare with stored results 
by reading punched results. 

b. Compute sum check as punched 
tape is read and compare 
with previously computed and 
stored or punched sum check. 

a. Print twice and compare by 
visual inspection. 

b. Prepare tape by hand from 
printed page and check as 
in 4b. 

c. Use a modulo check on printed 
numbers and compute by hand. 

a. Use a sum check routine such 
as X 10. 

b. Combine with arithmetic check. 

Use a sum check such as the one 
included in N 7> 

a. Use an identity which relates 
the various quantities computed 

b. Use an approximate method which 
is faster than the method being 
used in the actual calculation. 


Table 13.2 

Checking Methods for Various Operations 
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CHAPTER 14 


THE PROGRAM LIBRARY 

The problem of planning and carrying out a large-scale 
computation is often almost entirely a matter of organization in 
which the detailed parts of the computation may be subcomputations 
which are common to many problems. Consequently, it is of great 
importance for the efficient use of a large-scale computer to have 
readily available to the programmer means for carrying out these 
subcomputations. For example, it would be very wasteful if each 
programmer had to prepare his own routine for converting binary 
numbers in the machine to decimal numbers. 

The Illiac library is a collection of routines which have 
been prepared to make programming easier. They have all been 
machine tested and copies punched on tape are kept available for 
use. All are prepared for use with the Decimal Order Input (See 
Chapters 4 and 5). 

l4.1 LIBRARY ORGANIZATION . Each Illiac routine bears 
a label consisting of a letter (or possibly two letters) followed 
by a number. The library has been divided into categories according 
to the types of operations performed, and the letter in a label 
denotes the category to which the routine belongs. Occasionally a 
routine, though belonging to one category, is useful only in 
connection with operations in another category (as, for example, a 
floating decimal auxiliary), and in such cases the letter of the 
second category has been added to the label. The number in the label 
is simply a serial number within the category. 
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An accession serial number is also to be found on each 
routine; this is a Laboratory reference number and has no connection 
with the categories. 

li+ ' 2 PROGRAM DESCRIPTIONS AND SUMMARY SHEETS , Each library 
routine is written up in two ways. One is a complete description of 
the method used with the list of orders of the routine. The other is 
a symmary sheet which gives enough information to enable the routine 
to be used correctly and to be compared for usefulness with another 
routine of the same type. Programmers can obtain copies of summary 
sheets by inquiring in the Teletype Room. Complete program descriptions 
are available for reference in the Teletype Room and a user who must 
make frequent reference to a particular program or who needs to modify 
it for his own use may obtain a copy from the Teletype Operator. 

.LIBRARY TAPES. All library programs are prepared for 
use with the Decimal Order Input (Library Routine X l) and each is 
punched on heavy gray Teletype tape and filed in a cabinet accessible 
to all programmers. To use a library program the programmer need 
only help himself to the tape, copy it himself or ask the Teletype • 
Operator to assist him, and replace the library copy. 

ILLIAC LIBRARY CA TEGORIES . The categories of the 
Illiac library are as follows: 

Programmed Arithmetic 

A'. Floating Point 

B. Other Programmed Arithmetic 
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Code Checking 


c. 

Post mortem checks 

D. 

Dynamic code checks 

Integration 

E. 

Quadrature 

F. 

Ordinary differential equations 

G. 

Partial differential equations 

Operations 

on Functions 

H. 

Zeros and minima 

I. 

Interpolation 

J. 

Operations on polynomials and power series 

K. 

Approximations and statistics 

Linear Algebra 

L. 

Simultaneous linear equations 

M. 

Other operations on matrices and vectors 

Input and Output 

N. 

Number input 

0. 

Scope output 

P. 

Printing and punching 

Mathematical Logic 

Q- 

Mathematical logic 

Particular 

Functions 

R. 

Roots and fractional powers 

s. 

Logarithmic, exponential and hyperbolic functions 

T. 

Trigonometrical functions 

V. 

Other special functions 

Organizational 

w. 

Countingj sorting and selecting 

x. 

Program preparation 

Miscellaneous 

z. 

Miscellaneous complete programs 



1^-5 ILLIAC LIBRARY Library routines which are of most 


interest or most heavily used are placed in the active library, while 
seldom-used routines are placed in the inactive library. On April 1, 
1955# the active Illiac library contained the following routines: 


LABEL 

SERIAL 

A1 

63 

A2 

128 

A3 

125 

A4 

87 

A5 

138 

A6 

154 

Cl 

48 

C3 

94 

C4 

97 

C5 

103 

c6 

146 

C7 

163 

C8 

164 

C9 

174 

CIO 

175 

Cll 

176 

D1 

95 

d4 

70 

E2 

158 

FI 

114 

F2 

115 

F3 

129 

FA1 

122 

G1 

78 

G2 

98 

HI 

71 

H2 

72 


TITLE 


Floating Decimal Arithmetic Routine 
Floating Decimal with Standard Auxiliary 
Convert a Number from Floating Decimal 
Representation to a Normal Machine Form 
1.7 Precision Floating Binary Arithmetic 
with Floating Decimal Conversion 
Complex Number Arithmetic 
Floating Decimal Routine and Auxiliaries 

Post Mortem of the Decimal Order Input 

Post Mortem for Fractions and their Locations 

Post Mortem for Integers and their Locations 

Print Sets of Order Pairs and their Locations 

Address Search Routine 

Drum Sexadecimal Post Mortem 

Drum Fraction Post Mortem 

Print Sets of Floating Decimal Numbers and 

their Decimal Locations 

Combined Post Mortem 

Drum Order Pair Post Mortem 

Check Point Routine 
Control Transfer Check 

Integration by Simpson's Rule 

Solution of a System of Ordinary Differential 
Equations 

Solution of a System of Differential Equations 
by Milne's Method 

Integration of n Simultaneous Second Order 
Differential Equations with Initial Conditions 
Specified 

Second Order Linear Differential Equation with 
Two Point Boundary Condition 

LaPlace's Equation - Liebmann Method 
Poisson’s Equation - Liebmann-Frankel Method 

Inverse Interpolation 

A Search for the Real Roots of f(x) = 0 
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LABEL 

SERIAL 

TITLE 

H5 

80 

Minimization of a Function of Two Variables 

H4 

84 

Minimization of a Function of Four Variables 

H5 

85 

Minimization of a Function of n Variables 

H6 

86 

Minimization of a Function of n Variables 

Treating One Variable at a Time 

HP1 

111 

Zero of a Solution of a Differential Equation 

11 

67 

Interpolation 

JJ. 

159 

Roots of a Polynomial 

K2 

135 

Product Moment Correlations, Means, Standard 
Deviations, Variances and Covariances 

K5 

137 

Least Squares 

K5 

182 

Autoc orr elations 

LI 

112 

Solution of a Set of Simultaneous Linear Algebraic 
Equations 

L2 

113 

Automatic Linear Equation Solver 

L5 

100 

Complete Linear Equation Solver 

L4 

166 

Half-precision Solution of a Set of n Simultaneous 
Linear Algebraic Equations 

L5 

165 

Half-precision Automatic Linear Equation Solver 

MO 

l4i 

Eigenvalues and Eigenvectors of a Symmetric Matrix 

M2 

118 

Automatic Inversion of a Symmetric Matrix 

M3 

117 

Automatic Eigenvalue-Eigenvector Program 

M4 

136 

Closed Eigenvalues and Eigenvectors 

M5 

139 

Solution of the Determinantal Equation j A - XB | = C 

m6 

145 

New Matrix Multiplication 

MLO 

153 

Estimation of Communalities by the Method of 

Maximum Likelihood 

Mil 

160 

Matrix Multiplication (Closed) 

M12 

173 

Triangularization of a Matrix and its Determinant 

M15 

183 

Linear Programming 

MAI 

104 

Matrix Multiplication with Floating Decimal Auxiliary 

N1 

6 l 

Input One Number from Tape, Integer or Fraction 

N2 

88 

Input a Sequence of Decimal Fractions 

N5 

167 

Drum Fraction Input 

n6 

168 

Drum Integer Input 

N7 

172 

Transfer Blocks of Words from the Memory to the 

Drum or from the Drum to the Memory 

N8 

181 

Read One Number from Tape as Integer or Fraction 

02 

143 

Routine to Plot Points and Axes 

03 

147 

Display Numbers and Letters on the Cachode Ray Tube 
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LABEL 

SERIAL 

04 

l6l 

05 

162 

06 

177 

Pi 

75 

P2 

52 

P3 

53 

P4 

55 

P5 

37 

p6 

25 

P10 

133 

Pii 

148 

R1 

116 

R2 

105 

R3 

106 

RA1 

92 

RA2 

99 

SI 

^3 

S2 

64 

S3 

130 

SA2 

127 

SA3 

155 

T4 

l4o 

T5 

157 

TA1 

126 

TA2 

156 

VI 

82 

V2 

120 

V4 

178 

V5 

184 


TITLE 


Fast Fraction Display 4x6 Raster 

Fast Character Display 4x6 Raster 

Punch to Cathode Ray Tube Conversion Program 

Print One Number Fractional or Integer in a 
Manner Determined by a Program Parameter 
Print (A) with or without Sign to n Places as 
Determined by a Program Parameter 

Print n Digits of an Integer with or without 
a Sign 

Zero Suppression Integer Print 

Print One Number in a Parameter Set Layout 

Single Column Print 

Print Headings Interlude 

Hounded or Unrounded Print 

Square Root Routine 
Integral Root Al/t> 

Fractional Power Routine 

Floating Decimal Square Root Auxiliary 

Floating Decimal Cube Root Auxiliary 

Natural Logarithm 

Exponential 

Logarithm 

Exponential Auxiliary for Floating Decimal 
Natural Logarithm 

New Arctan X Routine 
Sine Cosine 

Sine Auxiliary for Floating Decimal 

Arc Tangent Auxiliary to the Floating Point Routine 

Legendre Polynomials 
Tchebyscheff Polynomials 
Fourier Analysis 
Spherical Bessel Functions 


XI 

X3 

X7 

X8 

X9 

X10 

XA1 


18 

90 

142 

144 

169 

170 
121 


Decimal Order Input 
Constant Listing Auxiliary 
Sum Check 
Assemble Blocks 

Program Interruption Routine (Without Sum Check) 
Program Interruption Routine with Sum Check 
Floating Point Constant Listing Auxiliary 



CHAPTER 15 


DEFINITION OF TERMS 


* A 2-5 

The accumulator register or the contents of the accumulator 
register. 

Access Time 

The time to get a number from the memory to a register or the 
time to get a number from a register to the memory. 

Accumulator 2-5 

The register into which numbers can be put by addition; same 
as A register. 

Action Cycle 

The type of memory cycle in which the memory is connected to 
the arithemtic unit for the transfer of a number. 

Addend 2-5 

The number in the memory to be added to the accumulator. 


Address 3-1 

The right-hand ten binary digits of a twenty-digit order; 
usually the location of a number in the memory. 

Analog 

Indicates a continuous operation (as opposed to digital). 
Arithmetic Unit 2-5 

The part of the computer which is required to carry out 
arithmetic operations. 

AS 3-9 

The double-length register including A and Q but excluding 
the first digit q Q of the Q register. 

* Numbers refer to pages. 
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Augend 2-5 


The number in the accumulator to which another number is 
to be added. 


-B- 


Binary Switch 7-4 

A portion of a program which is used to cause two separate 
sequences of orders to be followed alternately. 

Bit 2-1 


One binary unit of information or one binary "digit." 
Blocking Order 8-6 


A transfer order inserted 
sequence of orders so that 


in a program to stop the "normal" 
some check can be made on the program. 


Blunder 8-1 


A mistake in programming, or more usually a mistake in pro¬ 
gramming which is most obvious even without running the program 

Boolean 


Pertaining to the symbolic logic of George Boole, the 
mathematician. , 


Bootstrap Start 5-13 

QOoSfaS 4ornS' Wi i h °1 y the aid 0f the orl £ inal order pair 
4 °° 00 ' makes Possible for the machine to continue 
inputting and storing under the control of orders brought into 
the machine from the tape. ® 


C- 


Carriage Return 9-3 


The paper tape code which will cause the 
Printer to go to the left-hand margin of 


typing carriage of 
the paper. 


a 
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Cathode Ray Tube Output 12-1 

The cathode ray tube unit which can be used to display results 
from the computer directly by output orders. 

Check Point 8-6 

A point in a routine where some check is to be made. 

Circular Right Shift 13*7 

A right shift of AQ in which use is made of the digits shifted 
into Q. 

Clear 2-4 

The operation of changing all the digits of a register to a 
common state, usually 0. 

Code 

A group of more or less arbitrary symbols used to represent 
some other group of things. 

Coding Error 8-1 

An error in the detailed preparation of a program. 

Collate 3-27 

Digitwise logical product, same as extract. 

Comparer 9-16 

A device for checking the identity of two tapes. 

Complement 2-2 

See one's complement and two's complement. 

Conditional Transfer 3“5 

An instruction which will cause a transfer from the pattern 
of taking orders sequentially if the sign digit of the 
accumulator is 0. 


Control 1-4 


That part of the machine which serves to "control" the arith¬ 
metic unit, the memory and the input-output. 


15-3 



Control Counter 8-4 


Usually the order counter. 

Control Transfer 3-5 

An instruction which may cause a transfer from the sequential 
pattern for handling orders; see unconditional control transfer 
and conditional control transfer. 


Counter 


Usually a short routine which counts the number of iterations 
carried out in some program. 


-D- 


Diagnostic Routine 

A routine for finding and diagnosing a machine fault. 

Digital 1-2 

Pertaining to digits; discrete units. 

Digitwise-complement 2-2 

Same as one's complement; a binary number is the digitwise 
complement of another binary number if and only if each digit 
of the number disagrees with the corresponding digit of the 
other number. 

Directive 5-4 


An order, usually on a tape, which specifies the location at 
which a subroutine is to be stored. 


D.O.I. 5-1 

The decimal order input routine. Library Routine XI; used as 
a general purpose input routine for almost all programs. 

Double Precision 

Pertaining to numbers of approximately two 40 binary digit 
number lengths; hence requiring two registers or locations 
for a single number. 



Drum 


A rotating cyclinder with a ferromagnetic coating used as a 
memory. 


Error 

Mathematically the difference between the correct result and 
the computed result; often used in place of blunder. 

Error Frequency 13-2 

The average number of machine errors per unit time. 

Even Order 3“1 

The left-hand order of an order pair, thus using digits 2 
through 2“^9. 


Extract Order 3~27 

An operation which puts a 1 In the Q register wherever the 
number In the designated memory location and the number In Q 
are both 1, leaving all other digits zero; also called a 
logical multiply order because the resulting number in the 
q register is the digitwise binary product of the Q register 
and the number in the memory; same as collate. 


-F- 


Fixed Address 5-2 

The numerical address remains unaltered when put into the 
memory; usually in connection with D.O.I. where a fixed 
address is followed by an F. 


Fixed Point 2-1 

The binary point is always in the same position of a register 
Illiac is a fixed point machine. 
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Floating Point 

Pertaining to operations in which numbers are represented 
by a number multiplied by a power of a base; thus numbers 
may have different multiplying factors; floating point 
operations in the Illiac must be obtained by programming. 

Fraction 

A number with an absolute value less than 1. 

Function Digits 3-2 

The first eight binary digits of an order; the T and V digits. 


-H- 


Hang up 

An unplanned stop of the machine due either to a machine 
fault or a coding error. 


-I- 


Input Routine 

A routine for inputting other routines, usually aiding in 
conversion from binary-decimal to binary. 

Instruction 3-2 

The operation in the machine designated by the first eight 
binary digits of an order. 


Integer 


Pertaining to integers; although Illiac is a fixed point machine 
with the binary point between the first two digits it is possible 
to carry out operations using integer numbers less than 2+^°. 

Interlude 5-11 

A routine which carries out operations and is then destroyed 
by overwriting as the problem continues, usually carried out 
during input of information. 
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Instruction Code 3-1 


The order code. 

Interpretive Routine 4-9 

A routine in which a sequence of operations (instead of a 
single operation) may he carried out on a number, the sequence 
designated by parameters carried with the number. 

Iteration 

A sequence of orders, usually to be executed more than once, 
and arranged to converge to some analytic result. 


-L- 


Leapfrog 2-17 

An engineering routine which may be used to test the machine. 


Left-hand Order 3-1 

The left-hand or even order of a pair of orders, thus using 
digits 2^ through 2”19. 

Line Feed 9-3 

The paper tape code which causes a Teletype to advance the 
paper one line; always in conjunction with carriage return. 


Link 4-4 

The part of the routine used to bring a subroutine into 
operation and designate the point of return. 

Location 1-3 

The designation of a number location, address or position in 
the memory. 

Logical Product 

When applied to numbers of more than one binary digit, the 
digitwise logical product; the same as collation of two 
numbers and the same as extract. 
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8-4 


A sequence of orders which may be carried out more than once 
automatically. 


-M- 


Machine Error 

An error caused by a fault in Illiac. A malfunction. 
Malfunction 13-2 

A mistake in a computer function. A machine error. 
Memory 1-2 


A device which stores numbers; 
or Williams memory in Illiac. 


usually the cathode ray tube 


-N- 


Number Register 2-5 

The temporary location in which the addend, subtrahend, 
multiplicand and divisor are automatically placed by the 
control while the corresponding order is being executed; R 3 . 


— 0 — 


Odd Order 3-1 

The right-hand order of an order pair consisting of the 
digits 2 through 2~39. 

One 1 s Complement 2-2 

Same as digitwise complement; each digit of a binary number 
is changed to get its binary complement. 
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Order 

A set of 20 binary or 5 sexadecimal digits which is used to 
define a machine operation. 

Orde r Code 3-1 

A set of sexadecimal characters used to describe the operations 
in the Illiac. 

Order Counter 8-4 

The counter which keeps a record of the location of the next 
order pair if no transfer of control is required. 

Order Pair 3-2 

Two orders which are stored together in one location of the 
memory. 

Order Register 3“^ 

The register R_ into which an order pair is transferred from 
the memory just before either of the orders can be obeyed; R^- 


Output 

A display by oscilloscope, paper tape or Teletype printer of 
the contents of some part of Illiac. 

Overflow Digits 2-3 

The digits which, as a result of a computation, would require 
digits to the left of 2° in a register. 


-P- 


Page Printer 9-9 

An automatic typewriter. 

Paper Tape 

The common ll/l6 inch width paper tape used for handling 
numbers in and out of Illiac. 
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Parameter 4-7 


A number which may be altered from one computation to another 
but is held fixed during a single run. 

Partial Substitution Order 3-17 

An order which allows the address digits only to be stored 
in the memory; a 42 or 46 order. 


Plant 4-4 


The act of putting some number into an order or routine; 
usually with reference to an address. 

Playback 

The number coming from the magnetic drum memory. 

Position 

Referring to one of the 2 to 2 ^ places of the memory or 
of a machine register. 

Post Mortem 8-8 

The checking of the routine in the memory after it has been 
used, usually by an automatic routine which compares the 
contents of the memory with the original input tape. 

Program 1-4 

The plan of a calculation. 


Punch 9-1 


A unit to punch holes in ll/l6 inch paper tape, usually the 
punch providing output from Illiac. 


-Q- 


B or ^ Register 2-5 

The register which must be used for the multiplier during a 
multiplication; also the register into which the quotient is 
placed during division; the contents of the Q register. 
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Quotient Register 2-5 


Same as Q register. 

Quadrature 

The process of evaluating a definite integral by numerical means. 


-R- 


h 


3 -^ 

The order register; bottom row of 40 neon lights on Illiac. 


Range 

The set of numbers which may be handled in an Illiac register 
directly, from and including -1, up to and including 1-2 


Raster 


The array on a memory tube. 

Read-around 

An index of the interference of one spot with another in 
the Williams memory. 


Reader 9-1 

The device to take data from a punched paper tape for insertion 
into Illiac. 


Record 

The process of storing a number on the magnetic drum. 
Regeneration 

The process of refreshing the stored information on a 
Williams tube. 

Register 2-5 

A row of ^0 flipflops which can hold a binary number. 
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Relative Address 5-2 


An address which is relative to a directive and hence must 
be added to the directive address to get the true memory 
location. J 

Remainder 2-12 


The quantity left from the dividend after a division process 
(without regard tc Illiac). 

Reperforator 9-11 

A standard Teletype machine for making paper tapes from 
electrical signals. 

Residue 2-13 

The quantity in A after a division order has been obeyed 
(in the Illiac). 3 

Right-hand Order 3-1 

The order using positions 2 -20 through 2 _39 ; same as odd order. 
Round-off 2-12, 13 

The process of adding 2 in multiplication and making q = l 
in division when the term is used with Illiac. 39 


Routine 4-1 


A completed sequence of orders in coded form. 

Run 13-2 

A period of time during which a program is being executed 
by the computer. 


-13- 


Scale 6-1 


The adjustment of a number to come within range of Illiac. 
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Sexadecimal 3 - 2 

A number system with base l6. 


Shift 

The process of moving a number to the right or left in a 
register. 

Sign Digit 2-2 

The first position in a register} 2 . 

Single Address Code 

A machine order code which has one address only with each 
order; Illiac has a single address code. 

Slave Tube 8-4 

A device for displaying the entire contents of one memory 
cathode ray tube. 


Store 

The memory, usually Williams memory. 


Subroutine 4-1 

A routine arranged according to a standard pattern so that 
it can be easily used as part of other routines-. 


Sum Check 13-6 

A check on computer reliability made by adding a set of 
quantities twice and comparing the two results. 


-T- 


Tape Code 9-1 

The hole patterns in a paper tape that are used to represent 
numbers and instructions. 
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Tape Comparer 9-16 


A device to compare two paper tapes. 

T Disit 3-16 

The first sexadecimal digit of an instruction. 

Temporary Storage 

Locations in the memory used briefly during a calculation 
and not assigned to any number to be retained for results. 

Terminating Symbol 

A symbol on the tape indicating the end of a section of 
code or tape. 

Two r s Complement 2-2 

The difference between 2 and the number whose two's complement 
is to be found; in the Illiac identical, because it is modulo 
2, to reversing each digit and adding 2“39. 


-U- 


Unconditional Transfer 3-5 

A transfer out of the ordinary sequential pattern of handling 
orders, regardless of the sign of A. 

Useful Time 13-2 

The average useful time per run is the probability that no 
error will occur multiplied by the length of the run. 


-V- 


V Digit 3-6 


The second sexadecimal digit of an instruction. 
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-w- 


Waste Order ^-3 


An order that serves no computational purpose but is 
inserted because of ease in coding or limitation of the 
machine because of odd and even pairing o or ers. 


Williams Tube 


ft rav tube used for storing binary digits 


Word 

Forty binary digits in a register or a single location- 
in the memory. 


Working Space 

See temporary storage. 
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CHAPTER 16 


THE DRUM STORAGE UNIT 

16.1 GENERAL DESCRIPTION . Auxiliary storage of 12,800 
words is provided for the Illiac by the magnetic drum storage unit. 
Some knowledge of the physical characteristics and mode of operation 
of the drum storage unit is necessary for its most efficient use. 

The magnetic drum is, physically, a rotating cylinder whose 
surface is coated with magnetic material. Lengthwise, the cylinder 
is subdivided into 200 tracks; associated with each track is a re¬ 
cording and reading head capable of recording digital information 
on the magnetic surface and capable of reading data previously 
recorded (Page 16-2). The periphery of the cylinder is subdivided 
into 64 sectors; as a sector passes under one head, the digits of 
one word are recorded or read. Thus, the selection of any one of the 
12,800 words requires a spatial selection of a track and a temporal 
selection of a sector; that is, access to one word requires selection 
of one of 200 heads and the selection of one of 64 time intervals 
during which the digits of the desired word pass under the selected 
head. 

16.2 COMMON DRUM INSTRUCTIONS . The record instruction 
86 11 OOp places the contents of A in drum location p and performs 
eleven left shifts of AQ. The playback (or read) instruction 

85 11 OOp performs eleven left shifts of Q and places the contents of 
drum location p in A. Either type of drum instruction must occupy a 
complete 40 digit word. 
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The address p of either drum instruction is specified by 

_O/ Q 

the rightmost 14 binary digits (2 to 2 ) of the 40 digit 

instruction and hence is obeyed modulo 2^ = 16,384. It is restricted 
to the range 0 < p < 12,800 for playback and 2560 < p < 12,800 for 

14 

record, with the interpretation in both cases being made modulo 2 
Of the 14 binary digits specifying p, the least significant six 
(2 - ^ to 2“ 2 ^) determine in which of sixty-four sectors the desired 
word lies; the remaining eight (2“ to 2~ JJ ) specify the track. 

Except in special circumstances, the programmer need not be concerned 
with the breakdown of the address p into track and sector addresses, 
since the track switching is automatic and imposes no special timing 
restrictions on drum use. 

The section of the drum 0 < p < 2560 is used for storage of 
the drum bootstrap and commonly used routines; no recording by drum 
instructions is possible in this range. If either the record or 
playback instruction has an address in the range 12,800 < p < 16,384 
(mod. 16,384), the instruction will cause the computer to stop. 

For most efficient use of the drum storage unit, the desired 
sector of the drum should be approaching the reading and recording head 
when a drum instruction is to be executed. Subroutines for transfer of 
blocks of words between the drum and electrostatic memories require time 
for execution of other Illiac instructions between successive drum 
accesses; to facilitate such use, sectors whose Illiac addresses differ 
by one are physically five sectors apart on the drum surface (Fig. 16.1). 
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The circuits for switching from one track to another are 
so designed that all words on the same sector are equally accessible. 

For this to be possible, a delay of one word time precedes any drum 
consultation. 

As an example, suppose that the digits of word 0 have been 
drum track 0, sector 0. As the drum rotates (Fig. 16.1), 
sectors 13, 26, and 39 pass under the reading heads. This interval 
Of three word times is available for calculation and for address 
modification. The drum instruction requesting consultation of 
sector 1 should be in the Illiac order register before sector 52 begins 
to pass under the reading heads. The enforced delay of one word time 
ccurs as sector 52 passes under the head, whereupon the digits of any 

word in sector 1 can be read into the Illiac accumulator as sector 1 
passes under the heads. 

16-3 SALem-ATloti OF ACCESS TIM,. The timing data necessary 
for efficient use of the drum storage unit can be deduced from the 
following facts: 

1. The time required for one revolution of the drum is 16.9 milliseconds. 

2. Since the digits of 64 words are stored serially around the periphery 

of each track, the time required for the digits of one word (one 

sector interval) to pass a reading or recording head is 0.264 
milliseconds* 

3. The drum instruction must be in the Illiac order register at least 
one sector interval (0.264 milliseconds) before the digits of the 
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word can be transferred between the drum and the Illiac accumulator. 
Thus, the minimum time required for execution of a drum instruction 
is two sector intervals or 0.53 milliseconds. 

4# words corresponding to successive addresses in drum instructions 
are spaced five sectors apart on the drum surface. 

An isolated drum instruction with no definite timing relation¬ 
ship to other drum instructions may be regarded as occurring when the 
drum is in a random position. In this case, the time taken by the drum 
instruction will be equally likely to lie anywhere in the range between 
the minimum of 0.53 milliseconds and the maximum of one revolution time 
plus one sector interval or 17*43 milliseconds. The average time would 
thus be 9*0 milliseconds. 

A program does not usually have isolated drum instructions, 
but instead it will contain a sequence of drum instructions separated 
by other Illiac instructions. Under these circumstances, three timing 
calculations are of importance. 

1. The time t between the completion of one drum instruction and the 
completion of the next drum instruction. 

2. The time T required by all the Illiac instructions which come 
between the two drum instructions. 

3. The maximum time T available for execution of Illiac instructions 

^ max 

if the second drum instruction is to be obeyed as quickly as 
possible after the first. 

The commonest case is that in which successive drum instructions 

have successive addresses. The time T is then equal to three sector 

max 
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intervals or 0.79 milliseconds. If T < T then the time t is 

max * 

equal to five sector intervals, or 1.32 milliseconds. If the time T 

required for execution of Illiac instructions exceeds T , at least 

max 

one additional drum revolution will be required for each access, 
giving t = 18.2 milliseconds as the time between completion of drum 
instructions. 

A general calculation of t and for arbitrary addresses 

of two successively obeyed drum instructions can be made in the following 

way. If the address of one drum instruction is and that of the next 

is P 2 , °ne first calculates 5 (p 2 ~ P-^) - 2 — 3*8T and forms its residue 

modulo 64. (All times are in milliseconds.) Call this quantity R. In 

other words, R is formed by successively adding or subtracting 64 from 

the above expression so as to make R lie in the range 0 < R < 64 . The 

time t between completion of the two drum instructions is then given by 

the formula t = 0.53 ♦ T + .264R. The time T is chosen as the 

max 

shortest time T to make R = 0. Thus one obtains 

T max = * 26/f {Residue of [5 (p 2 - p^) - 2 ] (modulo 64 )} 

16.4 CALCULATION OF TIME OF EXECUTION OF ILLIAC INSTRUCTIONS . 
If access time to the drum storage is to be minimized, it is important 
that accurate calculations be made of the time of execution of Illiac 
instructions obeyed between successive drum instructions. The most 
accurate way of performing this calculation is to determine the number 
of clock periods of the Illiac electrostatic storage unit. The clock 
periods required for instructions commonly used in drum transfer routines 
are given in Table 16.4. 
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Readout of each order pair 2 

Order types K, S, F, and L 4 

Order types 2, 4, 5, J, and 3 when 

control is transferred. 2 

Order types 85, 86, and 3 when control 

is not transferred. 0 


Table 16.4 

Two clock periods must be allotted for each order pair, 
including the drum instruction itself and regardless of whether 
one or two orders are obeyed. The remaining clock periods listed 
in Table 16.4 are in addition to those allotted for order pair 
readout. The clock period of the Illiac storage unit is 
0.0185 + 0.0002 milliseconds; the maximum time T required for 
executing instructions between successive drum orders is found by 
multiplying the total number of clock periods required for the 
instructions by 0.0187. As an example, the time required for the 
record loop in library routine Y1 is calculated below: 
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Location 


Order 


10 

40 

L5 

11 L 
(n) F 

11 

86 

11 

F 


00 


F 

12 





F5 

11 

L 


40 

11 

L 

13 





F5 

10 

L 


42 

10 

L 


14 

FO 35 L 
32 10 L 


2 readout order pair 
4 order type L 


2 readout order pair 
0 order type 86 


2 readout order pair 
4 order type F 
2 order type 4 


2 readout order pair 
4 order type F 
2 order type 4 


2 readout order pair 
4 order type F 
2 control transfer 

32 Total 


The time T required for the record loop is therefore 32 x 0.0187 - 0.598 
milliseconds, which is less than the time of three sector intervals 
(0.792 milliseconds). Successive words will therefore be recorded at 
the rate of 1.32 milliseconds per word, after the first word has 
been recorded. 

16.5 OTHER TYPES OF DRUM INSTRUCTIONS . A more general 
type of drum instruction my be written in the form 8 11 T V ^ P* 

If T is 0, 1, 8, or 9 the drum instmction will be obeyed as described 
in the first section. If T is any other digit type the right hand 
half of the instruction will be obeyed as a separate Illiac instruction 
after the drum instruction has been obeyed. In such a case the drum 
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address is interpreted modulo 8,192 - 2 13 while the drum instruction 
is being obeyed and for this reason no hangup is possible. The 
rightmost binary digit of V 2 is uniquely deterained by the drum address 

p but V 0 is otherwise arbitrary. 

Possible choices for Vj. include 5, 6, 7, J and L. The 
playback instructions have ^ equal to 5 or J and both produce the 
same result. It is customary to take ^ - 5 for playback. The 
record instructions have V x equal to 6, 7, and L, and they record H(A), 
zero, and l/2 respectively as may be deduced from the characteristics 

of the variant V^. 

If V, is taken equal to h a playback is executed but since A 
is not cleared before carrying out the playback the result wall also 
depend on the initial contents of A. The result of this instruction is 
to shift A and Q left eleven places but during the shift to insert 
zeros into right hand eleven digits of A rather than digits from Q. 

A digitwise inclusive or (or logical sum) of this shifted quantity in 

A is formed with the word taken from the drum. 

A number of other instructions exist which affect the drum 

but which also have little utility. They include the cases: 

a. Drum instructions with the left hand address different from 11. 

b. Drum instructions with the left hand type digit equal to 9 
rather than 8. 

c. Drum instructions lying entirely within the right half of a 
word. (These instructions will produce the conventional effect 
only if the address equals eleven mod 6A.) 
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Since applications of these types of instructions to 
actual programs are so limited we give no description of their 
properties here* 

16 • 6 MODIFICATION OF ILLIAC . The initial instruction 
pair 8002840000 which has been used to start the tape bootstrap was 
rewired so that the instruction pair 8500S40000 is used instead. 

This change will have no effect in the operation of the computer 
except when a bootstrap start is used, and even then the computer will 
usually respond in the same way now as it did with the original 
instruction pair. The differences may be noted as follows: 

1. No D.O.I. is now required on the front of your tape since if the 
D.O.I. is omitted it is taken from the drum automatically and 
placed in locations 3F7 to 3LL. The additional locations 3F5, 

3F6, 000, and 001 are used in this process. 

2. If a tape bootstrap such as that appearing on the D.O.I. tape is 
at the beginning of your tape, the computer will handle it in 
the same way as it has in the past except that locations 3L6 
through 3LL will be used before the bootstrap is read. 

3. If your program requires the reading of a jump instruction by means 
of a bootstrap start, it must be altered. A jump instruction 
written as 2406400000 must be replaced by 00K24100N, while if the 
instruction 263F700000 is to be read on a bootstrap start this 
instruction may be omitted entirely. In any such case, the 
locations 3F5, 3F6, 000, 001 will be used and the D.O.I. will be 
replaced in the Williams memory. 
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The program which provides the facilities described above 
consists of a set of routines which are stored on the drum. The 
initial instruction pair 8500S40000 plays back a drum bootstrap which 
occupies locations 000, 001, 400, 401 on the drum and 000, 001 in the 
high speed memory. It, in turn, plays back routine 1 from 3L6-3L9 
on the drum to 3L6-3L9 in the high speed memory. Routine 1 plays 
back routine 2 from 3LK-3LL on the drum to 3LK-3LL in the high speed 
memory. Routine 2 records the contents of the high speed memory 
locations 002-01L at locations 31F2-31LL on the drum to preserve them 
and returns control to routine 1. Routine 1 plays back routines 3 
and 4 from locations 002-013 on the drum to 002-013 in the high speed 
memory an* jumps to routine 3. Routine 3 reads one sexadecimal 
character from the tape to determine which routine should be played 
back next and prepares routine 4 accordingly. If the character 0 is 
read from the tape, routine 4 will play back part of the D.O.I. and 
routine 5 and will jump to routine 5. Routine 5 replaces the words 
in 31F2-31LL in locations 002-01L in the high speed memory and then 
plays back the rest of the D.O.I. Entry is made to the D.O.I. in 
such a way that the reading of the initial function digit 0 of the 
directive is suppressed the first time. This is necessary since 
routine 3 has already read the initial 0. Routine 5 occupies locations 
3P5 , 3 f 6, 3L9-001. All except locations 3F5, 3F6, 000, and 001 are 

covered by the D.O.I. 
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If routine 3 reads the character 8 from the tape a tape 
bootstrap is indicated and routine k is prepared so that it plays 
back routine 6 into locations 3L8-3LL and jumps to routine 6. All 
words of routine 6 are contained in routine 5, but entry is made at 
a different point. Routine 6 replaces the words in 31F2-31LL in 
locations 002-01L and causes the next nine characters to be read 
from the tape and stored in 000 after faking the initial 8 which was 
already read. Routine 6 then jumps to 000 and the usual tape boot¬ 
strap sequence has been started. 

Other initial characters have been assigned as follows: 

(1) Leapfrog 

(3) Flying Leapfrog 

(F) Drum Post Mortem Routines 

(L) Store Post Mortem Routines 

USE OF T HE DRUM . The drum will most commonly be 
used by means of library subroutines. When this is done no special 
knowledge of drum instructions is required, and all problems of 
reducing access time to the "minimum" of 1.32 milliseconds are handled 
automatically by the subroutine. 

Y-l is a typical subroutine for transferring information 
back and forth from the drum. It may be used to record a block of r 
words on successive locations of the drum starting at location q 
provided these words are present in successive locations in the 
Williams memory starting at location n. One must enter the subroutine 
by means of the instructions 
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p 


JO n 
50 p 


p+1 26 — 

00 q 

p+2 00 r 

in order to cause the recording to occur. Control will be returned 

to the right hand side of p+2 after the recording. 

The same subroutine may be used to play back r words from 
successive locations of the drum starting at location q provided these 
words were previously recorded by the subroutine. They wall be 
stored in successive locations of the Williams memory starting at any 
arbitrary location n. In order to perform this playback and store 
operation one enters the subroutine by means of the instructions 

50 n 

p 50 P 

p+1 26 — 

00 q 

p+2 00 r. 

In addition to recording and playing back a block of words 
in minimum time, this subroutine computes a sum check for the block 
of words, which is stored on location q + r of the drum. Upon play¬ 
back this sum check is recomputed and compared with the recorded sum 
check. If any word has been recorded or played back incorrectly, the 
two sum checks will not agree and the computer will stop on the 
instruction FF010* 

In some complex problems one may wish to play back words 
which do not occupy successive locations on the drum, and if this is 
the case library subroutines are not practicable. For example, it 
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may be desirable to record a matrix in such a way that either a 
row or a column can be played back with minimum access time. If 
the matrix has no more than 65 columns, one may achieve minimum 
access for either rows or columns by recording the elements of each 
row m successive locations but recording successive rows 65 locations 
apart. Thus, we see that T mx (as described in 16.3) is 1.32 milliseconds 
for both rows and columns and corresponds to the usual "minimum" access 
time. If the number of columns N lies in the range 65< N < 129 one 
must record the rows 129 locations apart, etc. 

If a technique of this sort is to be used with the drum, it is 
up to the programmer to write his own record and playback loops with 
due attention to access time and to sum checking. 
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INDEX 


-A- 


A register, 2-5 
Accumulator, 2-5 

shifting without shifting Q, 7-1° 
Accuracy, 6-7 
Add from Q, 2-8, 3-26 
Addend, 2-5 

Addition, 2-5* 2-7, 3-29 
Address, 1-3* 3-1 

digits, 3-8 
fixed, 5-1 

relative, 5-1 

Aft, 3-9 

Arithmetic checking, 13-15 
Arithmetic unit, 1-2, 2-5 
Augend, 2-5 


■B- 


Binary 

chopping, 7-15 
-coded decimal, 9-3 
digit, 2-1 
point, 2-1 
point shifted, 6-1 
switch, 7-4 
Blocking order, 8-6 

example of use, 11-37 
Blunders, 8-1 

list of typical, 8-3 
discovery of, 8-2 
Bootstrap, 5-12 

start, 5-13 


-C- 


Camera operation, 12-4 
Cathode ray tube display, 12-1 

orders for controlling, 12-3, 4 
programming for, 12-5 
example of result, 12-9 



Check, 

programmed, 13-1 
sum, 13-6 
Checking, 8-1 

arithmetic, 13-15 
drum, 13-14 
input, 13-6 
memory, 13-12 
methods, table of, 13-17 
of a particular program, II -36 
Check point, 8-6 
Coding tricks, 7-1 
Combination Teletype, 9-11 
Comparing of tapes, 9-16 
Complement, 2-2 
Conditional transfer, 3-5 » l6 
Constants, orders and addresses as, 7-7 
Control counter, 8-k 
Control transfer, 3-5 
reversed, 7-3 
Control unit, 1-4 
Convergence criteria, 7-12 
Correction tape, 8-5 
Cyclotron trajectories, 11-1 
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Data tape, 13-9 
Decimal order input, 5-1 
Definition of terms, 15-1 
Differential equations, 11-17 

auxiliary subroutine, 11-20 
Digital, 1-2 

Digitwise complement, 2-2 
Directive, 5-^ 

Division, 2-12, 3-19 

"double precision", . 2-18 
hangup, 2-16 
remainder, 2 - 15 , 7-1^ 
special cases, 2-16 

Drum, 

checking, 13-15 
playback, 3-21, 13-1^ 
record, 3-22, 13 - 1 ^ 
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Error frequency, 13-2 
Extract, 3-27 

-F- 


FF order, 3-28, 13-8 
Fifth hole, 9*^ 

Fixed point, 2-1 
Format characters, 9-3 
Fractions, input of, 5-8 


-G- 


Glossary, 15-1 


-H- 


Hand punch, 9-17 


-I- 


Increment add, 2-8, 3-28 

from Q, 2-8, 3-25 

Information sources, 8-4 

Input, 9-1 , 

Input checking, 13-8 
Input-output orders, 3-21, 29 
Institute for Advanced Study, 1-1 
Instruction, 3-1 
Integer operations, 2-20 
Interchange of numbers, 7-H 
Interlude, 5-l 1 
Interpretive routine, 4-9 



-K- 


Keybcar^ perforator, 9-9 


-L- 


Leapfrog, 2-17 

Left shift, 2-4, 2-9, 3-12, 7-l8 
Letter printing, 9-4 
Library, l4-l 
Link, 4-4, 

Location, 1-3 
Logical product, 3-27 
Loop, 8-4 


-M- 


Machine error, 13-2 

Magnetic drum memory, 16-1 

Malfunction, 13-2 

Marking, 7-13 

Master tape, 13-9 

Mean time between errors, 13-2 

Memory, 1-2 

checking, 13-12 
to Q, 2-17, 3-18 
Modifications to program, 8-4 
Multiplications, 2-11, 3-20 
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Negative numbers, 2-2 
Number register, 2-5 


- 0 - 


Qne-address code, 3-1 



Order, 1-4 

as a constant, 7-7 
blocking, 8-6 
code, 3-1 

execution time, 10-2 
pair, 3-2 
types, 3-6, 7 
variants, 3-8 
Ordvac, 1-1 
Output, 9-1 
Overflow, 2-3 
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Fage printer, 9-9 
Parameters, 4-7 

preset, 4-11, 5-7 
program, 4-7 
use of, 11-16 

Photo-electric tape reader, 9-17 
Playback from drum, 13-14 
Poisson distribution, 13-3 
Polynomial evaluation, 7-16 
Post mortem routine, 8-8 
Product summation, 7-1 
Program, 1-4 

complete, 11-1 
modification, 8-4 
tape, 13-9 

Programmed check, 13-1 


-Q- 


Q register, 2-3 
Quotient register, 2-5 

for interchanges, 7-H 



Ry 3-4 
R 3 , 2-5 



Record on drum, 13-14 
Register, 2-5 
Remainder, 2-12 

in integer division, 7-14 
Reperforator, 9-11 
Resetting of cycles, 7-8 
Residue, 2-13 

Right shift, 2-4, 2-9, 3-13 
Routine, 1-4 

address search, 8-11 
check point, 8 - 8 , 12 
closed, 4-2 
entry to, 4-3 
exit from, 4-3 
interpretive, 4-9 
post mortem, 8 - 8 , 9 
sequence checking, 8 - 8 , 11 
Run, 13-2 

Running time, 10-1 

simple formula, 10-4 


-S- 


Scaling, 6-1 

adjustable, 6-10 
continuous, 6-10 
of a complete problem, 11-10 
Sexadecimal, 3-2 
Sign digit, 2-2 
Starting a program, 5-6 
Stop orders, 3-6 
Store, 2-17, 3-17 
Subtraction, 2 - 6 , 2 - 7 , 3-29 
Sum check, 13-6 


-T- 


Tape, 

data, 13-9 
master, 13-9 
program, 13-9 



Tape code, 9-1 

complete, 9-5 
orders for, 9-6 
sexadecimal, 9-1 
Tape preparation, 9-1 
Teletype equipment, 9-7 
Testing for 0 and -1, 7-6 
for 1/2, 7-H 

Transfer of control, 1-5# 5-5 
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Unconditional transfer, 3-5, 3-1^ 



UNIVERSITY OF ILLINOIS 
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ILLIAC ORDER CODE SUMMARY 

Input output and drum orders are not included. 
N = Storage location specified in order. 


Add, Subtract or Copy into A 75 microseconds 


; J ' L->' <-4 
^ -J i 

- Ht+<\ 


Inst 

Replace A by 

IA- 

A 

4- 

N 

LO 

A 

Ml 

N 

Ii6 

A 

4 

1*1 

L2 

A 


|N| 

f4 

A 

+ 

(N+2~ 59 ) 

FO 
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•m 

(N+2’ 39 ) 

f6 

A 

4* 

|N+2* 39 | 

F2 

A 

M 

|N+2** 39 j 


Inst Replace A by 



last Replace A by 

LJ 1/2 + N 
L9 1/2 - N 

LL 1/2 + |N| 

LS 1/2 - |N| 

FJ 1/2 + (N+2“ 39 ) 
F9 1/2 - (N+2 -39 ) 

FL 1/2 + |N+2~ 39 | 
FS 1/2 * |N+2~ 59 I 


ah 

A ■+ 

Q 

S 5 

Q 

SO 

A - 

Q 

si 

•'Q 

s 6 

A + 

|Q| 

S7 

W 

S2 

A * 

Kl 

S3 

-mi 


SJ 1/2 + Q 
S9 l/2 - Q 

SL 1/2 + |Q j 
SS 1/2 * |Q| 


Kb A + (Q+2" 59 ) 

RO A - (Q+2 -39 ) 

K6 A + |Q+2~ 39 | 

K2 A - |Q+2“ 39 | 


N5 


(Q+ 2 - ^ 9 ) 


KL «(Q+ 2 ” 59 ) 

K 7 |Q+ 2 ~ 59 | 
K3 »lQ+2“ 59 | 


KJ l/2 + (Q+2* 39 ) 
I© 1/2 - (Q+2*“ 39 ) 

KL 1/2 + lQ+2 - " 39 1 
KB 1/2 . |Q+2“ 59 | 
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0 y ^ 

j) 
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,/ o d r J 
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Transfer from A 


55 microseconds 


f a 

a 


. 6 


Inst Effect 


(4°) N = A 
46 1 .h «a. N = A 

42 r.h.a. N = A 


Inst Effect 


41 

(45) 

47 

45 


N = A = 0 

l.h.a. N = 0,1 
A = 0 

r.h.a. N = 0,] 
A = 0 


Inst Effect 

& ■ - A - 1/2 

4L l.h.a. N = 0, 
A = 1/2 

4S r.h.a. N = 0, 
A = l/2 


Copy into Q 


Instruction 

50' (52,54,56) 
51 (53,55,57) 
59 (5S,5J,5L) 

JO (J2,j4,j6) 
J1 (J3,J5,J7) 

J9 (JS,JJ,JL) 


55 microseconds 

Effect 

Copy N, into Q 

« *» H ” and clear A 

». m « »i an( j pu t l/2 in A 

Logical product N and Q to O' 


If 

If 


If 

TT 


" " ,r " tr an( ^ dear A 

" " w ” " anri put l/2 in A 


Multiply setting 


= 0 


700 microseconds 


Inst 

HOLD 

Replace AQ by 

Inst 

unrounded 

Replace AQ by 

Inst 

74 

Q.N .+ 2 Jy k 

75 

Q *N * 

7J 

70 

-39. 

-Q.N + 2 A 

7i 

-Q.N. 

79 

76 

q|n|+ 2" 59 A 

77 

Q|N| 

7L 

72 

-q|n|+ 2“ 59 A 

73 

-Q|Nl 

7s 


ROUNDED 
Replace AQ by 

Q.N. + 2"^° 
-Q.N. + 2" 40 
Q|N| + 2~ k0 
—Q |N I + 2^° 


- 2 - 



Divide leaving residue in A 


800 microseconds 


Inst 

66 

6j 

■ (6s) 


Replace Q by 

(aq)/n 

2 39 Q/N if Q > 0} 2 -39 (1+Q)/N if Q < 0 

[l/2 + 2. 39 Q1 /n if Q > 0; [l/2 + 2~ 39 (l+Q)] /n if Q < 0 


Shift Aq Address = Number of shifts (mod 6k) i 0 
l6n microseconds 


Instructions 

00 

(02,04,06) 

01 

(03,05,07) 

09 

(OS,OJ,OL) 


Effect 


Shift left 
Clear A 

and shift left 

Put 1/2 in A 
and shift left 


Instruction Effect 


10 

(12.14.16) 
11 

(13.15.17) 

, 19 

(1S,1J,1L) 


Shift right 
Clear A 

find shift right 

Put 1/2 in A 
and shift right 


Transfer Control UNCONDITIONAL TRANSFER 55 microseconds 


Inst 

Effect 

26 

Jump to l.h. 

22 

Jump to r.h. 

24 

Stopj to l.h 

20 

Stopj to r.h 


CONDITIONAL TRANSFER 

Inst 

36 
32 
3b 
30 


Inst 

Effect 

27 

(37) 

Jump to l.h. 
A = 0 

23 

(33) 

Jump to r.h. 
A = 0 

iaIa 

OJ fA 

Stopj to l.h. 
A = 0 

21 

(31) 

Stopj to r.h. 
A = 0 

55 

microseconds 


Effect 

Jump to l.h. if A > 0 
Jump to. r.h. if A > 0 
Stop* to l.h. if A > 0 
Stop; to r.h. if A > 0 


lust Effect 


2L Jump to l.h., 

(3L) A = 1/2 


2S 

m 

2J 


(3J) 


Jump to r.h,, 
A = 1/2 

Stopj to l.h., 
A = 1/2 


29 Stopj to r.h., 
(39) A = 1/2 

(15 if no transfer) 


Final Stop 
OF 


S.' bill, 2/26/54 

rt/ 9/26/56 
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